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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development  Center 
(AEDC),  Air  Force  Systems  Command  (AFSC),  under  AEDC  Project  Number  CP52PW 
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produced  by  Calspan  Corporation/ AEDC  Operations,  an  Air  Force  contractor  operating 
the  aerospace  flight  dynamics  testing  program  at  AEDC,  Arnold  Air  Force  Base,  TN  37389. 
The  SKE  Project  Manager  was  Capt.  D.  Smith,  and  the  AEDC/DOFA  Project  Manager 
was  Mr.  M.  O.  King,  D.  R.  Carlson  was  the  Calspan  Project  Manager.  The  work  was 
conducted  during  the  period  March  1,  1990  through  June  30,  1991,  and  the  manuscript  was 
submitted  for  publication  September  26,  1991. 

The  authors  acknowledge  the  contributions  of  Dr.  J.  C.  Erickson,  Jr.,  for  his  editorial 
assistance  during  the  preparation  of  this  report.  Also,  the  authors  acknowledge  the  assistance 
provided  during  the  design  and  testing  of  the  the  code  from  W.  E.  Dietz,  T.  L.  Donegan, 
M.  J.  Rist,  Dr.  J.  L.  Jacocks,  and  J.  K.  Jordan. 
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1.0  INTRODUCTION 

Computational  fluid  dynamics  (CFD)  calculations  to  support  ground  testing  are  an 
important  aspect  of  the  AEDC  mission.  Aerodynamic  problems  related  to  wind  tunnel  testing 
present  two  important  challenges  to  CFD:  the  need  to  produce  timely  solutions  in  response 
to  testing  requirements,  and  the  solution  of  problems  related  to  complex  three-dimensional 
configurations.  The  main  obstacle  to  overcoming  these  challenges  is  the  difficulty  of  generating 
computational  meshes.  Although  the  time  required  to  generate  meshes  is  decreasing  with 
the  development  of  interactive  grid  generation  software  (Ref.  1),  it  is  still  desirable  to  keep 
the  region  that  requires  a  mesh  to  a  simple  topology  that  is  easy  to  modify.  In  addition, 
a  single  mesh  that  is  fine  enough  to  resolve  the  desired  aerodynamic  features  of  a  complex 
flow  field  may  be  too  large  for  available  computer  memory. 

To  address  these  problems,  Benek,  et  al.  (Refs.  2-5}  have  developed  a  method  of  domain 
decomposition  called  “chimera”  which  allows  a  system  of  relatively  simple  grids,  each 
describing  a  component  of  a  complex  aerodynamic  configuration,  to  be  combined  into  a 
composite  grid  to  yield  solutions  for  complex  flow  fields.  The  chimera  scheme  has  greatly 
increased  the  range  of  aerodynamic  problems  that  can  be  addressed  by  CFD,  and  has  been 
applied  to  store  separation  (Refs.  6-8),  cavity  flow  (Ref.  9),  transonic  tunnel  wall  interference 
calculations  (Refs.  10  -  1 1),  space  shuttle  ascent  (Ref.  12),  dynamic  body  movement  (Ref. 
13),  and  flow  in  an  artificial  heart  (Ref.  14). 

•  The  general  concept  behind  chimera  is  illustrated  in  Fig.  1 ,  which  depicts  two  independently 
generated  meshes  representing  a  flapped  airfoil.  The  flap  mesh  is  embedded  within  the  airfoil 
mesh.  Clearly,  the  flap  mesh  outer  boundary  can  receive  flow-field  information  interpolated 
from  appropriate  mesh  points  (often  referred  to  as  interpolation  stencils)  of  the  airfoil  mesh. 
However,  a  reverse  process  must  occur  as  well;  the  airfoil  mesh  must  receive  flow-field 
information  from  the  flap  mesh.  An  artificial  boundary  must  be  defined  within  the  airfoil 
mesh,  since  certain  points  within  the  airfoil  mesh  are  interior  to  the  flap  and  thus  do  not 
lie  within  the  airfoil/flap  flow  field.  The  artificial  boundary  points  of  the  airfoil  mesh  that 
are  fully  contained  within  the  computational  region  of  the  flap  mesh  can  be  updated  by 
interpolation  from  the  appropriate  mesh  points  of  the  flap  mesh.  Generally,  any  mesh  can 
receive  information  from  other  meshes  through  outer  boundary  and  artificial  boundary  points. 

The  interpolation  process  is  further  illustrated  in  Fig.  2,  which  depicts  a  portion  of  the 
overlap  region  between  the  airfoil  and  flap  meshes.  Airfoil  mesh  points  that  are  contained 
within  a  certain  region  surrounding  the  flap  are  excluded  from  the  computational  domain 
of  the  airfoil  mesh.  (In  chimera  terminology,  they  are  "blanked”  points,  or  hole  points.) 
The  exclusion  of  points  is  accomplished  by  defining  a  hole  creation  boundary  within  the  flap 
mesh  that  will  define  the  region  within  which  all  airfoil  points  are  to  be  blanked.  The  points 
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in  the  airfoil  mesh  surrounding  the  blanked  points  are  hole  boundary  points,  and  they  receive 
flow-fleld  information  interpolated  from  mesh  points  within  the  flap  mesh.  Correspondingly, 
points  on  the  outer  boundary  of  the  flap  mesh  receive  flow-field  information  interpolated 
from  mesh  points  within  the  airfoil  mesh. 

Application  of  the  chimera  scheme  requires  two  steps:  (1)  a  description  of  how  each  mesh 
is  to  communicate  flow-lleld  information  to  other  meshes,  and  (2)  execution  of  a  flow  solver 
that  uses  the  communication  information  generated  in  step  1 .  Step  1  is  performed  by  the 
code  named  PEGSUS.  The  processes  accomplished  by  PEGSUS  include  establishing  which 
boundary  points  in  a  mesh  must  be  updated  by  interpolated  flow  variables  from  other  meshes 
and  by  calculating  the  required  interpolation  coefficients  for  the  mesh  points  that  send 
information  to  an  interpolated  boundary  point  in  another  mesh.  The  relationship  of  PEGSUS 
to  the  chimera  domain  decomposition  scheme  is  depicted  schematically  in  Fig.  3.  The  individual 
meshes  and  user-defined  mesh  connection  data  are  input  into  PEGSUS.  PEGSUS  produces 
(1)  a  composite  mesh,  i.e.,  a  single  file  consisting  of  the  concatenation  of  all  meshes  in  the 
multiple-mesh  configuration,  and  (2)  an  interpolation  file,  which  is  a  table  associating  all 
interpolated  boundary  points  in  the  composite  mesh  with  mesh  points  that  supply  the 
interpolated  flow-field  values.  In  step  2,  the  composite  mesh  and  interpolation  file  are  input 
into  the  flow  solver.  The  flow  solver  can  be  any  finite-difference  or  finite-volume  code  that 
has  been  modified  for  the  chimera  interpolation  data  structure.  Several  compressible  flow 
codes,  namely  XMER3D  (Ref.  5),  F3D  (Ref.  12),  XLIM3D  (Ref.  15),  and  EAGLE  (Ref. 
16),  and  an  incompressible  flow  code  (Ref.  14)  can  currently  use  the  chimera  scheme. 

The  chimera  method  has  evolved  in  terms  of  its  capabilities  and  applicability  to  complex 
three-dimensional  flow  configurations.  Much  of  the  evolution  has  been  manifested  in  several 
versions  of  the  PEGSUS  program.  PEGSUS  1.x  required  the  definition  of  a  hierarchical  mesh 
structure  which  restricted  which  meshes  could  embed  solid  surfaces  in  other  meshes.  PEGSUS 
2.x  allowed  much  more  generality  in  relationships  among  meshes  by  eliminating  any 
hierarchical  relationships.  However,  PEGSUS  2.x  was  limited  in  the  generality  in  which  mesh- 
to-mesh  communications  could  be  defined,  was  restricted  to  simple  hole  boundary  shapes, 
and  required  different  mesh  topologies  to  be  treated  as  special  cases.  PEGSUS  3.x  (Refs. 
17  -  18)  added  several  enhanced  features  that  increased  the  flexibility  and  applicability  of 
the  code  to  multiple  mesh  problems.  Most  importantly,  PEGSUS  3.x  provided  greater 
flexibility  in  the  definition  of  hole  boundaries  and  mesh  interactions.  Complex  mesh 
interactions  and  hole  boundary  deHnitions  were  handled  through  input,  usually  without 
requiring  program  modification. 

PEGSUS  4.0  is  the  latest  version  of  the  PEGSUS  series  of  mesh  interpolation  codes. 
It  is  a  fully  three-dimensional  code.  PEGSUS  4.0  adds  further  flexibility  to  the  user  inputs 
and  control  of  the  overall  process.  Major  features  are  as  follows: 
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—  Meshes,  boundaries,  and/or  regions  can  be  added  and/or  subtracted  using  the 
restart  capability 

—  Single  or  double  level  hole  boundaries  can  be  specified 

—  Hole  boundary  points  on  a  reflection  plane  can  be  ignored 

—  The  restrictions  for  obtaining  a  valid  interpolation  stencil  can  be  relaxed 

—  Holes  can  be  defined  using  a  number  of  methods 

—  Regions  can  be  directly  blanked  for  interior  boundaries 

—  The  code  has  been  written  to  be  portable  (The  code  has  been  successfully  run  on 
Apollo  and  Silicon  Graphics  workstations  and  on  a  Cray  supercomputer). 

A  glossary  of  PEGSUS  terminology  is  provided  in  Appendix  A.  Appendixes  B  through 
F  contain  information  about  PEGSUS  program  variables,  program  structure,  subroutines, 
direct  access  flies,  and  output  file  structure,  respectively.  Appendixes  G  and  H  contain  complex 
aerodynamic  configuration  examples  and  sample  Cray  job  control  language  required  for  the 
examples,  respectively. 

The  purpose  of  PEGSUS  4.0  is  identical  to  that  of  earlier  versions,  namely  to  produce 
a  composite  mesh  and  interpolation  file  for  the  flow  solver  in  question.  To  produce  the  required 
inputs  for  the  flow  solver,  PEGSUS  4.0  has  four  basic  steps  to  execute.  These  steps  are; 
(1)  process  the  meshes  and  user  inputs,  (2)  identify  the  hole  and  interpolation  boundary  points, 
(3)  determine  the  interpolation  stencil  and  interpolation  coefficients  for  each  interpolation 
boundary  point,  and  (4)  supply  diagnostic  information  on  the  execution  and  output  the  results 
for  input  to  the  flow  solver. 

The  details  of  PEGSUS  4.0  are  given  in  the  subsequent  sections.  Section  2.0  describes 
the  overall  approach  adopted  by  PEGSUS  4.0.  The  methods  used  to  identify  holes  and 
interpolation  boundaries  are  described  in  Sec.  2.1,  whereas  the  intricacies  of  the  interpolation 
process  and  the  different  controls  that  can  be  placed  on  this  process  by  user  inputs  are  described 
in  Sec.  2.2.  Section  2.3  concludes  the  description  of  the  overall  approach  with  a  discussion 
of  the  additional  features  available  for  the  control  of  PEGSUS  results.  Section  3.0  describes 
the  implementation  of  PEGSUS  4.0.  The  program  parameters  are  described  in  Sec.  3.1.  Section 
3.2  contains  a  description  of  the  user  inputs  referred  to  earlier.  In  Sec.  3.3,  the  format  for 
the  input  meshes  is  given.  The  files  output  by  PEGSUS  4.0  are  described  in  Sec.  3.4.  Finally, 
six  examples  demonstrating  many  of  the  capabilities  of  PEGSUS  4.0  are  given  in  Sec.  4.0. 
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Since  the  descriptions  in  Secs.  2.0  and  3.0  are  given  in  relatively  general  and  abstract  terms, 
careful  study  of  the  examples  in  Sec.  4.0  and  Appendix  G  will  aid  the  reader  in  acquiring 
a  working  knowledge  of  the  concepts  and  their  implementation. 

2.0  GENERAL  APPROACH 

This  section  discusses  the  procedures  used  to  determine  which  points  in  a  mesh  will  require 
interpolation  and  how  data  for  these  points  are  interpolated.  Also  presented  are  the  special 
features  of  PEGSUS  4.0  that  allow  additional  control  of  the  results.  A  glossary  of  terms 
is  given  in  Appendix  A  to  assist  the  reader,  in  grasping  the  concepts  and  their  practical 
implementation. 

2.1  INTERPOLATION  BOUNDARY  POINT  SPECIFICATION 

There  are  two  ways,  hole  boundaries  and  outer  boundaries,  in  which  information  is 
transferred  from  one  mesh  to  another.  The  user  must  identify  these  interpolation  boundaries 
so  that  PEGSUS  can  find  the  interpolation  stencil  and  calculate  interpolation  coefEcients 
for  each  boundary  point.  The  identiHcation  of  interpolation  boundaries  is  made  through 
the  user  inputs  SBOUNDARY,  SSURFACE,  $BOX,  SREGION,  and  SVOLUME  (see  Sec. 

3.2  for  format  details).  These  user  inputs  are  supplied  in  a  NAMELIST-like  format  that  is 
read  by  an  interpreter  within  PEGSUS  and  thus  is  not  restricted  to  any  particular  FORTRAN 
compiler  (since  NAMELIST  is  not  standard  FORTRAN  77).  The  term  NAMELIST  that  is 
used  throughout  this  report  refers  to  this  modified  NAMELIST-like  format.  Using  these  inputs, 
three  methods  for  specifying  hole  boundaries  and  one  for  specifying  outer  boundaries  are 
defined. 

2.1.1  Blanking 

Before  considering  the  different  methods  for  identifying  interpolation  boundary  points, 
the  definition  of  blanking  must  be  understood.  The  identification  of  a  point  as  a  hole  point, 
an  interpolation  boundary  point,  or  a  field  point  must  be  passed  on  to  the  flow  solver  so 
that  the  point  will  be  handled  correctly.  The  status  of  each  point  is  contained  within  an  integer 
array  called  IBLANK.  The  IBLANK  array,  which  is  dimensioned  identically  to  the  number 
of  points  in  each  mesh,  has  a  value  for  each  point  within  each  computational  mesh.  The 
values  contained  within  the  IBLANK  array  are  either  0  for  a  point  that  is  a  hole  or  interpolation 
boundary  point,  or  1  for  a  point  that  is  a  field  point.  Within  the  flow  solver,  points  with 
IBLANK  =  0  are  either  considered  to  be  completely  out  of  the  computational  domain  or  are 
to  be  updated  by  interpolation,  while  points  with  IBLANK  =  1  are  updated  by  the  solution 
algorithm  and  associated  boundary  conditions.  The  IBLANK  array  is  supplied  to  the  flow 
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solver  through  the  interpolation  file  output  of  PEGSUS  (see  Appendix  F).  References  2  and 
4  describe  how  the  IBLANK  array  is  incorporated  into  the  flow  solution  algorithm. 

2.1.2  Hole  Boundary  Specification 

Three  methods  for  creating  holes  are  presented.  Two  of  these  methods  use  an  indirect 
means  (not  specifying  the  points  in  the  hole  directly),  while  the  third  method  directly  specifies 
the  points  that  make  up  the  hole.  The  first  method  uses  the  SSURFACE  and  SBOUNDARY 
user  inputs  to  indirectly  specify  hole  locations.  The  user  input  for  this  method  is  defined 
in  the  computational  domain  but  utilizes  surface  normals  in  the  physical  domain  to  determine 
the  hole  points.  The  second  method,  also  for  indirect  hole  specification,  uses  the  $60X  and 
SBOUNDARY  user  inputs.  These  input  sets  describe  the  hole  in  the  physical  domain  and 
employ  the  physical  domain  to  find  the  hole  points.  The  third  method,  which  is  a  direct  means 
for  defining  holes,  uses  the  SVOLUME  and  SREGION  user  inputs.  These  input  sets  define 
the  hole  in  the  computational  domain  and  are  applied  in  the  computational  domain  to  find 
hole  points.  The  two  indirect  methods  are  described  in  Sec.  2. 1.2.1  while  the  direct  method 
is  described  in  Sec.  2.1. 2^.2.  Once  the  holes  are  found,  a  hole  boundary  is  placed  around 
each  hole.  A  description  of  the  method  used  to  determine  hole  boundaries  is  given  in  Sec. 
2. 1.2.3.  Finally,  in  Sec.  2. 1.2.4  special  problems  which  occur  in  the  hole  creation  process 
are  presented  along  with  their  solutions. 

2.1.2. 1  Indirect  Hole  Specification 

Typically,  a  hole  is  defined  through  indirect  means,  since  the  actual  points  within  the 
hole  are  not  known  beforehand,  and  the  actual  hole  boundary  is  irregular  in  shape  (not  a 
continuous  J,  K,  or  L  plane,  where  the  grid  points  in  each  mesh  are  indexed  by  separable 
indices  J,  K,  and  L).  The  first  method  uses  the  SSURFACE  user  input  (identical  to  the  method 
in  PEGSUS  3.x)  and  requires  the  use  of  a  surface  or  group  of  surfaces  from  one  mesh  to 
create  a  hole  in  another  mesh.  An  example  of  this  method  is  given  in  Sec.  4. 1 . 1 .  A  hole  creation 
boundary  consists  of  a  surface  or  a  group  of  surfaces.  The  purpose  of  a  hole  creation  boundary 
is  to  identify  points  that  are  within  this  boundary.  A  mesh  point  is  considered  to  be  inside 
a  hole  creation  boundary  if  it  is  inside  all  surfaces  that  define  the  boundary.  The  method 
used  to  determine  whether  a  point  is  inside  or  outside  a  surface  is  illustrated  in  Fig.  4.  A 
mesh  point  (P)  is  considered  to  be  inside  a  surface  if  the  dot  product  between  R  (the  vector 
from  the  closest  point  on  the  surface  to  the  mesh  point),  and  N  (the  normal  vector  on  the 
surface  at  the  closest  point,  directed  outward  from  the  hole  region)  is  negative  or  zero.  If 
the  dot  product  is  positive,  the  mesh  point  is  considered  to  be  outside  the  surface. 

The  hole  location  process  for  a  triangular  hole  creation  boundary  consisting  of  the  three 
surfaces  SI,  S2,  and  S3  is  illustrated  in  Fig.  S.  PEGSUS  first  puts  the  indices  of  all  points 
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of  the  mesh  in  which  the  hole  is  to  be  created  into  a  list.  Then,  all  points  that  fall  outside 
surface  SI  are  eliminated  from  the  index  list  (as  indicated  by  the  shaded  region  in  Fig.  5a). 
Then  all  points  contained  in  the  shortened  list  that  are  outside  the  surface  S2  (Fig.  Sb)  are 
eliminated,  which  shortens  the  index  list  further.  When  the  process  is  repeated  for  surface 
S3  (Fig.  Sc),  the  points  remaining  in  the  list  are  the  hole  points  of  the  mesh  in  which  the 
hole  is  created. 

This  method  of  finding  hole  points  has  one  important  restriction,  which  is  illustrated  in 
Fig.  6.  In  Fig.  6a,  a  concave  hole  creation  boundary  is  defined.  The  candidate  point  P  will 
be  found  to  lie  outside  Surface  1 ,  and  therefore  outside  the  entire  hole  boundary.  As  a  result, 
the  point  will  not  be  identified  as  a  hole  point,  which  clearly  is  an  error.  Generally,  hole 
creation  boundaries  (as  viewed  from  the  outside)  must  be  convex  to  obtain  correct  results. 
However,  PEGSUS  4.0  allows  holes  to  be  defined  by  multiple  hole  creation  boundaries.  Any 
hole  creation  boundary  that  contains  concavities  can  be  divided  into  one  or  more  entirely 
convex  regions,  each  of  which  can  create  holes  in  one  or  more  meshes.  If  the  hole  creation 
boundary  is  redefined  as  two  separate  convex  boundaries,  as  in  Fig.  6b,  the  candidate  point 
will  be  inside  one  of  the  hole  creation  boundaries,  and  will  therefore  be  correctly  identified 
as  a  hole  point.  The  candidate  points  within  the  two  boundaries  will  all  be  labeled  as  hole 
points  and  will  therefore  merge  into  a  single  hole. 

Now  that  the  method  for  defining  a  hole  using  surfaces  has  been  described,  the  relationship 
of  a  surface  or  group  of  surfaces  to  a  hole  creation  boundary  and  how  this  hole  creation 
boundary  relates  to  meshes  must  be  described.  The  relationship  between  a  hole  creation 
boundary  and  a  mesh  is  particularly  important  since  the  hole  creation  boundary  must  have 
a  connection  to  the  mesh  that  contains  the  surfaces  deHning  the  hole  as  well  as  a  connection 
to  the  the  mesh  or  meshes  that  will  have  a  hole  made  in  them.  The  connections  between  hole 
creation  boundaries  and  meshes  are  depicted  in  Fig.  7.  Boundaries  are  defined  as  collections 
of  level  surfaces  (i.e.,  surfaces  that  share  a  common  J,  K,  or  L  value)  within  meshes.  Shown 
in  Fig.  7  are  H-Boundary  a  and  H-Boundary  b,  which  are  hole  creation  boundaries  defined 
by  Mesh  n  and  Mesh  n  +  2,  respectively.  In  NAMELIST  nomenclature,  H-Boundary  a 
"ISPARTOF”  (IS  PART  OF)  Mesh  n  and  H-Boundary  b  “ISPARTOF"  Mesh  n+2.  Hole 
creation  boundaries  also  have  relationships  to  other  meshes,  i.e.,  they  make  holes  in  other 
meshes  (see  Fig.  2).  This  relationship  is  denoted  as  “MHOLEIN”  (Makes  a  HOLE  IN)  in 
the  NAMELIST  input.  According  to  Fig.  7,  H-Boundary  a  makes  holes  in  both  Mesh  m 
and  Mesh  n-1  while  H-Boundary  b  makes  holes  in  Mesh  n-1,  Mesh  n,  and  Mesh  n  + 1 .  The 
maximum  number  of  meshes  that  a  boundary  can  make  a  hole  in  is  set  by  the  program 
parameter  MHDIM  (see  Sec.  3.1).  The  examples  in  Sec.  4.0  and  Appendix  G  will  clarify 
these  relationships. 
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Figure  8  illustrates  the  relationship  between  surfaces  and  hole  creation  boundaries.  Hole 
creation  boundaries  are  comprised  of  collections  of  surfaces,  and  therefore,  surfaces  have 
“ISPARTOF”  relationships  to  boundaries  in  much  the  same  way  as  boundaries  do  to  meshes. 
Hole  creation  boundaries  are  normally,  but  not  always,  composed  of  level  surfaces  from 
the  meshes  to  which  they  are  connected.  This  situation  is  graphically  depicted  in  Fig.  8  with 
Surface  1  being  the  exception.  As  was  shown  in  Fig.  7,  H-Boundary  a  is  part  of  (ISPARTOF) 
Mesh  n  and,  therefore,  the  two  surfaces  composing  H-Boundary  a,  namely  Surface  1  and 
Surface  q-2,  would  normally  be  connected  to  Mesh  n.  In  this  case,  though,  one  of  the 
connections  must  be  overridden  by  the  input  “MESH  NAME”  (see  Sec.  3.2).  Surface  1  has 
a  “MESH  NAME”  specified  for  Mesh  n-2  in  the  NAMELIST  input  (Sec.  3.2)  and  is  therefore 
composed  of  a  level  surface  in  Mesh  n-2.  That  is,  a  level  surface  from  another  mesh  is  required 
to  successfully  complete  the  hole  creation  boundary.  An  example  demonstrating  hole  creation 
using  surfaces  from  two  separate  meshes  is  given  in  Sec.  4.1.4.  The  remaining  surfaces  are 
connected  to  meshes  through  the  ISPARTOF  specification  for  each  boundary  (See  Fig.  7). 

So  far,  the  hole  creation  boundaries  that  have  been  demonstrated  completely  enclose  the 
hole  region.  The  fact  that  most  hole  creation  boundaries  are  completely  closed  has  been  used 
to  accelerate  the  process  of  finding  holes  (hole  point  determination  was  one  of  the  most 
computationally  intensive  processes  in  PEGSUS  3.0).  The  acceleration  is  obtained  by  creating 
a  parallelepiped  based  on  the  minimum  and  maximum  X,  Y,  and  Z  values  of  the  surfaces 
specified  to  create  the  hole.  The  minimum/maximum  parallelepiped  is  used  to  quickly 
determine  if  a  point  can  possibly  be  in  the  hole.  If  a  point  is  outside  this  parallelepiped  then 
it  is  immediately  removed  from  the  list  of  candidate  hole  points.  This  results  in  a  greatly 
reduced  list  of  candidate  hole  points.  Those  points  on  the  inside  of  the  minimum/maximum 
parallelepiped  must  still  be  checked  to  see  if  they  are  actually  in  the  hole  or  not.  As  discussed 
earlier,  the  final  determination  of  whether  a  point  is  in  a  hole  is  done  with  surfaces.  The 
minimum/maximum  ranges  of  the  parallelepiped  are  given  in  the  PEGSUS  output  for  each 
hole  boundary  description.  However,  under  certain  circumstances  when  the  hole  creation 
boundary  does  not  clearly  distinguish  between  points  inside  and  outside  the  hole,  such  as 
when  the  user  has  not  closed  the  hole  creation  boundary,  the  minimum/maximum 
parailelepiped  may  not  produce  the  hole  that  the  user  wanted.  To  allow  for  these  special 
cases,  the  user  input  “CLOSED”  can  be  specifled  to  be  equal  to  .FALSE,  (default  is  .TRUE.). 
Figure  9a  illustrates  the  effects  of  such  a  specification  with  the  shaded  area  being  the  area 
within  the  hole.  In  some  instances,  for  example  when  the  boundary  is  not  closed  but  intersects 
an  outer  plane  of  the  domain,  the  specification  of  “CLOSED"  equal  to  .TRUE,  or  .FALSE, 
gives  the  same  results  (although  specifying  .FALSE,  increases  the  computational  time).  This 
case  is  demonstrated  in  Fig.  9b.  The  “CLOSED”  user  input  is  used  only  as  part  of  the 
$BOUNDARY  user  input  when  a  hole  creation  boundary  is  defined  by  a  collection  of  surfaces. 
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The  second  method  for  defining  holes,  also  an  indirect  method,  uses  the  $BOX  user  input. 
The  $BOX  user  input  is  connected  to  a  boundary  and  can  be  used  in  conjunction  with  surfaces 
or  by  itself.  An  example  that  demonstrates  the  $BOX  user  input  is  given  in  Sec.  4.1.2.  A 
box  is  described  by  a  minimum-to-maximum  range  of  X,  Y,  and  Z  coordinates.  This  box 
is  used  to  create  holes  in  those  meshes  to  which  the  boundary  is  connected  by  the  MHOLEIN 
input.  For  those  meshes  in  question,  any  mesh  point  located  within  the  box  vhill  be  considered 
to  be  a  hole  point.  The  box-hole  creation  boundary  connection  is  depicted  in  Fig.  10.  H- 
Boundary  a  has  three  box  definitions.  Box  m,  Box  n-2,  and  Box  n-1.  H-Boundary  b  has 
no  box  connections.  H-Boundary  c  has  two  box  connections,  Box  n  and  Box  n  + 1 .  Similarly, 
H-Boundary  d  has  two  box  connections.  Box  n + 2  and  Box  p.  Therefore,  similar  to  the  surface- 
hole  creation  boundary  connection,  the  box-hole  creation  boundary  connection  makes  holes 
in  meshes  specified  with  the  “MHOLEIN"  input. 

2.1.2.2  Direct  Hole  Specification 

The  third  method  for  defining  holes,  a  direct  method,  is  accomplished  through  the 
$REGION  and  $  VOLUME  user  inputs.  An  example  demonstrating  this  method  is  given  in 
Sec.  4.1.3.  The  SREGION  user  input  differs  from  the  $BOUNDARY  user  input  in  that  it 
is  connected  directly  to  only  one  mesh,  the  mesh  in  which  the  hole  will  be  located.  The 
ISPARTOF  user  input  specifies  this  mesh.  Typical  connections  between  regions  and  meshes 
are  shown  in  Fig.  11.  For  this  example,  Region  a  is  connected  to  and  makes  a  hole  in  Mesh 
n-2,  and  likewise.  Region  b  and  Region  c  are  connected  to  and  make  holes  in  Mesh  n  and 
Mesh  p,  respectively. 

Specification  of  the  ranges  of  the  J,  K,  and  L  indices  that  will  become  a  hole  is  accomplished 
through  the  SVOLUME  user  input.  In  Fig.  12  the  volume-region  connection  is  demonstrated. 
Similar  to  the  surface-boundary  connection  and  box-boundary  connection,  any  number  of 
volumes  can  be  used  to  create  a  region  (limited  by  the  parameter  NVDIM,  see  Sec.  3.1). 
For  this  case,  three  volumes  are  connected  to  Region  a.  Volume  m.  Volume  n-2,  and  Volume 
n-1 .  Region  b  and  Region  c  are  connected  to  one  volume  each.  Volume  n  and  Volume  n  -f- 1 , 
respectively.  Finally,  Region  d  has  two  volume  connections,  Volume  n-t-2  and  Volume  p. 

2.1.2.3  Determining  the  Hole  Boundary  (Fringe) 

Once  all  holes  in  all  meshes  have  been  created,  mesh  points  surrounding  holes  (i.e.,  hole 
boundary  points,  also  called  fringe  points  in  chimera  terminology)  must  be  identified.  Hole 
boundary  points  are  by  definition  the  point  or  points  separating  hole  points  from  field  points. 
Field  points  are  points  within  the  computational  domain  for  which  data  are  updated  by  the 
flow  solver  rather  than  by  interpolation  from  another  mesh.  To  rind  hole  boundary  points, 
each  point  in  a  mesh  is  checked.  If  the  point  is  a  hole  point,  nothing  is  done.  If  the  point 
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is  a  field  point  and  is  next  to  a  hole  point,  it  is  identified  as  a  hole  boundary  point.  This 
procedure  results  in  a  hole  boundary  with  only  one  interpolated  boundary  point  between 
a  hole  point  and  a  field  point .  This  type  of  hole  boundary  is  termed  a  single  fringe  boundary 
and  is  demonstrated  in  Fig.  13a.  Alternatively,  two  interpolation  points  can  be  specified 
between  a  hole  point  and  a  field  point.  This  type  of  hole  boundary  is  termed  a  double  fringe 
boundary  and  is  demonstrated  in  Fig.  13b. 

When  using  a  single  fringe  in  conjunction  with  a  flow  solver  which  uses  4th-order 
smoothing,  the  4th-order  terms  require  a  S-point  difference  stencil.  The  difference  stencil 
must  be  either  converted  to  2nd  order  (3-point  difference  stencil)  or  otherwise  reduced  at 
the  hole  boundaries.  This  reduction  to  2nd  order  is  not  necessary  for  a  double  fringe  boundary. 
Experience  has  shown  that  using  a  double  fringe  gives  better  results,  particularly  in  regions 
of  high  gradients.  However,  double  fringes  require  more  mesh  points  and  make  it  more  difficult 
to  obtain  the  required  overlap  between  hole  boundaries  and  outer  boundaries.  Overlap  is 
defined  as  those  field  points  lying  between  interpolation  boundary  points  contained  in  adjacent 
meshes. 

2.1.2.4  Common  Problems  in  Creating  Hoie  Boundaries 

Because  of  the  nature  of  the  algorithms  employed  in  PEGSUS  4.0,  undesirable  mesh 
features  such  as  internal  discontinuities  or  extreme  skewness  may  result  in  erroneous  results 
unless  proper  care  is  exercised.  Almost  all  errors  resulting  from  undesirable  mesh  features 
will  be  manifested  as  spurious  hole  points.  In  these  cases,  a  review  of  the  diagnostic  maps, 
described  in  Sec.  3.4  and  Sec.  4.0,  will  reveal  hole  points  occurring  in  regions  that  should 
not  be  blanked  out.  For  instance.  Fig.  14a  depicts  a  hole  creation  boundary  that  has  been 
defined  as  a  single  level  surface  within  a  mesh.  In  this  case,  the  level  surface  contains  a 
discontinuity  in  slope.  PEGSUS  4.0  will  attempt  to  identify  points  belonging  to  other  meshes 
that  fall  inside  the  hole  creation  boundary.  The  hole  location  procedure  will,  for  the  candidate 
point  (P)  shown,  identify  the  point  on  the  vertex  of  the  level  surface  as  the  closest  point 
to  the  candidate  point,  calculate  the  dot  product  of  the  normal  at  the  vertex  and  the  vector 
defined  by  the  vertex  and  candidate  points,  and  decide  that  the  candidate  is  inside  the  hole 
creation  boundary,  because  the  dot  product  found  in  the  previous  step  will  be  negative.  This 
determination  is  clearly  incorrect;  overall,  the  effect  of  the  discontinuity  in  the  normal  vectors 
will  be  a  “cloud”  of  hole  points  surrounding  the  vertex. 

Fortunately,  the  flexible  surface  and  boundary  definitions  incorporated  into  PEGSUS 
4.0  allow  such  problems  to  be  solved  by  modifying  the  input.  If  the  hole  creation  boundary 
is  divided  at  the  discontinuity  in  slope  into  two  separate  surfaces  (Fig.  14b),  then  the  candidate 
point  will  be  found  to  be  inside  surface  1  (solid  line),  outside  surface  2  (dashed  line),  and 
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therefore,  outside  the  hole  creation  boundary.  It  is  good  practice,  in  general,  to  construct 
hole  creation  boundaries  such  that  each  surface  comprising  the  boundary  is  smooth. 

Another  common  problem  occurs  in  discontinuous  regions  of  meshes.  Figure  ISa  depicts 
a  hole  creation  boundary  defined  as  a  level  surface  in  an  O-mesh  in  which  the  edges  of  the 
boundary  are  coincident  in  the  region  of  the  cut.  That  is,  the  first  and  last  points 
circumferentially  coincide  on  the  cut.  This  leads  to  a  discontinuity  in  slope  of  the  same  nature 
as  that  shown  in  Fig.  14.  Here  again  the  closest  points  to  a  candidate  point  (P)  are  the  coincident 
points  on  the  cut.  PEGSUS  will  simply  choose  one  of  the  points  as  the  closest  to  the  candidate 
(i.e.,  the  first  one  it  Hnds).  A  problem  arises  if  the  point  chosen  is  associated  with  a  normal 
vector  pointing  away  from  the  candidate  point;  in  this  case  the  candidate  will  be  determined 
to  be  inside  the  hole  creation  boundary.  Typically  this  problem  manifests  itself  as  a  funnel- 
shaped  **cloud”  of  hole  points  extending  from  the  cut  region  of  the  hole  creation  boundary 
into  the  interior  of  the  mesh  containing  the  candidate  point.  The  solution  of  this  problem 
is  similar  to  that  employed  in  the  previous  example.  If  the  hole  creation  boundary  is  divided 
into  two  surfaces  (solid  line  and  dashed  line),  as  shown  in  Fig.  ISb,  then  the  candidate  point 
will  be  found  to  be  outside  the  hole  creation  boundary. 

Another  example  is  illustrated  in  Fig.  16a  and  displays  the  concavity  problem  discussed 
in  Sec.  2. 1,2.1  and  shown  in  Fig.  6.  Two  meshes  are  included  in  this  configuration.  The  first 
mesh  models  a  body  with  an  attached  fin.  The  second  mesh  models  a  strut  attached  to  the 
body.  The  mesh  topology  of  the  body/fin  mesh  is  such  that  the  fin  and  body  surfaces  lie 
on  level  surfaces  that  are  normal  to  different  computational  coordinates.  It  is  necessary  to 
blank  out  points  in  the  strut  mesh  that  lie  within  the  body  and  fin.  Initially,  a  hole  creation 
boundary  is  created  from  the  two  level  surfaces  comprising  the  body  and  fin.  Application 
of  the  hole  location  algorithm  will  have  an  undesired  effect;  all  points  within  the  body  will 
be  inside  Surface  2  and  outside  Surface  1,  whereas  ail  points  within  the  fin  will  be  outside 
Surface  2  and  inside  Surface  1.  Since  a  point  must  be  inside  all  surfaces  to  be  considered 
inside  a  boundary  comprised  of  those  surfaces,  no  points  will  be  blanked  out. 

The  solution  to  this  problem  is  depicted  in  Fig.  16b,  where  the  original  hole  creation 
boundary  has  been  divided  into  two  separate  boundaries,  each  comprised  of  a  single  level 
surface.  Now  the  points  inside  the  fin  will  be  inside  Boundary  1,  and  points  inside  the  body 
will  be  inside  Boundary  2.  Note  that  the  fact  that  the  fin  points  are  outside  Boundary  2  has 
no  overall  effect,  since  a  point  is  blanked  if  it  is  inside  any  boundary. 

2.1.3  Outer  Boundary  Specification 

Outer  boundary  points  are  directly  specified  by  giving  the  J,  K,  and  L  index  ranges  of 
the  boundary  with  the  user  input  SSURFACE.  The  surface  or  group  of  surfaces  is  connected 


16 


AEDC-TR-91-8 


to  a  boundary  that  points  to  the  mesh  that  is  to  have  these  outer  boundary  points.  The  user 
input  ISPARTOF  points  to  the  mesh  that  contains  the  outer  boundary  points.  The  outer 
boundary  input  is  distinguishable  from  hole  boundary  input  in  that  no  MHOLEIN  input 
is  specified.  Shown  in  Fig.  17  are  outer  boundary-to-mesh  connections,  namely  O-Boundary 
a  and  O-Boundary  b  with  the  "ISPARTOF”  relationship  to  Mesh  n  and  Mesh  p,  respectively. 
Outer  boundaries  are  always  composed  of  level  surfaces  from  the  meshes  to  which  the 
boundaries  are  connected  (see  Fig.  18).  To  specify  a  double  interpolation  outer  boundary 
(used  when  a  double  fringe  is  specified),  two  level  surfaces  are  deGned  in  SSURFACE.  Outer 
boundary  specifications  are  demonstrated  in  the  examples  detailed  in  Sec.  4.1. 

2.2  INTERPOLATION 

Once  all  hole  boundary  and  outer  boundary  points  have  been  identified,  they  are  combined 
into  a  list  of  points  from  each  mesh  that  require  interpolation  of  flow-field  data.  At  this 
point,  PEGSUS  searches  other  meshes  in  the  composite  mesh  for  interpolation  stencils  that 
may  be  used  to  interpolate  flow-field  information  to  the  boundary  points  requiring  data. 

2.2.1  LINK  Lists 

Each  mesh  in  the  composite  mesh  receives  information  from  other  meshes  in  one  of  two 
ways;  (1)  through  its  outer  boundaries,  or  (2)  through  hole  boundaries.  To  determine  which 
meshes  are  to  communicate  with  each  other,  the  user  input  LINK  must  be  specified  within 
the  SMESH  user  input  (see  Sec.  3.2).  The  LINKs  are  priority  lists  that  indicate  the  order 
in  which  other  meshes  will  be  searched  for  suitable  interpolation  stencils.  A  typical  LINK 
for  a  mesh  is  shown  in  Fig.  19.  In  this  case  Mesh  n  will  search  Mesh  n  + 1 ,  Mesh  p.  Mesh 
n- 1,  and  Mesh  n  +  2,  in  that  order,  for  interpolation  stencils  that  can  provide  information 
to  Mesh  n’s  interpolation  boundary  points.  Each  mesh  in  the  composite  mesh  will  be  related 
to  the  other  meshes  by  similar  linkages.  The  maximum  number  of  LINKs  allowed  per  mesh 
is  set  by  the  program  parameter  LNDIM  (see  Sec.  3.1). 

The  LINK  definitions  for  each  mesh  comprise  a  priority  list  containing  the  names  of  other 
meshes  that  are  to  be  searched  for  interpolation  stencils  that  can  provide  information  for 
boundary  points.  The  meshes  in  the  priority  list  are  searched  in  the  order  in  which  they  appear 
in  the  input.  If  an  interpolation  stencil  is  found,  then  the  interpolation  information  for  the 
point  is  stored  in  a  set  of  arrays.  If  the  boundary  points  cannot  be  interpolated  by  a  stencil 
in  the  first  mesh  in  the  priority  list,  then  the  second  mesh  is  searched  for  interpolation  stencils. 
The  process  continues  until  either  interpolation  stencils  for  all  boundary  points  have  been 
found,  or  the  mesh  priority  list  is  exhausted.  Any  boundary  points  that  remain  after  the  mesh 
priority  list  is  exhausted  are  termed  "orphans.”  Orphaned  boundary  points  are  treated  as 
hole  points  (IBLANK  =  0).  The  user  must  decide  whether  to  accept  the  orphan  points  and 
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Update  them  in  the  flow  solver  (e.g.,  average  them  from  nearby  points)  or  to  change  the 
meshes  to  eliminate  the  orphan  points. 

To  reduce  the  computational  time  required  for  an  execution  of  PEGSUS  4.0,  the  LINK 
list  for  each  mesh  is  preprocessed  before  being  used  as  the  interpolation  priority  list.  First, 
a  parallelepiped  formed  by  the  minimum  and  maximum  X,  Y,  and  Z  values  of  each  mesh 
is  created.  Using  the  minimum/maximum  parallelepiped  for  each  mesh,  each  LINK  input 
of  a  mesh  is  checked  to  determine  if  it  intersects  with  the  mesh.  If  the  LINKed  mesh  does 
not  intersect  with  the  mesh,  the  LINKed  mesh  is  dropped  from  the  mesh*s  LINK  list.  For 
example,  if  ‘MESH  B’  is  a  LINK  of  ‘MESH  A’  and  the  parallelepiped  of  ‘MESH  A'  does 
not  intersect  with  the  parallelepiped  of  ‘MESH  B‘,  then  ‘MESH  B’  is  dropped  from  the  LINK 
list  of  ‘MESH  A’. 

Since  hole  location  and  boundary  point  interpolation  occur  as  two  separate  functions, 
a  mesh  that  causes  a  hole  in  another  mesh  does  not  necessarily  have  to  be  the  mesh  that 
provides  interpolated  information  to  the  resultant  boundary  points.  This  allows  meshes  to 
be  used  for  no  other  function  than  to  generate  holes  in  other  meshes;  this  feature  is  used 
extensively  in  certain  instances,  such  as  the  representation  of  a  fuselage  and  attached  wing 
(see  Sec.  2.3.3). 

2.2.2  Stencil  Searching 

Identifying  interpolation  stencils,  i.e.,  mesh  points  that  can  be  used  to  supply  interpolated 
information  to  boundary  points,  requires  a  search  procedure.  An  exhaustive  search  of  a  mesh 
will  always  find  a  permissible  interpolation  stencil  if  one  exists;  however,  exhaustive  searches 
are  computationally  intensive.  PEGSUS  uses  a  heuristic  search  mechanism  that  generally 
can  find  an  interpolation  stencil  in  relatively  few  steps. 

The  steps  required  for  obtaining  valid  interpolation  points  (or  a  stencil)  for  a  candidate 
boundary  point  are  described  below.  Briefly,  the  steps  are:  (1)  And  a  point  in  the  LINKed 
mesh  near  the  candidate  boundary  point,  (2)  find  a  stencii  with  satisfactory  interpolation 
coefficients,  and  (3)  check  that  the  stencil  is  valid.  A  candidate  boundary  point  is  given  three 
attempts  at  steps  2  and  3.  The  first  attempt  uses  the  stencil  found  for  the  previous  candidate 
boundary  point  as  the  near  point.  The  second  and  third  attempts  use  subsets  of  the  LINKed 
mesh  to  find  a  near  point.  The  near  point  is  determined  by  finding  the  smallest  distance  between 
the  candidate  point  and  a  point  in  the  subset  mesh.  The  first  subset  mesh  is  every  fourth 
point  in  each  of  the  J,  K,  and  L  directions.  The  second  subset  mesh  is  also  every  fourth  point 
but  shifted  by  two  points  in  each  of  the  J,  K,  and  L  directions.  Once  a  near  point  has  been 
found,  the  “stencil-walking”  procedure,  described  in  the  next  paragraph,  is  performed  on 
the  candidate  stencil. 
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The  interpolation  algorithm  used  in  PEGSUS  is  a  trilinear  interpolation  scheme,  and  is 
fully  documented  in  Ref.  4.  The  trilinear  interpolation  algorithm  yields  interpolation 
coefficients  for  a  point  either  inside  or  outside  the  trial  stencil  (the  latter  case  is  actually  an 
extrapolation).  The  interpolation  coefficients  will  all  be  between  0,0  and  1 .0  if  the  point  is 
inside  the  stencil;  a  point  outside  the  stencil  will  cause  at  least  one  coefficient  to  he  outside 
the  range  of  0.0  to  1.0.  The  interpolation  stencil  in  effect  defines  a  local  coordinate  system, 
while  the  interpolation  coefficients  define  the  location  of  the  point  in  the  local  system.  In 
PEGSUS,  the  interpolation  coefficients  are  used  to  determine  the  direction  (in  computational 
coordinates)  in  which  the  search  should  proceed  to  eventually  find  a  stencil  that  surrounds 
the  point.  For  instance,  a  given  boundary  point  and  a  stencil  with  lowest  index  comer  point 
at  J,  K,  and  L  may  yield  interpolation  coefficients  of  5.2,  —2.8,  and  1.3.  Pegsus  will  then 
move  to  the  stencil  corresponding  to  the  point  J  +  2,  K — 2,  and  L  + 1 .  (Note  that  the  J  index 
is  only  changed  by  2  since  the  procedure  limits  the  “walking”  to  2  indices  during  the  first 
10  iterations  and  1  index  for  the  final  10  iterations.)  The  new  stencil  will  usually  be  closer 
to  the  boundary  point.  The  procedure  is  repeated  until  all  interpolation  coefficients  are  between 
0.0  and  1.0,  which  is  considered  a  satisfactory  interpolation.  Although  this  searching 
mechanism  can,  in  principle,  fail  in  certain  cases  (e.g.,  highly  curved  or  distorted  meshes), 
it  is  highly  reliable  in  practice.  This  search  mechanism  is  referred  to  in  chimera  terminology 
as  “stencil  walking.”  The  method  of  "stencil  walking”  in  conjunction  with  finding  a  near 
point  in  the  subset  mesh  was  found  to  be  much  more  computationally  efficient  than  the 
“stencil-jumping”  procedure  used  in  PEGSUS  3.x  (Ref.  17). 

Once  a  satisfactory  interpolation  has  been  found,  the  stencil  from  which  this  interpolation 
is  made  must  be  checked  to  determine  if  the  stencil  is  valid.  The  difference  between  valid 
and  invalid  interpolations  is  illustrated  in  Fig.  20.  In  the  illustration,  the  outer  boundary 
of  a  mesh  (dotted  lines)  is  to  receive  interpolated  information  from  another  mesh  (solid  lines). 
(We  will  refer  to  the  first  mesh  as  a  “recipient”  mesh,  and  the  second  mesh  as  the  “donor” 
of  interpolated  information).  An  interpolation  is  valid  if  an  interpolation  stencil  can  be  found 
in  the  donor  mesh  for  which  none  of  the  corner  points  making  up  the  stencil  are  a  hole  or 
boundary  point,  or  a  boundary  point  is  coincident  with  a  corner  point  of  an  interpolation 
stencil  in  the  donor  mesh,  and  the  comer  point  is  not  a  hole  or  boundary  point.  The  hole 
boundary  point  in  Fig.  20  comprises  a  corner  of  two  stencils  that  contain  Points  A  and  B; 
as  a  result.  Points  A  and  B  cannot  receive  interpolated  information  from  the  donor  mesh. 
Point  C  is  coincident  with  a  corner  point  of  the  donor  mesh.  Since  the  corner  point  is  not 
itself  a  boundary  or  hole  point.  Point  C  can  receive  information  from  the  donor  mesh.  None 
of  the  corner  points  of  the  stencil  containing  Point  D  is  a  boundary  or  hole  point;  therefore. 
Point  D  will  be  permitted  to  receive  interpolated  information  from  the  donor  mesh.  In  Fig. 
20,  Points  A  and  B  will  be  orphans  unless  the  recipient  mesh  has  a  link  to  another  mesh 
that  can  provide  a  valid  interpolation  stencil.  PEGSUS  prints  a  list  of  orphan  points  if  any 
are  found;  orphans  are  also  identified  in  graphical  form  (See  Sec.  3.4)  for  convenience. 
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2.2.3  Quality 

The  preceding  rules  for  determining  a  valid  stencil  can  be  relaxed  by  use  of  the  QUALITY 
user  input.  QUALITY,  which  can  be  set  globally  or  input  for  each  mesh,  is  a  value  between 
0.0  and  1.0,  with  0.0  being  bad  and  1.0  being  good  (the  previous  discussion  in  Sec.  2.2.2 
corresponds  to  a  QUALITY  =  1.0  for  a  valid  stencil).  QUALITY  is  specified  by  three  values. 
The  first  value  is  the  initial  value,  the  second  value  is  the  ending  value  (which  must  be  less 
than  the  initial  value),  and  the  third  value  is  the  increment  (which  must  be  a  negative  value). 
Using  these  three  values,  ail  LINKed  meshes  are  searched  for  valid  interpolations  based  on 
the  initial  QUALITY  value.  If  any  interpolation  stencils  are  not  found,  the  minimum  acceptable 
value  for  QUALITY  is  changed  by  the  increment  value  and  the  search  for  valid  interpolations 
is  repeated  for  the  link  list.  The  process  is  continued  until  all  interpolation  stencils  are  found 
or  QUALITY  has  reached  the  ending  value. 

The  method  in  which  QUALITY  is  implemented  is  illustrated  in  Fig.  21 .  For  demonstration 
-purposes,  a  two-dimensional  stencil  is  used  with  interpolation  coefficients  in  the  xi  and  eta 
directions.  As  part  of  the  check  for  a  valid  interpolation,  each  stencil  point  is  checked  to 
determine  if  it  is  an  interpolation  point.  When  a  stencil  point  is  found  to  be  an  interpolation 
point,  the  coefficients  of  the  interpolation  are  checked  to  determine  if  they  are  valid  based 
on  QUALITY.  Figure  21a  demonstrates  this  check  for  QUALITY  =  0.8.  In  this  case,  the 
lower  left  corner  is  an  interpolation  boundary  point.  The  shaded  region  is  the  area  where 
interpoiation  coefficients  are  valid.  If  a  second  point  is  also  an  interpolation  boundary  point, 
then  the  region  for  valid  interpolation  coefficients  is  reduced  as  shown  in  Fig.  21b.  Finally, 
in  Fig.  21c  the  region  is  further  narrowed  for  three  interpolation  boundary  points  within 
the  stencil.  Of  course,  other  combinations  are  possible  when  three  dimensions  are  considered. 

2.3  ADDITIONAL  FEATURES 

Several  additional  features  are  implemented  in  PEGSUS  4.0  to  allow  greater  user  control 
of  PEGSUS  results.  The  features  discussed  in  the  following  sections  are  interior  blanking, 
inclusion,  phantom  meshes,  and  restarting  by  addition  and  subtraction. 

2.3.1  Interior  Blanking 

When  complex  configurations  are  modeled,  there  sometimes  arises  a  need  to  use  the 
blanking  capability  of  the  chimera  scheme  to  allow  points  within  the  interior  of  a  mesh  to 
be  updated  by  the  application  of  boundary  conditions  on  a  surface.  Typically  a  solid  wall 
geometry  internal  to  a  mesh  is  handled  in  this  way.  For  example,  the  fins  on  a  store  or  the 
corner  of  a  wall  can  be  handled  with  a  single  mesh  and  minimal  mesh  stretching  (see  Fig. 
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22).  This  interior  region,  which  is  referred  to  as  an  interior  blanked  region,  will  have  no 
interpolation  boundary  points  (fringe  points)  placed  around  the  blanked  area.  To  specify 
these  computational  regions  for  interior  blanking,  the  SREGION  input  is  used  with 
TYPE  =  TNTR’  (See  Sec.  3.2).  The  ranges  of  J,  K,  and  L  indices  of  the  blanked  region  are 
set  with  the  $VOLUME  user  input  (See  Sec.  3.2).  The  user  must  specify  the  correct  boundary 
conditions  on  the  surfaces  of  the  region  within  the  flow  solver  itself. 

Within  PEGSUS,  interior  blanked  points  are  handled  differently  from  the  points  blanked 
for  holes,  particularly  when  an  interpolation  stencil  is  being  checked  to  determine  if  it  is 
valid.  Valid  and  invalid  interpolation  stencils  near  interior  blanked  regions  (see  Sec.  2.2.2 
for  discussion  of  valid  and  invalid  interpolation)  are  shown  in  Fig.  23  with  two  interpolation 
points  that  are  considered  valid  (open  circles)  and  two  points  that  are  considered  invalid  (shaded 
squares).  As  can  be  seen,  an  interpolation  boundary  point  must  be  completely  enclosed  by 
interior  blanked  points  to  be  considered  invalid.  When  the  validity  of  an  interpolation  stencil 
in  or  near  interior  blanked  regions  is  checked,  QUALITY  (See  Sec.  2.2.3)  is  never  invoked. 
Interior  blanked  points  are  also  handled  differently  in  the  IBLANK  array  for  plotting  (see 
Sec.  3.4). 

2.3.2  Inclusion 

To  help  reduce  the  number  of  interpolations  and  possible  orphan  points,  user  inputs  can 
be  specified  to  enclose  the  physical  domain  globally;  in  addition  the  physical  and/or  the 
computational  domains  may  be  enclosed  for  each  mesh.  The  enclosure  of  the  domains  will 
exclude  boundary  points  outside  the  specified  ranges.  The  specification  of  the  domain  is 
accomplished  through  the  inputs  XINCLUDE,  YINCLUDE,  ZINCLUDE,  JINCLUDE, 
KINCLUDE,  and  LINCLUDE.  These  inputs  specify  either  the  physical  or  computational 
coordinates  within  which  ail  boundary  points  must  be  included.  Boundary  points  outside 
the  prescribed  coordinates  are  dropped  from  the  interpolation  point  lists.  This  capability 
allows  points  on  a  reflection  plane  to  be  excluded  from  the  interpolation  boundary  point 
list  and,  therefore,  reflection  plane  points  will  not  show  up  as  orphan  points. 

2.3.3  Phantom  Meshes 

Phantom  meshes  are  used  to  improve  and  simplify  the  creation  of  holes,  as  mentioned 
in  Sec.  2.2.1.  An  example  of  phantom  mesh  use  occurs  in  representing  a  wing/body  junction. 
A  complex  wing/body  junction  can  result  in  incorrect  hole  points  if  the  junction  is  not  properly 
represented.  Figure  24a  depicts  the  junction  of  a  wing  and  a  body.  The  boundary  of  the  wing 
mesh  extending  from  the  wing  root  is  coincident  with  the  surface  of  the  fuselage.  As  a  result, 
the  fuselage  surface  does  not  create  a  hole  in  the  wing  mesh.  The  wing  does,  however,  create 
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a  hole  in  the  fuselage  mesh.  A  hole  creation  boundary  defined  in  the  wing  mesh  cannot  be 
completely  closed;  an  aperture  in  the  hole  creation  boundary  will  exist  at  the  intersection 
of  the  wing  surface  and  the  hole  creation  boundary.  As  a  result,  some  points  close  to  the 
fuselage  may  be  considered  by  the  hole  location  algorithm  to  be  outside  the  hole  creation 
boundary;  these  points  will  not  be  blanked  as  required  (see  Fig.  24a}. 

A  solution  to  this  problem  is  illustrated  in  Fig.  24b.  The  structure  of  PEGSUS  4.0  allows 
a  mesh  to  serve  only  as  a  hole  creation  boundary.  Such  meshes  are  called  phantom  meshes 
and  do  not  interact  in  any  other  way  with  the  other  meshes;  i.e.,  they  do  not  provide 
interpolated  information  to  boundaries  of  other  meshes  nor  do  they  receive  information 
interpolated  from  other  meshes.  In  Fig.  24b,  the  hole  creation  boundary  shown  in  Fig.  24a 
has  been  used  to  define  the  outer  boundary  of  a  phantom  mesh.  However,  the  phantom  mesh 
extends  into  the  interior  of  the  fuselage  mesh.  The  outer  boundary  of  the  phantom  mesh 
is  defined  as  a  hole  creation  boundary  and,  therefore,  creates  a  hole  in  the  fuselage  mesh. 
In  effect,  the  phantom  mesh  acts  to  create  the  proper  hole  in  the  fuselage  mesh  for  the  wing 
mesh.  Since  the  hole  creation  boundary  of  the  phantom  mesh  is  completely  closed,  hole  points 
are  correctly  identified  in  the  fuselage  mesh. 

The  surfaces  of  the  fuselage  and  wing  do  not  create  holes  in  the  phantom  mesh,  since 
the  phantom  mesh  is  not  used  in  flow  computations.  PEGSUS  4.0  does  not  include  phantom 
meshes  in  the  composite  mesh  file.  However,  all  phantom  meshes  are  stored  for  restart 
executions  within  the  PEGSUS  restart  file.  Phantom  meshes  are  identified  by  PEGSUS  when 
no  LINK  input  is  specified  for  a  mesh. 

2.3.4  Restarting  —  Addition  and  Subtraction 

Applications  of  the  chimera  scheme  to  more  complex  configurations  have  increased  the 
number  of  meshes  required  as  well  as  the  computational  time  to  perform  the  functions  of 
PEGSUS.  PEGSUS  3.x,  being  a  batch  operation,  required  all  corrections  to  be  recalculated 
when  a  small  change  or  correction  was  made  to  the  user  inputs  or  meshes.  PEGSUS  4.0  includes 
a  restart  capability  to  allow  subtraction  and/or  addition  of  meshes,  boundaries,  and/or  regions. 
A  composite  mesh  and  associated  interpolation  data  can  therefore  be  incrementally  developed 
for  a  given  configuration. 

The  philosophy  used  in  PEGSUS  4.0  for  restart  conditions  is  to  recalculate  as  little  as 
possible.  This  requires  setting  up  an  internal  data  structure  that  allows  the  reuse  of  previous 
boundary  interpolation  data.  Previous  boundary  interpolation  stencils,  and  coefficients 
associated  with  subtracted  meshes  are  recalculated.  Interpolation  boundary  points,  stencils, 
and  coefficients  that  are  no  longer  required  are  removed  from  the  appropriate  lists.  Thus, 
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88  much  as  possible  of  the  previously  interpolated  boundary  results  is  used  when  calculating 
new  and/or  changed  boundaries. 

The  determination  of  the  interpolation  boundary  points  is  performed  differently  when 
adding  and  subtracting  a  component,  particularly  for  hole  boundaries.  The  addition  or 
subtraction  of  outer  boundaries  is  easily  determined  since  they  are  a  direct  specification  of 
the  boundary  points.  The  difficulty  comes  with  redefining  hole  boundaries  because  they  are 
created  by  several  different  methods  and  can  be  generated  by  one  or  more  meshes  or 
boundaries. 

The  process  used  for  adding  a  hole  is  depicted  in  Fig.  25  for  a  case  with  a  single  fringe. 
The  previous  hole  and  hole  boundary  shown  in  Fig.  25a  are  brought  into  PEGSUS  from 
the  previous  interpolation  file.  The  new  hole  and  hole  boundary  that  have  been  specified 
for  this  restart  execution  are  shown  in  Fig.  2Sb.  The  new  hole  and  hole  boundary  are  created 
in  a  temporary  array  within  PEGSUS.  Next,  the  new  hole  and  hole  boundary  are  superimposed 
on  the  previous  hole  and  hole  boundary.  The  result  is  illustrated  in  Fig.  2Sc.  Note  that  some 
of  the  previous  hole  boundary  points  are  unchanged,  while  others  have  become  hole  points 
and  must  be  removed  from  the  interpolation  data  file.  Also,  there  are  new  boundary  points 
requiring  interpolation. 

Subtraction  of  a  hole  and  hole  boundary  is  done  differently,  as  depicted  in  Fig.  26.  The 
hole  defined  in  Fig.  25c  is  shown  in  Fig.  26a.  The  hole  and  hole  boundary  defined  in  Fig. 
25a  are  to  be  subtracted  from  the  hole  and  hole  boundary  in  Fig.  26a.  The  new  resulting 
hole  and  hole  boundary  are  found  by  creating  a  list  of  points  that  includes  all  points  of  the 
previous  hole  and  hole  boundary,  searching  this  list  of  points  to  determine  which  points  are 
still  in  the  hole  from  hole  definitions  not  subtracted  (See  Sec.  2.1.2),  and  creating  a  fringe 
around  the  resulting  hole.  The  resulting  hole  and  hole  boundary  are  shown  in  Fig.  26b,  which 
is  the  same  hole  and  hole  boundary  shown  in  Fig.  2Sb.  Just  like  addition,  there  are  previous 
boundary  points  that  must  be  removed  from  the  interpolation  data  Hie  and  new  boundary 
points  that  require  interpolation. 

The  process  of  subtraction  and  addition  of  a  hole  and  hole  boundary  can  be  performed 
within  the  same  PEGSUS  run.  This  process  is  basically  a  combination  of  the  two  earlier 
processes  with  subtraction  being  performed  first  and  the  superimposition  of  the  added  hole 
coming  second.  Examples  of  addition  and  subtraction  are  given  in  Sec.  4.2  and  Appendix  G. 
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3.0  PROGRAM  IMPLEMENTAHON 

The  input  and  output  flies  for  PEGSUS  4.0  are  shown  in  Fig.  27.  The  following  is  a 
description  of  the  parameters,  NAMELIST  user  inputs,  input  mesh  format,  and  outputs  for 
PEGSUS  4.0. 

3.1  PARAMETERS 

The  following  constants  are  defined  in  PARAMETER  statements  at  the  beginning  of 
all  program  modules.  The  values  of  the  constants  must  be  equal  to  or  exceed  the  expected 
values  for  a  given  composite  mesh. 


Name 

Description 

MDIM 

maximum  number  of  meshes 

LNDIM 

maximum  number  of  LlNKs 

ICHAR 

maximum  number  of  characters  in  a  name  (must  be  40  characters 
or  less).  Note  that  this  parameter  also  controls  the  length  of  the 
mesh  name  read  from  the  mesh  input  file  (See  Sec.  3.3) 

NBDIM 

maximum  number  of  boundaries 

MHDIM 

maximum  number  of  MHOLEIN 

NSDIM 

maximum  number  of  surfaces 

MSLEN 

maximum  number  of  points  defining  a  single  surface 

NXDIM 

maximum  number  of  boxes 

NRDIM 

maximum  number  of  regions 

NVDIM 

maximum  number  of  volumes 

INCORE 

=  0  external  storage  used  (i.e.,  SSD  or  disk) 

=  1  internal  storage  used  (i.e.,  CRAY  2) 

MLEN 

maximum  number  of  points  in  a  single  mesh 
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MLEMAX  if  INCORE  =  0,  MLEMAX  =  MLEN 

if  INCORE  =  1,  total  number  of  points  in  all  meshes 

MILEN  maximum  number  of  interpolation  or  stencil  points  in  a  single  mesh 

NBYTEI  number  of  bytes  for  an  integer  value 

NBYTER  number  of  bytes  for  a  real  value 

(NBYTEI  and  NBYTER  are  used  to  set  up  the  record  lengths  of  the  direct  access  files  when 
INCORE  =  0.) 

The  parameters  MSSMAX,  MIMAXl,  and  MIMAX2  are  calculated  based  on  earlier  defined 
values  and  should  not  be  changed. 

As  may  be  noted  in  the  parameters,  PEGSUS  4.0  has  the  ability  to  keep  all  run  time 
memory  within  core  or  to  utilize  the  disk  (i.e.,  for  computers  with  limited  memory).  To  keep 
all  run  time  memory  within  core  the  INCORE  parameter  must  be  set  equal  to  one,  MLEMAX 
must  be  set  to  the  total  of  all  mesh  points,  and  MLEN  equal  to  the  maximum  sized  mesh. 
To  utilize  disk  memory,  set  INCORE  equal  to  zero,  MLEN  equal  to  the  maximum  sized 
mesh,  and  MLEMAX  equal  to  MLEN.  Using  INCORE  =  1  is  preferable  (though  typically 
not  obtainable)  since  the  I/O  of  the  program  is  reduced. 

3.2  NAMELIST  GROUPS  AND  RECORD  DEFINITIONS 

The  NAMELIST  input  is  divided  into  seven  groups:  SGLOBAL,  SMESH,  $BOUNDARY, 
SSURFACE,  $BOX,  SREGION,  and  $VOLUME.  Note  that  all  user  input  record  names  are 
capitalized  and  all  alphanumeric  names  (i.e.,  NAME,  MHOLEIN,  LINK)  are  case  sensitive. 
The  group  names  and  associated  records  are  defined  as  follows: 

SGLOBAL 

This  input  sets  global  values  for  all  meshes.  Each  record  can  be  overridden  for  any  individual 
mesh  by  using  the  SMESH  input.  However,  the  mesh  override  is  limited  to  that  individual 
mesh.  For  subsequent  meshes,  the  global  values  are  used  unless  specifically  overridden.  The 
SGLOBAL  inputs  are  recognized  and  implemented  during  an  initial  execution  of  PEGSUS 
4.0  but  not  for  a  restart  execution  (i.e.,  SGLOBAL  cannot  be  changed  for  a  restart  execution). 

EPS  —  Tolerance  of  interpolation  coefficient.  Value  of  interpolation  coefficients  less 
than  0  and  greater  than  1.0  that  is  allowed  for  a  converged  stencil  (default  =  0.001). 
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FRINGE  —  Set  equal  to  I  for  single  fringe  (default)  or  equal  to  2  for  double  fringe  (See 
Sec.  2. 1.2.3). 

QUALITY  —  Value  deflning  quality  of  interpolation  stencil  ranging  from  0.0  (bad)  to 
I.O  (good),  (default  =  1.0)  (See  Sec.  2.2.3).  Three  values  are  entered  in  the  foliowing  order: 
the  initial  value,  the  ending  value,  and  the  increment  (which  must  be  negative). 

XINCLUDE,  YINCLUDE,  ZINCLUDE  —  Range  of  X,  Y,  and  Z  values,  respectively, 
to  be  included  as  boundary  points  (See  Sec.  2.3.2).  Composed  of  two  values  (a  minimum 
and  a  maximum  value)  in  each  of  the  X,  Y,  and  Z  directions,  respectively. 

$MESH 

ADDHOLB  —  List  of  hole  creation  boundaries  defined  in  the  previous  PEGSUS  execution 
that  are  to  make  a  hole  in  this  added  mesh  (must  be  any  alphanumeric  characters  with  no 
special  symbols).  This  input  is  used  only  for  restart  when  MODE  =  ‘ADD*.  Each  boundary 
name  is  a  string  not  exceeding  ICHAR  characters.  No  default  values  are  defmed. 

ADDLINK  —  List  of  meshes  defined  in  the  previous  PEGSUS  execution  that  must  be 
linked  to  this  added  mesh  (must  be  any  alphanumeric  characters  with  no  special  symbols). 
The  named  mesh  (specified  with  NAME)  is  put  at  the  top  of  the  LINK  list  of  each  mesh 
specified  by  ADDLINK.  This  input  is  used  only  for  restart  when  MODE  =  *ADD’.  Each  mesh 
name  is  a  string  not  exceeding  ICHAR  characters.  No  default  values  are  deHned. 

ALFA,  BETA,  GAMA  —  Euler  angles  (degrees)  determining  rotation  of  input  meshes. 
ALFA,  BETA,  and  GAMA  are  rotations  about  the  Z,  Y,  and  X  coordinates,  respectively. 
Default  values  for  ALFA,  BETA,  and  GAMA  are  0.0.  The  coordinates  of  a  mesh  are 
transformed  in  the  following  sequence  of  operations:  scaling,  translation,  and  rotation. 

EPS  —  Tolerance  of  interpolation  coefficient.  Value  of  interpolation  coefficients  less 
than  0  and  greater  than  1.0  that  is  allowed  for  a  converged  stencil.  Default  is  {GLOBAL  value. 

FRINGE  —  Set  equal  to  1  for  single  fringe  or  equal  to  2  for  double  fringe.  Default  is 
SGLOBAL  value  (See  Sec.  2. 1.2. 3). 

JINCLUDE,  KINCLUDE,  LINCLUDE  —  Range  of  J,  K,  and  L  indices,  respertively, 
to  be  included  for  boundary  points  (See  Sec.  2.3.2).  Composed  of  two  values  (a  starting  and 
an  ending  value)  in  each  of  the  J,  K,  and  L  directions,  respectively. 
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LINK  —  Priority  list  of  meshes  that  are  to  be  searched  to  provide  interpolations  for 
boundary  points  (must  be  any  alphanumeric  characters  with  no  special  symbols).  Each  LINK 
is  a  string  not  exceeding  ICHAR  characters.  LINKs  are  specified  as  follows:  LINK  =  ‘mesh 
name  1’,  ‘mesh  name  2',  etc.,  (to  link  all  meshes  specify  LINK=,) 

MODE  —  Defines  a  mesh  that  is  added  (=  ‘ADD’),  or  subtracted  (  =  ‘SUB’).  MODE 
is  specified  only  for  restart. 

NAME  —  Name  of  mesh  to  which  NAMELIST  record  refers  (must  be  any  alphanumeric 
characters  with  no  special  symbols).  A  NAME  is  a  string  not  exceeding  ICHAR  characters. 
No  default  values  are  defined. 

QUALITY  —  Value  defining  quality  of  interpolation  stencil  ranging  from  0.0  (bad)  to 
1.0  (good),  (default  is  SGLOBAL  value)  (See  Sec.  2.2.3).  Three  values  are  entered  in  the 
following  order;  the  initial  value,  the  ending  value,  and  the  increment  (which  must  be  negative). 

SCALE  —  Mesh  scaling  factor.  The  original  mesh  coordinates  are  multiplied  by  the  scaling 
factor.  The  default  value  is  1 .0,  The  coordinates  of  a  mesh  are  transformed  in  the  following 
sequence  of  operations:  scaling,  translation,  and  rotation. 

XINCLUDE,  YINCLUDE,  ZINCLUDE  —  Range  of  X,  Y,  and  Z  values,  respectively, 
to  be  included  for  boundary  points  (See  Sec.  2.3.2).  Composed  of  two  values  (a  minimum 
and  a  maximum  value)  in  each  of  the  X,  Y,  and  Z  directions,  respectively. 

XR,YR,ZR  —  Point  (in  translated  and  scaled  coordinates)  about  which  input  mesh  is 
to  be  rotated.  Default  values  are  0.0.  The  coordinates  of  a  mesh  are  transformed  in  the 
following  sequence  of  operations:  scaling,  translation,  and  rotation. 

X0,Y0,Z0  —  Mesh  translation  factors.  These  factors  are  added  to  the  original  mesh 
coordinates.  Default  values  are  0.0.  The  coordinates  of  a  mesh  are  transformed  in  the  following 
sequence  of  operations;  scaling,  translation,  and  rotation. 

SBOUNDARY 

CLOSED  —  Either  .TRUE,  (default)  or  .FALSE.  Set  equal  to  .FALSE,  when  the  surfaces 
do  not  describe  a  completely  closed  boundary  (See  Sec.  2. 1.2.1). 

ISPARTOF  —  Mesh  which  contains  boundary  (must  be  any  alphanumeric  characters 
with  no  special  symbols).  ISPARTOF  is  a  string  not  exceeding  ICHAR  characters.  No  default 
is  defined. 
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MHOLEIN  —  Meshes  in  which  the  named  boundary  causes  holes.  MHOLEIN’s  are 
specified  as  follows:  MHOLEIN  *  ‘mesh  name  1*,  ‘mesh  name  2’,  etc.  Each  mesh  name 
in  MHOLEIN  is  a  string  not  exceeding  ICHAR  characters.  No  defaults  are  defined. 

MODE  —  Defines  a  boundary  that  is  added  (=  ‘ADD’),  or  subtracted  (=  ‘SUB’),  MODE 
is  specified  only  for  restart.  If  boundary  is  connected  to  added  mesh,  MODE  is  not  required. 
If  boundary  is  connected  to  subtracted  mesh,  SBOUNDARY  need  not  be  specified, 

NAME  —  Name  of  boundary  (must  be  any  alphanumeric  characters  with  no  special 
symbols).  A  boundary  name  is  a  string  not  exceeding  ICHAR  characters.  No  default  is  defined. 

$SURFACE 

ISPARTOF  —  Boundary  to  which  surface  belongs  (mast  be  any  alphanumeric  characters 
with  no  special  symbols).  ISPARTOF  is  a  string  not  exceeding  ICHAR  characters.  No  default 
is  defined. 

JRANGE,  KRANGE,  LRANGE  —  Ranges  of  indices  that  define  surface.  Composed 
of  two  values  (a  starting  and  an  ending  value)  in  each  of  the  J,  K,  and  L  directions,  respectively. 
Maximum  allowable  ranges  are: 

JRANGE:  1..JMAX 
KRANGE:  1..KMAX 
LRANGE:  1..LMAX 


No  defaults  are  defined. 

MESH  NAME  —  Name  of  mesh  that  contains  the  surface  (ISPARTOF)  (must  be  any 
alphanumeric  characters  with  no  special  symbols).  MESH  NAME  is  used  only  for  hole  creation 
boundaries.  Name  is  a  string  not  exceeding  ICHAR  characters.  Default  is  the  mesh  name 
given  in  ISPARTOF  for  $BOUNDARY. 

NVOUT  —  Direction  of  normal  outward  from  hole.  NVOUT  is  required  only  for  hole 
creation  boundaries.  Allowable  values  are  “  + J”,  “-J”,  “  +  K”,  “-K”,  “  +  L”,and“-L”. 

$BOX 

ISPARTOF  —  Boundary  to  which  box  belongs  (must  be  any  alphanumeric  characters 
with  no  special  symbols).  ISPARTOF  is  a  string  not  exceeding  ICHAR  characters.  No  default 
is  defined. 
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XRANGE,  YRANGE,  ZRANGE  —  Ranges  of  X,  Y,  and  Z  coordinates  that  deflne  box. 
Composed  of  two  values  (a  minimum  and  a  maximum  value)  in  each  of  the  X,  Y,  and  Z 
directions,  respectively. 

SREGION 

ISPARTOF  —  Mesh  which  contains  region  (must  be  any  alphanumeric  characters  with 
no  special  symbols).  ISPARTOF  is  a  string  not  exceeding  ICHAR  characters.  No  default 
is  defined. 

MODE  —  Defines  a  region  that  is  added  (=  ‘ADD’),  or  subtracted  (  =  ‘SUB*).  MODE 
is  specified  only  for  restart.  If  region  is  connected  to  added  mesh,  MODE  is  not  required. 
If  region  is  connected  to  subtracted  mesh,  SREGION  need  not  be  specified. 

NAME  —  Name  of  region  (must  be  any  alphanumeric  characters  with  no  special  symbols). 
A  region  name  is  a  string  not  exceeding  ICHAR  characters.  No  default  is  defined. 

TYPE  —  Defines  type  of  region  being  specified.  If  TYPE  =  ‘HOLE’,  then  hole  with  a 
fringe  is  created.  If  TYPE  =  ‘INTR’,  then  an  interior  region  will  be  blanked  out  but  with 
no  fringe.  Default  is  ‘HOLE’. 

SVOLUME 

ISPARTOF  —  Region  to  which  volume  belongs  (must  be  any  alphanumeric  characters 
with  no  special  symbols).  ISPARTOF  is  a  string  not  exceeding  ICHAR  characters.  No  default 
is  defined. 

JRANGE,  KRANGE,  LRANGE  —  Ranges  of  indices  that  define  volume.  Composed 
of  two  values  (a  starting  and  an  ending  value)  in  each  of  the  J,  K,  and  L  directions,  respectively. 
Maximum  allowable  ranges  are: 

JRANGE:  1..JMAX 
KRANGE:  1..KMAX 
LRANGE:  1..LMAX 


No  defaults  are  defined. 

3.3  INPUT  MESH  FORMAT 

The  format  of  the  mesh  input  Hie  is  identical  to  that  of  PEGSUS  3.x.  The  input  format 
uses  three  records.  These  records  are  repeated  for  each  mesh. 
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Record  Variables 

1  Mesh  Name 

2  JMAX,  KMAX,  LMAX 

3  (((X(J,K,L),  J  =  1 .  JMAX),K  =  1  ,KMAX),L  =  1  ,LMAX), 
(((Y(J,K,L),J=  1,  JMAX), K  =  1, KMAX), L=  1, LMAX), 
(((Z(J,K,L),  J  =  1 ,  JMAX),K  =  1  ,KMAX).L  =  1  ,LMAX) 

where  JMAX,  KMAX,  and  LMAX  are  the  maximum  mesh  indices  in  the  J,  K,  and  L 
directions,  respectively.  Mesh  Name  must  be  ICHAR  in  length. 

3.4  OUTPUT 

There  are  seven  output  files  generated  by  PEGSUS  4.0  during  each  execution  (see  Fig. 
27).  The  files  are;  (1)  composite  mesh,  (2)  interpolation  data,  (3)  execution  summary,  (4) 
diagnostic  maps,  (S)  PEGSUS  restart,  (6)  IBLANK  array  for  plotting,  and  (7)  orphan  point  list. 

The  composite  mesh  flle,  interpolation  data  file,  and  the  orphan  point  list  file  are  used 
by  the  flow  solver.  The  structure  of  these  files  is  given  in  Appendix  F.  The  composite  mesh 
file  and  the  interpolation  data  file  are  required  for  a  restart  of  PEGSUS. 

The  execution  summary,  which  is  written  to  logical  unit  6,  contains  (1)  a  summary  of 
the  meshes  read,  (2)  a  summary  of  user  inputs,  (3)  error  messages,  (4)  an  orphan  point 
summary,  and  (5)  flow  solver  parameters.  Items  1, 2, 4,  and  5  will  be  illustrated  in  Sec.  4.0. 
Error  checking  is  performed  by  PEGSUS  on  the  inputs  specified  by  the  user  and  on  array 
overruns.  First,  key  words  (user  inputs,  i.e.,  NAME,  FRINGE,  ISPARTOF,  etc.)  are  checked 
to  determine  if  they  are  recognized  by  the  program.  If  a  key  word  is  not  recognized,  an  error 
message  is  printed  and  the  value  is  not  read  by  PEGSUS.  The  second  level  of  checking 
determines  if  connections  (i.e.,  ISPARTOF  and  MHOLEIN)  specified  by  the  user  point  to 
existing  meshes,  boundaries,  or  regions.  Again,  if  an  error  is  found,  a  message  defining  the 
problem  is  printed.  The  third  level  of  checking  determines  if  an  array  is  overrun.  The  error 
message  for  these  errors  indicates  what  parameter  needs  to  be  changed  to  correct  the  problem. 
PEGSUS  tries  to  find  all  user  input  error  before  terminating  the  processing. 

The  diagnostic  maps  generated  by  PEGSUS  are  diagrams  of  computational  space  showing 
holes,  boundary  points,  and  interpolation  stencils.  These  maps  and  their  usefulness  will  be 
illustrated  in  Sec.  4.0. 
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The  PEGSUS  restart  file  contains  all  user  inputs  and  calculated  values  from  the  current 
PEGSUS  run.  This  file  is  required  to  restart  PEGSUS.  The  structure  of  this  file  is  given  in 
Appendix  F. 

The  IBLANK  array  for  plotting  is  a  file  used  in  conjunction  with  the  mesh  files  for  use 
with  PLOT3D  (Ref.  19)  and  FAST  (Ref.  20).  The  IBLANK  array  for  plotting  is  different 
from  the  IBLANK  array  sent  to  the  flow  solver  in  that  the  boundary  points  are  not  set  equal 
to  zero  but  equal  to  minus  the  mesh  number  of  the  donor  mesh  (see  Sec.  2.3.1).  Field  points 
are  still  equal  to  1  and  hole  points  equal  to  0.  For  interior  blanked  points  the  surfaces  (minimum 
and  maximum  index  planes)  of  the  region  are  set  equal  to  one  with  the  points  within  the 
region  set  equal  to  zero.  The  difference  in  setting  the  IBLANK  array  for  plotting  is  to  allow 
values  to  be  plotted  on  the  surfaces  where  the  boundary  conditions  are  set.  The  structure 
of  this  file  is  given  in  Appendix  F. 


4.0  EXAMPLES 

Six  examples  of  the  application  of  PEGSUS  4.0  to  representative  multiple-mesh  problems 
are  presented  in  Secs.  4. 1  and  4.2.  These  examples  are  not  intended  to  be  illustrations  of 
the  application  of  PEGSUS  4.0  to  actual  engineering  problems.  Instead,  they  are  intended 
to  demonstrate  most  of  the  basic  principles  of  PEGSUS  4.0  that  are  common  to  a  majority 
of  PEGSUS  4.0  scenarios.  A  set  of  examples  of  the  application  of  PEGSUS  4.0  to  an  actual 
engineering  problem  (an  empty  cavity,  a  store/sting  combination  in  a  cavity,  and  a  store 
in  a  cavity)  is  presented  in  Appendix  G. 

In  Sec.  4.1  four  examples  are  given,  the  first  three  of  which  illustrate  in  turn  each  of 
the  different  methods  available  for  hole  creation.  The  fourth  example  consists  of  a  hole  creation 
boundary  comprised  of  surfaces  from  two  different  meshes.  Section  4.2  consists  of  two 
examples  designed  to  illustrate  the  use  of  the  addition  and  subtraction  features  of  PEGSUS 
4.0.  Common  to  all  examples  is  the  use  of  a  double  fringe.  It  is  recommended  that  prospective 
PEGSUS  4.0  users  carefully  study  the  annotated  input  and  diagnostic  maps  included  for  all 
six  examples  to  gain  a  solid  understanding  of  the  process. 

4.1  HOLE  CREATION  EXAMPLES 

The  following  four  examples  illustrate  various  means  of  hole  creation  through  the  use 
of  simple  cylindrical  and  Cartesian  mesh  configurations. 
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4.1.1  Indirect  Hole  Creation  Using  SBOUNDARY  and  SSURFACE 

A  simple  two-mesh  problem  is  depicted  in  Fig.  28.  A  right  circular  cylinder  of  finite  length 
(surrounded  by  a  cylindrical  mesh)  is  embedded  in  a  Cartesian  mesh.  The  outer  boundary 
of  the  cylindrical  mesh  receives  data  interpolated  from  the  surrounding  Cartesian  mesh.  Since 
the  cylinder  is  embedded  in  the  Cartesian  mesh,  a  region  of  the  Cartesian  mesh  must  be 
excluded  from  the  computation.  The  exclusion  is  achieved  by  a  hole  creation  boundary  defined 
by  a  set  of  surfaces  within  the  cylindrical  mesh.  To  illustrate  the  flexibility  of  PEGSUS  4.0, 
a  hole  creation  boundary  is  defined  that  consists  of  two  half  cylinders  of  unequal  radii.  If 
the  hole  creation  boundary  were  to  be  defined  as  a  collection  of  four  surfaces  (the  two  half- 
cylindrical  surfaces  and  the  two  rectangular  surfaces  between  the  larger  and  smaller  radii), 
the  hole  boundary  would  contain  concavities.  For  reasons  described  in  Sec.  2. 1.2.1,  defining 
a  hole  boundary  with  concavities  will  result  in  an  incorrect  hole  in  the  Cartesian  mesh.  The 
correct  way  to  describe  the  desired  hole  boundary  is  as  two  separate  hole  boundaries,  neither 
of  which  contains  any  concavities.  In  this  case,  one  hole  boundary  is  defined  as  the  larger 
half-cylindrical  and  two  rectangular  surfaces  extending  to  the  surface  of  the  cylinder;  the 
other  hole  boundary  consists  of  the  smaller  half-cylindrical  and  two  corresponding  rectangular 
surfaces.  The  two  hole  boundaries  will  merge  into  a  single  hole  created  in  the  Cartesian  mesh. 

The  annotated  input  for  the  problem  of  Fig.  28  is  given  in  Fig.  29  (see  Sec.  3.2  for  a 
description  of  user  input).  Two  mesh,  three  boundary,  and  eight  surface  definitions  are 
required.  Maximum  index  values  for  both  the  cylindrical  and  Cartesian  meshes  are  41, 41, 
and  5  in  the  J,  K,  and  L  directions,  respectively.  The  hole  boundary  defined  by  the  half¬ 
cylinder  of  smaller  radius  is  specified  by  four  surfaces,  while  the  hole  boundary  defined  by 
the  half-cylinder  of  larger  radius  is  specified  by  three  surfaces.  The  outer  boundary  of  the 
cylindrical  mesh  is  defined  by  a  single  surface.  No  boundaries  or  surfaces  of  the  Cartesian 
mesh  are  defined  in  the  input,  since  no  predefined  surface  in  the  Cartesian  mesh  receives 
information  from  the  cylindrical  mesh. 

When  PEGSUS  4.0  is  run,  an  execution  summary  file  is  produced  as  well  as  diagnostic 
maps  which  are  intended  to  aid  the  user  in  determining  whether  the  mesh  communications 
are  being  performed  properly.  Figure  30  is  a  reproduction  of  a  map  produced  from  the  input 
in  Fig.  29.  The  maps  are  graphical  depictions  of  the  computational  meshes  which  show:  (1) 
which  points  are  blanked  out  hole  points,  (2)  which  points  are  boundary  points  and  which 
mesh  updates  the  boundary  points,  (3)  which  points  are  stencil  reference  points  (i.e.,  the 
point  in  the  interpolation  stencil  with  lowest  J,  K,  and  L  indices)  and  which  meshes  the 
interpolation  stencils  update,  and  (4)  which  points  are  “orphans,"  i.e.,  interpolation  boundary 
points  which,  for  some  reason,  cannot  be  interpolated  from  other  meshes.  At  the  beginning 
of  the  diagnostic  maps  is  a  legend  giving  a  key  to  the  letters  and  symbols  used.  The  map 
in  Fig.  30  clearly  shows  the  location  of  the  hole  in  the  Cartesian  mesh  as  well  as  indicating 


32 


AEOC-TR-91-8 


that  the  hole  boundary  points  are  being  updated  by  the  cylindrical  mesh.  Also,  the  location 
of  stencil  reference  points  in  the  Cartesian  mesh  that  are  interpolating  information  to  the 
outer  boundary  of  the  cylindrical  mesh  are  shown.  With  a  little  practice,  a  user  can  obtain 
much  useful  information  from  the  maps,  such  as  the  amount  of  mesh  overlap,  and  whether 
interpolation  stencils  have  been  found  for  boundary  points.  For  example,  it  can  be  seen  from 
the  map  depicted  in  Fig.  30  that  the  amount  of  mesh  overlap  (i.e.,  the  distance  from  the 
hole  boundary  points  to  outer  boundary  points  of  the  cylindrical  mesh)  is  marginal  near  the 
bottom  of  the  hole  boundary.  Shown  in  Fig.  31  is  a  portion  of  the  execution  summary  for 
this  problem.  The  execution  summary  is  only  included  for  this  example. 

4.1.2  Indirect  Hole  Creation  Using  {BOUNDARY  and  $BOX 

PEGSUS  4.0  provides  a  great  deal  of  flexibility  for  the  creation  of  holes.  Figure  32  depicts 
the  same  two-mesh  system  used  in  Sec.  4.1 .1.  In  this  instance,  the  hole  creation  boundaries 
have  been  replaced  by  a  box.  The  box  consists  of  ranges  of  X,  Y,  and  Z  values  (the  minimum 
and  maximum  X  and  Y  ranges  are  shown  in  Fig.  32)  which  define  a  region  within  the  cylindrical 
mesh.  All  points  in  the  Cartesian  mesh  contained  within  this  box  are  considered  hole  points. 
A  hole  boundary  (fringe)  is  placed  around  the  hole  to  specify  which  points  require  interpolation. 
The  remaining  points  are  considered  Held  points.  The  annotated  input  for  this  problem  is 
shown  in  Fig.  33.  Figure  34  shows  a  portion  of  the  diagnostic  maps  that  are  output  by  PEGSUS 
4.0  for  this  problem.  The  hole  in  the  Cartesian  mesh  created  by  the  box  is  clearly  visible. 

4.1.3  Direct  Hole  Creation  Using  {REGION  and  {VOLUME 

In  many  situations,  one  would  like  to  directly  exclude  a  specific  region  (or  regions)  in 
a  particular  mesh.  PEGSUS  4.0  allows  this  process  through  the  use  of  the  region/volume 
combination.  Figure  35  illustrates  just  such  a  situation  with  the  familiar  Cartesian/cylindrical 
mesh  system.  Here  the  points  to  be  blanked  consist  of  a  single  region,  with  the  volume  being 
the  J,  K,  and  L  values  of  the  points  in  the  Cartesian  mesh  which  are  to  be  blanked.  The 
NAMELIST  input  for  this  problem  is  shown  in  Fig.  36  and  a  sample  diagnostic  map  for 
the  Cartesian  mesh  is  shown  in  Fig.  37.  Examination  of  the  map  shows  that  a  hole  boundary 
has  been  placed  around  the  hole  created  by  the  region.  This  boundary  would  be  omitted 
if  the  input  for  the  region  had  included  the  TYPE  =  TNTR’  statement.  In  the  absence  of 
this  input,  the  region  type  was  defaulted  to  ^HOLE’. 

4.1.4  Hole  Creation  Using  Surfaces  from  Two  Separate  Meshes 

When  computing  the  flow  field  about  complex  geometric  configurations,  it  is  often 
desirable  to  break  the  mesh  surrounding  an  object  with  a  simple  topology,  such  as  a  cylindrical 
mesh,  into  separate  meshes.  The  need  to  do  this  frequently  arises  from  the  fact  that  a  single 
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mesh  would  be  too  large  for  the  computer  memory  available,  although  this  is  certainly  not 
the  only  cause.  Breaking  a  single  mesh  into  two  or  more  separate  meshes  can  lead  to  difficulties 
when  the  split  mesh  configuration  is  embedded  within  other  meshes.  In  earlier  versions  of 
PEGSUS,  the  surfaces  defining  a  hole  creation  boundary  were  required  to  be  part  of  the 
same  mesh  to  which  the  hole  creation  boundary  was  connected.  This  would  lead  to  problems 
for  a  split  mesh  system  as  the  separate  hole  creation  boundaries  from  each  mesh  would  not 
create  the  hole  that  was  needed.  A  phantom  mesh  would  then  be  required  to  obtain  the  correct 
hole.  This  shortcoming  has  been  alleviated  in  PEGSUS  4.0  by  allowing  the  surfaces  comprising 
a  hole  creation  boundary  to  come  from  meshes  other  than  the  one  to  which  the  hole  creation 
boundary  is  connected.  To  illustrate  this  feature,  we  again  turn  to  the  Cartesian  mesh  and 
cylindrical  mesh  combination.  In  this  case  the  cylindrical  mesh  surrounding  the  cylinder  has 
been  split  into  two  pieces  as  is  shown  in  Fig.  38.  The  PEGSUS  4.0  input  is  shown  in  Fig. 
39.  Examination  of  the  PEGSUS  input  reveals  that  a  single  hole  creation  boundary  has  been 
defined,  and  that  this  boundary  is  part  of  ‘CYLINDER  GRID  1*.  This  boundary  is  composed 
of  two  surfaces;  one  surface  from  ‘CYLINDER  GRID  T  and  one  surface  contained  in 
‘CYLINDER  GRID  2*.  The  maps  produced  by  PEGSUS  4.0  for  the  L  =  1  plane  are  shown 
in  Fig.  40. 

4.2  ADDITION  AND  SUBTRACTION  EXAMPLES 

A  major  enhancement  of  PEGSUS  4.0  over  earlier  versions  of  the  PEGSUS  code  Is  the 
ability  to  add  and  subtract  meshes,  boundaries,  and  regions  to  multiple-mesh  configurations 
which  have  already  been  processed  by  PEGSUS  4.0  in  a  previous  run.  With  earlier  versions 
of  PEGSUS,  the  only  means  by  which  objects  could  be  added  or  subtracted  from  a  multiple- 
mesh  configuration  was  to  rerun  PEGSUS  with  the  corresponding  lines  of  input  added  or 
removed.  The  entire  configuration  then  was  reprocessed  by  the  code.  This  was  very  costly 
and  time  consuming  when  only  minor  changes  were  required.  With  PEGSUS  4.0  the  situation 
is  different.  Information  from  an  earlier  execution  of  the  code  is  stored  in  such  a  way  that 
it  may  be  reused  as  input  for  later  runs.  In  this  manner,  when  one  wishes  to  make  changes 
to  a  multiple-mesh  configuration,  any  information  generated  by  an  earlier  run  of  the  PEGSUS 
4.0  code  is  recomputed  only  if  the  newly  added  or  subtracted  features  require  it  to  be 
recomputed.  This  capability  leads  to  substantial  savings  in  the  time  required  to  construct 
a  complex  mesh  configuration.  The  following  two  examples  illustrate  the  use  of  the  addition 
and  subtraction  features  of  PEGSUS  4.0. 

4.2.1  Addition 

To  demonstrate  the  addition  feature  of  PEGSUS  4.0,  the  example  in  Sec.  4.1.1  is 
reconsidered  (Figs.  28  -  31).  It  is  assumed  that  the  composite  mesh  Hie,  the  interpolation 
file,  and  the  PEGSUS  restart  file  are  available  for  input  to  PEGSUS  in  the  restart  mode 
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(see  Fig.  27).  In  this  example,  another  Cartesian  mesh  will  be  added  to  the  existing  two-mesh 
system.  The  new  three-mesh  configuration  is  illustrated  in  Fig.  41,  where  the  J  and  K 
coordinates  are  the  same  as  in  Fig.  28.  Since  the  new  Cartesian  mesh  overlaps  both  the  original 
Cartesian  mesh  and  the  cylindrical  mesh,  PEGSUS  4.0  must  have  some  way  to  incorporate 
the  new  mesh  interactions.  Examination  of  Fig.  41  shows  that  an  outer  boundary  must  be 
added  to  the  original  Cartesian  mesh  to  allow  information  to  be  transferred  from  the  new 
Cartesian  mesh  to  the  original  Cartesian  mesh.  It  must  be  recalled  that  in  PEGSUS  terminology 
an  outer  boundary  specifically  refers  to  those  points  in  a  mesh  other  than  hole  boundary 
points  which  receive  data  interpolated  from  other  meshes.  In  addition,  the  new  Cartesian 
mesh  must  be  added  to  the  MHOLEIN  list  of  the  hole  boundaries  created  by  the  cylindrical 
mesh  and  also  must  be  added  to  the  link  lists  of  the  other  two  meshes.  The  PEGSUS  4.0 
input  for  this  problem  is  shown  in  Fig.  42.  The  outer  boundary  added  to  the  original  Cartesian 
mesh  contains  the  MODE  =  ‘ADD’  statement  so  that  PEGSUS  will  add  the  boundary.  The 
new  Cartesian  mesh  is  added  to  the  MHOLEIN  list  for  the  hole  boundaries  created  by  the 
cylindrical  mesh  through  inclusion  of  the  line  ADDHOLE= ‘CYLINDER  HOLE 
BOUNDARY  r,  ‘CYLINDER  HOLE  BOUNDARY2’.  The  new  Cartesian  mesh  is  also  added 
to  the  link  lists  of  the  original  Cartesian  mesh  and  cylindrical  mesh  through  inclusion  of 
the  line  ADDLINK  = ‘CARTESIAN  GRID’,  ‘CYLINDER  GRID’.  The  maps  output  by 
PEGSUS  4.0  for  the  L  =  1  plane  for  all  three  meshes  are  shown  in  Fig.  43. 

4.2.2  Subtraction 

For  the  last  example,  the  cylindrical  grid  will  be  subtracted  from  the  previous  example 
of  Sec.  4.2.1  (Figs.  41  -  43).  The  PEGSUS  input  for  this  problem  is  shown  in  Fig.  44,  while 
the  maps  for  the  Cartesian  grids  for  the  L  =  1  plane  are  shown  in  Fig.  45.  An  examination 
of  the  maps  shows  that  the  cylindrical  grid  has  been  removed  as  well  as  the  boundaries 
contained  within  the  mesh.  Note  that  the  hole  has  been  removed  from  both  Cartesian  meshes. 
For  generality,  note  that  the  processes  of  addition  and  subtraction  may  be  performed 
simultaneously. 


5.0  CONCLUDING  REMARKS 

PEGSUS  4.0  is  a  code  used  in  the  chimera  multiple-mesh  flow-field  calculation  scheme. 
It  has  been  developed,  implemented,  and  applied  successfully  to  aerodynamic  conEgurations. 
PEGSUS  4.0  is  an  improvement  over  earlier  versions  in  that  it  contains  a  restart  capability 
that  allows  for  addition  and/or  subtraction  of  meshes,  boundaries  and/or  regions,  the  user 
has  more  control  of  the  process  through  additional  inputs,  the  code  was  written  to  be  portable, 
the  code  is  4  to  10  times  faster  than  PEGSUS  3.x,  and  requires  30  percent  less  memory  than 
PEGSUS  3.x. 
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FEGSUS  4.0  may  be  obtained^  by  contacting  AEDC/DOCS,  Arnold  Engineering 

Development  Center,  Arnold  AFB,  TN  37389. 
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Figure  2.  Overlap  region  between  meshes. 


40 


Individual 

Meshes 

PEGSUS 

Interpolation 

Data 

Flow 

Solver 

User 

Composite 

Inputs 

Mesh 

Flow  Solution 

- » 


Figure  3.  Chimera  scheme. 
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a.  Area  eliminated  by  SI 


b.  Area  eliminated  by  SI  and  S2 


c.  Area  eliminated  by  SI,  S2,  and  S3 
Figure  5.  Hole  location. 
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a.  Hole  creation  boundary  with  concavity. 


Hole  Creation  Boundary  No.  1 


Figure  6.  Restrictions  on  hole  creation  boundaries. 
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Figure  7.  Hole  creation  boundary-mesb  connections. 
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Figure  8.  Surface-hole  creation  boundary  connections. 


AEDC-TR-91-8 


CLOSED  =  TRUE 


CLOSED  =  FALSE 


b.  “CLOSED’*  has  no  effect 

Figure  9.  Effects  of  the  “CLOSED”  user  input  on  hole  boundaries. 
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Figure  10.  Box-hole  creation  boundary  connections. 
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Figure  11.  Regloii>mesh  connections. 
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Figure  12.  Volume-region  connections. 
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a.  Single  fringe 

Figure  13.  Hole  boundary  generation. 
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b.  Double  fringe 
Figure  13.  Concluded. 
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a.  Problem:  spurious  hole  points  near  sharp  comer 


Surface  2 


b.  Solution:  divide  hole  creation  boundary  into  two  surfaces 
Figure  14.  Hole  creation  boundary  with  sharp  comer. 
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a.  Problem:  spurious  hole  points  near  coincident  points 


b.  Solution:  divide  hole  creation  boundary  into  two  surfaces 
Figure  15.  Single  boundary  with  coincident  points. 
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b.  Solution:  divide  hole  creation  boundary  into  two  boundaries 
Figure  16.  Concave  hole  creation  boundary. 
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figure  17.  Outer  boundary-mesh  connections. 
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Figure  19. 
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QUALITY  =  0.8 


Valid 

Interpolation 
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c.  Three  blanked  points  in  stencil 
Figure  21.  Concluded. 
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b.  Comer 

Figure  22.  Interior  blanking. 
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Figure  23.  Valid  and  invalid  interpolation  near  an  interior  blanked  region. 
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a.  Problem:  unblanked  bole  points  near  wlng/body  junction 
Figure  24.  Compiex  wing/body  Junction. 
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M  Region  in  Fuselage  Mesh 
Outside  of  Hole  &eation 
Boundary 


of  Phantom  Mesh 

b.  SoluUon:  use  of  phantom  mesh  to  create  holes 
Figure  24.  Concluded. 
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a.  Previous  hole  and  hole  boundary 
Figure  25.  Addition  of  hole  and  hole  boundary. 
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b.  New  hole  uid  hole  boundary 
Figure  25.  Continued. 
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Combined  hole  and  hole  boundary 
Figure  25.  Concluded. 
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a.  Hole  and  hole  boundaiy 
Figure  26.  Subtraction  of  hole  and  hole  boundaiy. 
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b.  Resulting  hole  and  hole  boundary 
Figure  26.  Concluded. 
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Figure  27.  PEGSUS  4.0  input  and  output  files. 
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Figure  28.  Graphical  depiction  of  indirect  hole  creation  using 
SBOUNDARY  and  SSURFACE. 
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c  FEGSUS  4.0  INPUT 

C  CARTESIAN  GRID  /  CYLINDER 

C  EXAMPLE  4.1.1  -  INDIRECT  HOLE  SPECIFICATION 

C  DEFINITION  OF  GLCfflAL  PARAMETERS 
$CLOBAL 

FRINGE  -  2, 

$END 

C  MESH  definitions 

SHESH  name  -  'CARTESIAN  GRID', 

LINK  -  , 

SEND 

SME5H  NAME  •  'CYLINDER  GRID', 

LINK  -  'CARTESIAN  GRID', 

XO  -  20.0, 

XO  -  20.0, 

SEND 

C  BOUNDARY  DEFINITIONS  CONSIST  OF, TWO  HOLE 
C  CREATION  BOUNDARIES  AND  AN  OUTER  BOUNDARY. 

SBOUNDARY  NAME  -  'CYLINDER  HOLE  BODHDARYl', 
I5FARTOF  “  'CYLINDER  GRID', 

MHOLEIN  -  'CARTESIAN  GRID', 

SEND 

SBOUNDARY  NAME  -  'CYLINDER  BOLE  BOONDARY2', 
ISFARTOF  -  'CYLINDER  GRID', 

MHOLEIN  -  'CARTESIAN  GRID', 

SEND 

SBOUNDARY  NAME  =  'CYLINDER  OUTER  BOUNDARY', 
ISPARTOF  -  'CYLINDER  GRID', 

SEND 

C  SURFACE  DEFINITIONS 

C  FOUR  SURFACES  ARE  NEEDED  TO  DESCRIBE 
C  CYLINDER  HOLE  BOUHDARYl . 

C  THE  FIRST  TWO  SURFACES  DESCRIBE  THE  HALF 
C  CYLINDER  OF  SMALLER  RADIOS. 

C  THE  REMAINING  TWO  SURFACES  DESCRIBE  THE 
C  RECTANGULAR  SEGMENTS. 

SSURFACE  ISPARTOF  “  'CYLINDER  HOLE  BOUNDARYl', 
ORANGE  -  11,11, 

KRANGE  -  1,11, 

LRANGE  -  1,5, 

KVOUT  -  '  +  J'  , 

SEND 

SSURFACE  ISPARTOF  -  'CYLINDER  HOLE  BOUNDARYl', 
ORANGE  o  11,11, 

KRANGE  <•  31,41, 

LRANGE  -  1,5, 

NVOUT  -  '+J' , 

SEND 

SSURFACE  ISPARTOF  ■  'CYLINDER  HOLE  BOUNDARYl', 
ORANGE  -  1,11, 

KRANGE  -  11,11, 

LRANGE  -  1,5, 


a.  Page  1 

Figure  29.  Input  for  indirect  hole  creation  using 
SBOUNDARY  and  SSURFACE. 
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NVOUT  =  '+K'  , 

$EHD 

^SURFACE  ISPARTOF  a  'CYLINDER  BOLE  BOUMDARYl'/ 
JRAH6E  -  1,11, 

KRAH6E  ■>  31,31, 

LRAH6E  -  1,5, 

HVOUT  -  '-K', 

$END 

C  THREE  SURFACES  ARE  HEEDED  TO  DESCRIBE 
C  CYLINDER  BOLE  BODNDARY2 . 

C  THE  FIRST  SURFACE  DESCRIBES  THE  HALF 
C  CYLINDER  OF  LARGER  RADIUS. 

C  THE  REMAINING  TWO  SURFACES  DESCRIBE  THE 
C  RECTANGULAR  5E«IEHTS. 

$ SURFACE  ISPARTOF  »  'CYLINDER  HOLE  BOUNDARY2' , 
JRANGE  -  21,21, 

KRANGE  -  11,31, 

LRAHGE  -  1,5, 

HVOUI  *•  '+J', 

SEND 

SSURFACE  ISPARTOF  -  'CYLINDER  HOLE  BOUNDARYB' , 
ORANGE  B  1,21, 

KRANGE  -  11,11, 

LRAHGE  B  1,5, 

NVOUT  -  '-K', 

SEND 

SSURFACE  ISPARTOF  >  'CYLINDER  HOLE  BOUNDARY!' , 
JRAHGE  B  1,21, 

KRANGE  B  31,31, 

LRANGE  -  1,5, 

NVOUT  -  '+K', 

SEND 

C  THE  LAST  SURFACE  IS  THE  OUTER  BOUNDARY  OF 
C  THE  CYLINDER  GRID.  NO  NORMAL  VECTOR  IS 
C  SPECIFIED  FOR  SURFACES  WHICH  ARE  CONNECTED 
C  TO  OUTER  BOUNDARIES  THROUGH  THE  ISPARTOF 
C  USER  INPUT. 

SSURFACE  ISPARTOF  ■>  'CYLINDER  OUTER  BOUNDARY', 
ORANGE  B  40,41, 

KRANGE  B  1,41, 

LRAHGE  B  1,5, 

SEND 

C  END  OF  PEGSDS  4.0  INPUT  FOR  EXAMPLE  4.1.1 

b.  Page  2 

Figure  29.  Concluded. 
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F^ure  30.  Diagnostic  map  for  indirect  hole  creation  using 
SBOUNDARY  and  SSURFACE. 
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PFPPFFPPP  EEEFP-EEEE  ggggggggg  sssssssss  uon  uou  sssssssss 
PPPPPPPPP  EEEEEEEBE  GGGGGGGGG  SSSSSSSSS  DUG  UDD  SSSSSSSSS 
PPP  PPP  EEE  GGG  SSS  UDD  DUD  SSS 

FPP  PPP  EEE  GGG  SSS  UUU  UDD  SSS 

PPPPPPPPP  EEEEEEEEE  GGG  SSSSSSSSS  DUD  DUD  SSSSSSSSS 

PPPPPPPPP  EEEEEEEEE  GGG  SSSSSSSSS  DUD  DUD  SSSSSSSSS 

PPP  bee  GGG  GGGG  SSS  UDD  DDU  SSS 

PPP  EEE  GGG  G6  SSS  DDU  UDD  SSS 

PPP  EEEEEEEEE  GGGGGGGGG  SSSSSSSSS  DDDUDDUUD  SSSSSSSSS 

PPP  EEEEEEEEE  GGGGGGGGG  SSSSSSSSS  DDDUDUDDD  SSSSSSSSS 


VERSION  4.0 

1991 


GLOBAL  DESCRIPTION  READ 

2  MESH  DESCRIPTIONS  READ 

3  BOUNDARY  DESCRIPTIONS  READ 
B  SURFACE  DESCRIPTIONS  READ 
0  BOX  DESCRIPTIONS  READ 

0  REGION  DESCRIPTIONS  READ 
0  VOLUME  DESCRIPTIONS  READ 


READING  MESH:  CARTESIAN  GRID 
READING  HESBi  CYLINDER  GRID 


GLOBAL  PARAMETER  SUMMARY 


MIN/HAX  OF  DOMAIN 
MIN  X  -  O.OOOOE-fOO  MAX  X 

MIN  Y  -  O.OOOOE+00  MAX  Y 

MIN  Z  -  O.OOOOE-l-00  MAX  Z 


0.4000E+02 
0 . 4000E+02 
0.4000E+01 


NO  INCLUDE  RANGES  SPECIFIED 


HO.  OF  FRINGE  POINTS 
QUALITY (1) 

QUALITY (2) 

QUALITY  (3) 

EPS 


1.0000 

1.0000 

1.0000 

O.lOOOE-02 


a.  Page  1 

Figure  31.  Execution  summary  for  indirect  hoie  creation  using 
(BOUNDARY  and  (SURFACE. 
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MESa  DAME:  CARTESIAN  SRID 


REGSUS  MESH  HO.-  1 
COMPOSITE  HBSa  HO.-  1 


INDEX  RANGES:  JHAX 

TRANSLATIOH:  XO 
ROTATKH):  XR 

ALPHA  - 
SCALING:  SCALE  - 


41  KMAX 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+OO 

O.lOOOE+01 


41  LHAX  -  5 


YO  - 

Q.OOOOE4-00 

zo  - 

O.OOOOB+00 

YR 

O.OOOOE-^-OO 

ZR  - 

0,OOOOE4-O0 

BETA  - 

O.OODOE+00 

GAM  - 

O.OOOOE+OO 

HO.  or  FRINGE  POINTS 
Ot)ALITY41) 

QUALITY <2) 

QUALITY  (3) 

EPS 

MIH/HAX  OF  MESH 
HIH  X  i>  O.OOOOB+OO 

HIM  Y  -  0.0000E4'OO 

MIH  Z  -  O.OODOE'fOO 


1.0000 

1.0000 

1.0000 

O.lOOOE-02 


MAX  X  -  0.4000E'I-02 
MAX  Y  -  0.4000E4-02 
MAX  Z  -  0.4000E'I>01 


-  2 


HO  INCLUDE  RANGES  SPECIFIED 


LINKS 

NO.  NAME 

1  CYLINDER  GRID 


MESR  NAME:  CYLINDER  GRID 


PEG5US  MESH  NO.-  2 
COMPOSITE  MESH  HO.-  2 


INDEX  RANGES:  JHAX  - 

41  KMAX 

-  41 

LHAX  -  S 

TRANSLATION;  XO 

Q.20D0E+02 

YO 

0.2000E+02 

ZO  - 

O.OOOOE+OO 

ROTATION:  XR 

O.OOOOE+OO 

YR  - 

O.OOOOE+OO 

ZR  - 

O.OOOOE+OO 

ALPHA  - 

O.OOOOE+OO 

BETA  - 

O.OOOOE+OO 

GAM  - 

O.OOOOE+OO 

SCALING:  SCALE  - 

O.IOOOB+Ol 

NO.  OF  FRINGE  POINTS 

ES 

2 

tfJALITYd) 

- 

1.0000 

QUALITY (2) 

- 

1.0000 

QUALITY (3) 

■ 

1.0000 

EPS 

B 

O.IOOOE- 

02 

MIN/HAX  OF  MESH 

MIN  X  -  0.2000E+01 

HAX  X  - 

0.3800E+02 

MIN  Y  -  0.2000E+01 

MAX  Y  - 

0,3800E+02 

MIN  Z  -  O.OOOOE+OO 

MAX  Z  - 

0.4000E+01 

NO  INCLUDE  RANGES  SPECIFIED 


LINKS 

HO.  NAME 

1  CARTESIAN  GRID 


b.  Page  1 

Figure  31.  Continued. 
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*•**•  boundary  parameter  summary 


BOUNDARY  NAME:  CYLINDER  HOLE  BOUNDARYl 

BOUNDARY  TYPE  -  BOLE 
CLOSED  -  T 

IS  PART  OF  MESH  -  CYLINDER  GRID 

MIN /MAX  or  HOLE  BOUNDARY 

MIN  X  -  0.2000E+02  MAX  X  °  0.2825E+02 

MIN  Y  =  0,1175B+02  MAX  Y  =  0.2825E+02 

HIM  Z  °  O.OOOOE+00  MAX  Z  »  0.4000E'f01 

MAKE  HOLE  IN: 

NO.  NAME 
1  CARTESIAN  GRID 


BOUNDARY  HAKE:  CYLINDER  HOLE  BOUNDARYZ 

BOUNDARY  TYPE  =  HOLE 
CLOSED  ••  T 

IS  PART  OP  MESH  "  CYLINDER  GRID 

HIH/MAX  OF  HOLE  BOUNDARY 

HIM  X  n  0.8500E+01  MAX  X  -  0.2000E+02 

KIN  Y  ^  0.8SOOE4-01  MAX  Y  »  0.S150E+02 

MIN  Z  >  O.DOOOE-fOO  MAX  Z  0.4000E+01 

MAKE  HOLE  IN; 

HO.  NAME 
1  CARTESIAN  GRID 


BOUNDARY  HAHE:  CYLINDER  OUTER  BOUNDARY 

BOUNDARY  TYPE  -  OUTER 

CLOSED  =  T 

IS  PART  OF  HESB  -  CYLINDER  GRID 


SURFACE  NO.  »  1 

ISFARTOF  :  CYLINDER  HOLE  BOUNDARYl 

SURFACE  RANGES 

MIN  Jo  11  MAX  Jo  11 

MIH  K  -  1  MAX  K  o  11 

HIM  L  -  1  MAX  L  o  5 

NVOUT  =  +J 

CONNECTING  MESH  •>  CYLINDER  GRID 


SURFACE  NO.  -  2 

ISPARTOF  :  CYLINDER  HOLE  BOUNDARYl 

SURFACE  RANGES 

HIH  Jo  11  MAX  Jo  11 

MIN  K  o  31  MAX  K  ^  41 

MIN  L  o  1  MAX  L  o  5 

NVOUT  -  +J 

CONNECTING  MESH  o  CYLINDER  GRID 

c.  Page  3 

Figure  31.  Continued. 
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SURFACE 

NO.  B 

3 

ISPARTOF  :  CYLINDER 

HOLE 

BOnNDARYl 

SURFACE 

RANGES 

MIN  J  - 

1 

MAX 

J  = 

11 

MIH  K  B 

11 

MAX 

K  B 

11 

HIN  L  B 

1 

MAX 

L  = 

5 

HVODT  -  +K 

CONMECTIHG  HESH  -  CYLINDER  GRID 


SURFACE  HO.  -  4 

ISPARTOP  :  CYLINDER  HOLE  BOUNDARY! 

SURFACE  RANGES 

MIN  J  -  1  MAX  J  -  11 

MIM  K  =  31  MAX  K  -  31 

HIM  L  -  1  MAX  L  -  5 

NVOUT  «■  -K 

CONNECTING  MESH  -  CYLINDER  GRID 


SURFACE  MO.  -  5 

ISPARTOF  ;  CYLINDER  HOLE  BOUMDARY2 

SURFACE  RANGES 

MIN  J  -  21  MAX  J  -  21 

MIN  K  -  11  MAX  K  >  31 

MIN  L  ~  1  MAX  Lb  S 

NVOUT  -  +J 

CONNECTING  MESH  CYLINDER  GRID 


SURFACE  HO.  B  6 

ISPARTOF  ;  CYLINDER  BOLE  BOUHDARY2 

SURFACE  RANGES 

HIH  J  -  1  MAX  J  -  21 

MIN  K  -  11  MAX  Kb  11 

HIH  Lb  1  max  L  -  5 

HVODT  «  -K 

COHMECTING  MESH  >  CYLINDER  GRID 


SURFACE  HO.  -  7 

ISPARTOF  !  CYLINDER  HOLE  B0UnDARX2 

SURFACE  RANGES 

HIH  J  B  1  MAX  J  B  21 

HIM  X  B  31  max  K  -  31 

HIM  Lb  1  max  Lb  5 

NVOUT  B 

CONNECTING  MESH  b  CYLINDER  GRID 


SURFACE  NO.  -  8 

ISPARTOF  I  CYLINDER  OUTER  BOUNDARY 

SURFACE  RANGES 

HZN  J  B  40  MAX  J  -  41 

MIN  K  B  1  MAX  Kb  41 

MIN  Lb  1  max  Lb  5 

NVOUT  - 

CONNECTING  HESH  -  CYLINDER  GRID 

d.  Page  4 

Figure  31.  Continued. 
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INTERPOLATING  BOUNDARY  POINTS  FOR 
HESB;  CARTESIAN  GRID 
LINK:  CYLINDER  GRID 

QUALITY-  1.00000  STENCILS  FOUND:  630 

NO  ORPHAN  POINTS  LEFT  IN  MESH  -  CARTESIAN  GRID 


INTERPOLATING  BOUNDARY  POINTS  FOR 
MESH:  CYLINDER  GRID 
LINK:  CARTESIAN  GRID 

QUALITY-  1.00000  STENCILS  FOUND:  410 

HO  ORPHAN  POINTS  LEFT  IN  MESH  -  CYLINDER  GRID 


**•»*  quality  summary 


HIHIMUM 
QUALITY 
-1.0 
>-0.9 
>-0.8 
>-0.7 
>-0.6 
>-0.S 
>-D.4 
>-0.3 
>-0.2 
>-0.1 
>-0.0 

MESH I  CARTESIAN  GRID 

NUMBER  OF  BOUNDARY  POINTS;  630 

HUMBER  OF  INTERPOLATION  STENCILS:  4l0 

MESH:  CYLINDER  GRID 

number  of  boundary  points:  410 

NUMBER  OF  INTERPOLATION  STENCILS:  630 

e.  Page  5 

Fignre  31.  Condnded. 


NO.  OF  INTERPOLATED 
boundary  POINTS  FOUND 
1040 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
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Box  Used  to  Create  a  Hole  in 
the  Cartesian  Mesh 


Figure  32.  Graphical  depiction  of  indirect  hole  creation  nslng 
SBOUNDARY  and  $BOX. 
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FEGSOS  4.0  INPUT 
CARTESIAN  GRID  /  CYLINDER 
EXAMPLE  4.1.2  -  DIRECT  HOLE  SPECIFICATION 

C  DEFINITION  OF  GLOBAL  PARAMETERS 
SGLOBAL 

FRINGE  >’2, 

SEND 

C  MESH  DEFINITIONS 

SMESH  NAME  »  'CARTESIAN  GRID', 

LINK  °  , 

SEND 

SMESH  NAME  °  'CYLINDER  GRID', 

LINK  »  'CARTESIAN  GRID', 

XO  =  20,0, 

YO  ->  20.0, 

SEND 

C  BOUNDARY  DEFINITIONS  CONSIST  OF  ONE  HOLE 
C  CREATION  BOUNDARY  AND  AN  OUTER  BOUNDARY. 

SBOUNDARY  NAME  ••  ' CYLINDER  HOLE  BOUNDARY' , 
ISPARTOF  -  'CYLINDER  GRID', 

MHOLEIN  •*  'CARTESIAN  GRID', 

SEND 

SBOUNDARY  NAME  =  'CYLINDER  OUTER  BOUNDARY', 
ISPARTOF  =  'CYLINDER  GRID', 

SEND 

C  BOX  DEFINITION 

C  ONLY  ONE  BOX  IS  DEFINED.  THIS  BOX  IS 
C  PART  OF  THE  CYLINDER  BOLE  BOUNDARY 

SBOX  ISPARTOF  =  'CYLINDER  HOLE  BOUNDARY' 
XRANGE—  13.0,24.0, 

YRANGE  -  13.0,24.0, 

ZRANGE  =  0.0, 4.0, 

SEND 

C  SURFACE  DEFINITION 

C  THE  ONLY  SURFACE  IS  THE  OUTER  BOUNDARY  OF 
C  THE  CYLINDER  GRID.  NO  NORMAL  VECTOR  IS 
C  SPECIFIED  FOR  SURFACES  HHICH  ARE  CONNECTED 
C  TO  OUTER  BOUNDARIES  THROUGH  THE  ISPARTOF 
C  USER  INPUT. 

SSURFACE  ISPARTOF  =•  'CYLINDER  OUTER  BOUNDARY', 
ORANGE  -  40,41, 

E(RANGE  1,41, 

LRANGE  -  1,5, 

SEND 

C  END  OF  PEGSUS  4.0  INPUT  FOR  EXAMPLE  4.1.2 

Figure  33.  Input  for  indirect  hole  creation  using 
SBOUNDARY  and  SBOX. 
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DIAGNOSTIC  MAPS 

LEGEND 

BOUNDARY  POINTS  UPDATED  BY  ; 

A  -  CARTESIAN  GRID 

B  -  CYLINDER  GRID 

INTERPOLATION  STENCILS  OPDATING  POINTS  IS  ! 
a  -  CARTESIAN  GRID 

b  -  CYLINDER  GRID 

-  FIELD  POINT 
«  .  hole  POINT 
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Figure  34.  Dii^ostic  map  for  indirect  hole  creation  using 
SBOUNDARY  and  $BOX. 
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Figure  35>  Graphical  depiction  of  direct  hole  creation  using  $R£GION  and  $VOLlJME. 
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C  FEG30S  4.0  INPUT 

C  CARTESIAN  GRID  /  CTLIMDER 

C  EXAMPLE  4.1.3  -  DIRECT  HOLE  SPECIFICATION 

C  DEFINITION  OF  GLOBAL  PARAMETERS 
SGLOBAL 

FRINGE  -  2, 

SEND 

-  C  MESH  DEFINITIONS 

SMESH  NAME  -  'CARTESIAN  GRID', 

LINK  «  , 

SEND 

$M£SB  NAME  -  'CYLINDER  GRID', 

LINK  -  'CARTESIAN  GRID', 

XO  -  20.0, 

YO  =  20.0, 

SEND 

C  BOUNDARY  DEFINITION  CONSISTS  OF  ONE  OUTER  BOUNDARY. 

$ boundary  name  -  'CYLINDER  OUTER  BOUNDARY', 

ISPARTOF  -  'CYLINDER  GRID', 

Send 

C  REGION  DEFINITION 

C  ONLY  ONE  REGION  IS  DEFINED.  THIS  REGION 
C  WILL  MAKE  A  HOLE  IN  THE  CARTESIAN  GRID. 

SREGIOH  name  -  'CARTESIAN  HOLE  BOUNDARY', 

ISPARTOF  =»  'CARTESIAN  GRID', 

SEND 


C  SURFACE  DEFINITION 

C  THE  ONLY  SURFACE  IS  THE  OUTER  BOUNDARY  OP 
C  THE  CYLINDER  GRID.  NO  NORMAL  VECTOR  IS 
C  SPECIFIED  FOR  SURFACES  WHICH  APE  CONNECTED 
C  TO  OUTER  BOUNDARIES  THROUGH  THE  ISPARTOF 
C  USER  INPUT. 

$ SURFACE  ISPARTOF  =  'CYLINDER  OUTER  BOUNDARY', 
ORANGE  »  40,41, 

KRANGE  -  1, 41, 

LRANGE  -  1,5, 

SEND 

C  VOLUME  DEFINITION 

C  ONLY  ONE  VOLUME  IS  DEFINED.  THE  VOLUME 
C  CONSISTS  OF  THE  J  ,  K  AND  L  RANGES  OF 
C  THE  POINTS  IN  THE  CARTESIAN  GRID  WHICH 
C  ARE  TO  BE  BLANKED  OUT. 

SVOLUHE  ISPARTOF  -  'CARTESIAN  HOLE  BOUNDARY' 
ORANGE  -  14,25, 

KRANGE  -  14,25, 

LRANGE  1,5, 

SEND 


C  END  OF  PEGSUS  4.0  INPUT  FOR  EXAMPLE  4.1.3 

Figure  36.  Input  for  direct  hole  creation  using  SREGION  and  $VOLUME. 


86 


AEDC-TR-91-8 
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Figure  37.  Diagnostic  map  for  direct  bole  creation  using 
$REGION  and  SVOLUME. 
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a.  ^CYLINDER  GRID  1’  and  ‘CARTESIAN  GRID’ 
Figure  38.  Graphical  depiction  of  indirect  hole  creation  using 
surfaces  from  two  separate  meshes. 
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b.  ‘CYLINDER  GRID  2’  and  ‘CARTESIAN  GRID’ 
Figure  38.  Concluded. 
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C  PEGSUS  4.0  IHPUT 

C  CABTESIAH  GRID  /  CYLINDER 

C  EXAMPLE  4.1.4  -  BOLE  CREATION  BOUNDARY  COMPRISED 
C  OF  SURFACES  FROM  TWO  MESHS 

C  DEFINITION  OF  GLOBAL  -PARAMETERS 
SGLOBAL 

FRINGE  °  2, 

$END 


C  MESB  DEFINITIONS 

$MESH  NAME  =>  'CARTESIAN  GRID', 

LINK  =  , 

SEND 

SMESH  NAME  -  'CYLINDER  GRID  1', 

LINK  -  'CYLINDER  GRID  2',  'CARTESIAN  GRID', 
XO  ~  20.0, 

YO  -  20.0, 

SEND 

SMESH  NAME  »  'CYLINDER  GRID  2', 

LINK  -  'CYLINDER  GRID  1',  'CARTESIAN  GRID', 
XO  =>  20.0, 

YO  -  20.0, 

SEND 


c  BOUNDARY  DEFINITIONS  CONSIST  OF  ONE  HOLE 
C  CREATION  BOUNDARY  AND  TWO  OUTER  BOUNDARIES. 

SBOUNDARY  name  -  'CYLINDER  BOLE  BOUNDARY', 
ISPARTOF  =  ' CYLINDER  GRID  1'  , 

MHOLEIH  "  ' CARTESIAN  GRID' , 

SEND 


SBOUNDARY  NAME  -  'CYLINDER  1  OUTER  BOUNDARY', 

ISPARTOF  -  'CYLINDER  GRID  1', 

SEND 

SBOUNDARY  NAME  -  'CYLINDER  2  OUTER  BOUNDARY', 

ISPARTOF  *>  'CYLINDER  GRID  2', 

SEND 

a.  Page  1 

Figure  39.  Input  for  hole  creation  using  surfaces  from  two  separate  meshes. 
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C  SURFACE  DEFINITIONS 

C  TWO  SURFACES  ARE  NEEDED  TO  DESCRIBE ‘ 

C  THE  CYLINDER  HOLE  BOUNDARY.  ONE  SURFACE 
C  RESPECTIVELY  FROM  EACH  OF  THE  CYLINDER 
C  HESHS. 


$ SURFACE  ISPARTOF  -  'CYLINDER  HOLE  BOUNDARY', 
ORANGE  -  21,21, 

KBANGE  -  1,41, 

LRAHGE  -  1,5, 

NVOUT  -  '+J', 

SEND 

$ SURFACE  ISPARTOF  "  'CYLINDER  HOLE  BOUNDARY', 
MESH  NAME  b  'CYLINDER  GRID  2', 

JRAHGS  -  21,21, 

KRANGE  -  1,41, 

LRANGE  -  1,5, 

NVOUT  -  ' + J'  , 

SEND 


C  THE  LAST  SIX  SURFACES  ARE  THE  OUTER  BOUNDARIES 
C  OF  THE  CYLINDER  GRIDS.  NO  NORMAL  VECTORS  ARE 
C  SPECIFIED  FOR  SURFACES  WHICH  ARE  CONNECTED  TO 
C  OUTER  BOUNDARIES  THROUGH  THE  ISPARTOF  USER  INPUT. 

$ SURFACE  ISPARTOF  -  'CYLINDER  1  OUTER  BOUNDARY', 
ORANGE  o  40,41, 

KRANGE  -  1,41, 

LRAHGE  1,5, 

SEND 

S SURFACE  ISPARTOF  -  'CYLINDER  1  OUTER  BOUNDARY', 
ORANGE  -  1,41, 

KRANGE  -  1,2, 

LRAHGE  -  1,5, 

SEND 

SSURFACE  ISPARTOF  ■  'CYLINDER  1  OUTER  BOUNDARY', 
ORANGE  =  1,41, 

KRANGE  =  40,41, 

LRANGE  -  1,5, 

SEND 

SSURFACE  ISPARTOF  -  'CYLINDER  2  OUTER  BOUNDARY', 
ORANGE  -  40,41, 

KRANGE  «  1,41, 

LRANGE  -  1,5, 

SEND 

SSURFACE  ISPARTOF  ■  'CYLINDER  2  OUTER  BOUNDARY', 
ORANGE  °  1,41, 

KRANGE  -  1,2, 

LRANGE  o  1,5, 

SEND 

SSURFACE  ISPARTOF  »  ' CYLINDER  2  OUTER  BOUNDARY' , 
ORANGE  -  1,41, 

KRANGE  o  40,41, 

LRANGE  =  1,5, 

Send 

C  END  OF  PEGSUS  4.0  INPUT  EXAMPLE  4.1.4 

b.  Page  2 

Figure  39.  Concluded. 
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Figure  40.  Diagnostic  maps  for  hole  creation  using  surfaces  from 
two  separate  meshes. 
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Figure  40.  Continued. 
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Figure  40.  Concluded. 
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PEGSUS  4.0  INPUT 
CARTESIAN  GRID  /  CYLINDER 
EXAMPLE  4.2.1  -  ADDITION 


C  MESH  DEFINITIONS 

C  THE  NEW  CARTESIAN  MESH  IS  THE  ONLY  MESH  HERE. 

C  THIS  MESH  WILL  BE  ADDED  TO  THE  CONFIGURATION 
C  OF  EXAMPLE  4.1.1  . 

C  THE  'ADDLINK'  ENTRY  WILL  ADD  THE  NEW  MESH  TO 
C  THE  LINK  LISTS  OF  THE  OTHER  GRIDS. 

C  THE  'ADDHOLE'  ENTRY  WILL  ADD  THE  NEW  MESH  TO 
C  THE  MAKE  HOLE  IN  LISTS  OF  THE  HOLE  CREATION 
C  BOUNDARY. 

$MESH  NAME  -  'NEW  CARTESIAN  GRID', 

MODE  =  'ADD' , 

LINK  =  , 

ADDLINK  •=  'CARTESIAN  GRID' ,' CYLINDER  GRID', 

ADDHOLE  -  'CYLINDER  HOLE  BOUNDARYl' ,' CYLINDER  HOLE  BOONDARY2', 
XO  =  25.0, 

SEND 


C  BOUNDARY  DEFINITIONS  CONSIST  OF  TWO  OUTER  BOUNDARIES. 

C  AN  OUTER  BOUNDARY  WILL  BE  ADDED  TO  THE  ORIGINAL  CARTESIAN 
C  MESH  AND  AND  OUTER  BOUNDARY  WILL  BE  DEFINED  IN  THE 
C  NEW  CARTESIAN  MESH. 


SBOUNDARY  NAME  =  ' CARTESIAN  OUTER  BOUNDARY' , 
ISPARTOF  ■=  'CARTESIAN  GRID', 

MODE  -  ' ADD' , 

Send 

$ BOUNDARY  NAME  =  'NEW  CARTESIAN  OUTER  BOUNDARY', 
ISPARTOF  =  'NEW  CARTESIAN  GRID', 

SEND 

C  SURFACE  DEFINITIONS 


C  THE  LAST  TWO  SURFACES  ARE  THE  OUTER  BOUNDARIES 
C  OF  THE  CARTESIAN  GRIDS.  NO  NORMAL  VECTORS  ARE 
C  SPECIFIED  FOR  SURFACES  WHICH  ARE  CONNECTED  TO 
C  OUTER  BOUNDARIES  THROUGH  THE  ISPARTOF  USER  INPUT. 


SSURFACE  ISPARTOF  =  'CARTESIAN  OUTER  BOUNDARY', 
JRANGE  =  40,41, 

KRANGE  =1,41, 

LRANGE  “1,5, 

SEND 

SSURFACE  ISPARTOF  -  'NEW  CARTESIAN  OUTER  BOUNDARY', 
JRANGE  =1,2, 

KRANGE  =  1,41, 

LRANGE  =  1,5, 

SEND 


C  END  OF  PEGSUS  4.0  INPUT  FOR  EXAMPLE  4.2.1 


Figure  42.  Input  for  addition  of  mesh  and  boundary. 
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a.  Page  1 

Figare  43.  Diagnostic  maps  for  addition  of  a  mesh  and  boundary. 
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'  FOR 
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Figure  43.  Continued. 
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HAP  FOR  HEH  CARTESIAH  GRID 


PLANE  Lb  X 
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Figure  43.  Concluded. 
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C  PEGSUS  4.0  INPGT 

C  CARTESIAN  GRID  I  CYLINDER 

C  EXAMPLE  4.2.2  -  SnSTRACTlOM 


C  MESH  DEFINITIONS 

C  THE  CYLINDER  MESH  WILL  BE  SUBTRACTED 
C  FROM  THE  MULTIPLE -MESH  CONFIGURATION 
C  USED  IN  EXAMPLE  4,2.1  . 


$HESH  NAME  -  'CYLINDER  GRID', 
MODE  =  'SUB' 

SEND 


C  END  OF  PEGSUS  4.0  INPUT  FOR  EXAMPLE  4.2.2 

Figure  44.  Input  for  subtraction  of  a  mesh  and  boundary. 
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DIAGHOSTIC  MAPS 
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Figure  45.  Diagnostic  maps  for  subtraction  of  a  mesh  and  boundary. 
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HAP  FOR  NEW  CARTESIAN  GRID 
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Figure  45.  Concluded, 
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APPENDIX  A 
GLOSSARY  OF  TERMS 

Blanked  Point  A  mesh  point  completely  out  of  the  computational  domain 

(hole  point),  updated  by  interpolation  (interpolation  boundary 
point),  or  updated  by  boundary  conditions  (interior  blanked 
point). 

Boundary  Collection  of  one  or  more  surfaces  and/or  one  or  more  boxes. 

Box  Group  of  points  defined  by  X,  Y,  and  Z  coordinate  ranges. 

Collections  of  one  or  more  boxes  comprise  a  boundary. 

Candidate  Point  Boundary  point  for  which  an  interpolation  stencil  in  another 

mesh  is  to  be  found. 

Coincident  Points  Points  in  different  meshes  which  are  closer  than  a  predefined 

tolerance. 

Composite  Mesh  Single  file  containing  all  meshes  in  a  multiple-mesh 

configuration. 

Computational  Domain  Mesh  defined  in  the  coordinate  system  of  the  rectangular  box 

into  which  the  physical  domain  has  been  mapped. 

Donor  Mesh  Mesh  that  sends  information  by  interpolation  to  another  mesh. 

Field  Point  Point  within  computational  domain  (IBLANK=1)  that  is 

updated  by  the  flow  solver  or  boundary  conditions. 

Fringe  Point 

Hole 


Hole  Boundary  Collection  of  hole  boundary  points  (fringe  points)  surrounding 

a  hole. 

Hole  Boundary  Point  Point  on  a  hole  boundary,  by  definition  a  point  between  a  hole 

point  and  field  point. 


A  hole  boundary  point,  q.v. 

Collection  of  hole  points,  i.e.,  points  in  a  mesh  that  have  been 
blanked. 
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Hole  Creation  Boundary 

Hole  Point 

Interpolation  Boundary 
Point 

Interpolation  File 

Interpolation  Stencil 

Level  Surface 
Orphan  Point 

Outer  Boundary 
Outer  Boundary  Point 

Overlap 

Physical  Domain 

Recipient  Mesh 
Region 


Collection  of  surfaces  in  one  or  more  meshes  that  creates  a 
hole  in  another  mesh. 

A  mesh  point  that  has  been  specified  to  be  within  a  hole.  This 
point  is  considered  to  be  out  of  the  computational  domain 
(IBLANK=0). 

A  hole  boundary  point  or  outer  boundary  point  that  is 
interpolated  from  another  mesh  (IBLANK=0). 

File  consisting  of  (1)  boundary  points,  (2)  pointers  to 
interpolation  stencils,  and  (3)  interpolation  coefficients  for  each 
mesh.  The  interpolation  file  is  output  by  PEGSUS  and  is  input 
to  the  flow  solver. 

The  eight  mesh  points  that  surround  a  candidate  boundary 
point. 

A  constant  J,  K,  or  L  plane  of  a  mesh. 

Boundary  point  in  which  a  valid  interpolation  can  not  be  found 
in  any  LINKed  mesh. 

Collection  of  outer  boundary  points. 

A  point  in  a  mesh  that  has  been  specified  directly  by  a  surface 
to  be  interpolated  (typically  this  point  is  on  the  outer  boundary 
of  a  mesh,  not  the  extreme  boundary  of  the  composite  mesh). 

The  field  points  lying  between  interpolation  boundary  points 
contained  in  adjacent  meshes. 

Mesh  defined  in  a  coordinate  system  with  shape  and  topology 
that  is  natural  for  the  complex  geometry  of  the  configuration 
being  represented. 

Mesh  that  receives  information  from  another  mesh. 
Collection  of  one  or  more  volumes. 
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Surface 


Valid  Interpolation 


Volume 


Part  of  a  level  surface  of  a  mesh.  Collections  of  one  or  more 
surfaces  comprise  a  boundary. 

An  interpolation  where  all  eight  mesh  points  of  the 
interpolation  stencil  are  field  points  (not  hole  or  interpolation 
boundary  points). 

Group  of  points  defined  by  J,  K,  and  L  index  ranges. 
Collections  of  one  or  more  volumes  comprise  a  region. 
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APPENDIX  B 
PROGRAM  VARIABLES 

The  following  is  a  list  of  the  variables  used  within  the  program.  Listed  are  the  variable 
name,  type,  and  description.  The  CHARACTER  variable  type  is  abbreviated  with  the  letters 
CH  (i.e.,  CH*1  is  CHARACTER*  1).  The  variables  are  grouped  within  COMMON  blocks 
as  they  are  found  within  the  program. 

Name  Type  Deaciiptioii 

COMMON/MESHPARM/ 


NMESH 

INTEGER 

number  of  meshes 

NAME(ICHAR.MDIM) 

CH*1 

mesh  name 

LINK(ICHAR,MDIM,LNDIM) 

CH*1 

name  of  meshes  linked  to  this  mesh 

NLINK(MDIM) 

INTEGER 

number  of  meshes  linked  to  this  mesh 

LINKALL(MDIM) 

LOGICAL 

=  .TRUE,  if  all  mesh  links  are  to  be 
connected  to  this  mesh 

MMODE(MDIM) 

CH*3 

=  ‘ADD’  or  ‘SUB’ 

IFRINGE(MDIM) 

INTEGER 

number  of  fringe  points 
(=  1  or  2) 

QUALITY(3,MDIM) 

REAL 

values  of  interpolation  quality  for  mesh 

EPS(MDIM) 

REAL 

tolerance  of  interpolation  coefficient 

XINCLD(2,MDIM) 

REAL 

range  of  X  values  to  be  included  for 
interpolation 

YINCLD{2,MDIM) 

REAL 

range  of  Y  values  to  be  included  for 
interpolation 

ZINCLD(2.MDIM) 

REAL 

range  of  Z  values  to  be  included  for 
interpolation 

JINCLD(2,MDIM) 

INTEGER 

range  of  J  indices  to  be  included  for 
interpolation 

KINCLD(2,MDIM) 

INTEGER 

range  of  K  indices  to  be  included  for 
interpolation 

LINCLD(2,MDIM) 

INTEGER 

range  of  L  indices  to  be  included  for 
interpolation 

INCLD{2,MDIM) 

LOGICAL 

specifies  whether  include  parameter 
specified  for  mesh 

PHANTOM(MDIM) 

LOGICAL 

=  .TRUE,  if  phantom  mesh 
=  .FALSE,  if  not  a  phantom  mesh 
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Name 

Type 

Description 

COMMON/BOUNPARM/ 

BNAIVIE{1CHAR,NBDIM) 

CHM 

boundary  name 

NBOUN 

INT 

number  of  boundaries 

IPARTMB(ICHAR,NBDIM) 

CH*1 

name  of  mesh  th£^  this  boundary  is  part 
of 

=  ‘ADD’  or  ‘SUB’ 

BMODE(NBDIM) 

CH*3 

BTYPE(NBDIM) 

CH*5 

=  ‘HOLE’  or  ‘OUTER’ 

MHOLEIN(ICHAR.NBDIM. 

CH*1 

name  of  mesh  that  this  boundary 

MHDIM) 

makes  a  hoie  in 

NHOLEIN(NBDIM) 

INTEGER 

number  of  meshes  that  this  boundary 
makes  holes  in 

CLOSED(NBDIM) 

LOGICAL 

=  .TRUE. (default)  if  boundary 
described  by  surfaces  is  closed  (used 
only  for  hole  boundaries  described  by 
surfaces) 

COMMON/SURFPARM/ 


NSURF 

INTEGER 

total  number  of  surfaces 

JRANGE(2,NSDIM) 

INTEGER 

range  of  J  indices  of  surface 

KRANGE(2,NSDIM) 

INTEGER 

range  of  K  indices  of  surface 

LRANGE(2,NSDIM) 

INTEGER 

range  of  L  indices  of  surface 

NVOUT(NSDIM) 

CH*2 

outward  normal  vector  of  surface 
(required  for  hole  boundaries) 

IPARTBS(ICHAR,NSDIM) 

CH"! 

name  of  boundary  which  this  surface 
is  a  part  of 

ISMESH(ICHAR,NSDIM) 

CH"! 

name  of  mesh  that  this  surface  is 
created  from 

COMMON/BOXPARM/ 

NBOX 

INTEGER 

number  of  boxes 

XRANGE(2,NXDIM) 

REAL 

minimum  and  maximum  X  values  of 

box 

YRANGE(2,NXDIM) 

REAL 

minimum  and  maximum  Y  values  of 

box 
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Name  Type  Description 

COMMON/BOXPARM/  (Continued) 

ZRANGE(2,NXDIM)  REAL  minimum  and  maximum  Z  values  of 

box 

IPARTBX(ICHAR,NXDIM)  CH*1  name  of  boundary  that  this  box  is  a 

part  of 

COMMON/REGNPARM/ 

RNAME{ICHAR,NRDIM) 

IPARTMR(ICHAR,NRDIM) 

RMODEfNRDIM) 

RTYPE(NRDIM) 

NREGN 


COMMON/VOLPARM/ 

NVOL 

INTEGER 

number  of  volumes 

JRANGEV(2,NVDIM) 

INTEGER 

range  of  J  indices  defining  the  volume 

KRANGEV(2,NVDIM) 

INTEGER 

range  of  K  indices  defining  the  volume 

LRANGEV(2,NVDIM) 

INTEGER 

range  of  L  indices  defining  the  volume 

IPARTR(ICHAR,NVDIM) 

CH*1 

name  of  region  number  that  this 
volume  is  a  part  of 

COMMON/MESHLINK/ 

MLINK(MDIM,LNDIM) 

INTEGER 

mesh  number  that  is  linked 

COMMON/BOUNLINK/ 

IHNUM(MDIM) 

INTEGER 

number  of  holes  made  in  mesh 

IHOLE(MDIM,NBDIM) 

INTEGER 

boundaries  that  make  holes  in  mesh 

IONUM(MDIM) 

INTEGER 

number  of  outer  boundaries  in  mesh 

IOUTER(MDIM,NBDIM) 

INTEGER 

outer  boundaries  in  mesh 

CH*1  region  name 

CH*  1  name  of  mesh  that  this  region  Is  a  part 

of 

CH*3  =  'ADD*  or  ‘SUB’ 

CH*4  =  ‘HOLE’  or  ‘INTR’ 

INTEGER  number  of  regions 
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Name 

Type 

Description 

COMMON/SURFLINK/ 

NUSURF(NBD1M) 

INTEGER 

number  of  surfaces  that  make  up 
boundary 

IDSURF(NBDIM,NSDIM) 

INTEGER 

surfaces  that  make  up  boundary 

COMMON/BOXLINK/ 

NUBOX(NBDIM) 

INTEGER 

number  of  boxes  that  make  up 
boundary 

IDBOX(NBDIM,NXDIM) 

INTEGER 

boxes  that  make  up  boundary 

COMMON/REGNUNK/ 

IRNUM(MD1M) 

INTEGER 

number  of  regions  in  mesh 

IREGN(MDIM,NRDIM) 

INTEGER 

regions  in  mesh 

COMMON/VOLUNK/ 

NUVOUNRDIM) 

INTEGER 

number  of  volumes  that  make  up 
region 

IDVOL(NRDIM,NVDIM) 

INTEGER 

volumes  that  make  up  region 

COMMON/MESHCNTL/ 

NEW(MDIM) 

LOGICAL 

=  .TRUE,  for  initial  run  or  if  mesh  is 
added 

ADD(MDIM) 

LOGICAL 

=  .TRUE,  when  boundary  or  region  is 
added  to  mesh 

SUB(MD1M) 

LOGICAL 

=  .TRUE,  when  a  boundary  or  region 
is  subtracted  from  mesh 

COMMON/POINTERS/ 

IPIP(MDIM,MDIM) 

INTEGER 

array  pointer  for  interpolation  points 

IPMS(MDIM) 

INTEGER 

array  pointer  for  grids  (based  on  the 
grid  sizes) 
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Name  Type  Description 


COMMON/POINTERS/ 


(Continued) 

IPSF(NSDIM) 

INTEGER 

array  pointer  for  surfaces 

IPIN(MDIM) 

INTEGER 

array  pointer  for  points  to  be 
interpolated 

COMMON/GLOBAL/ 

IFRNGGL 

INTEGER 

global  value  for  FRINGE  no. 

QUALGL 

REAL 

global  value  for  QUALITY 

EPSGL 

REAL 

global  value  for  EPS 

XINCGL(2) 

REAL 

global  value  for  XINCLD 

YINCGL(2) 

REAL 

global  value  for  YINCLD 

ZINCGL(2) 

REAL 

global  value  for  ZINCLD 

COMMON/RSTART/ 

RESTART 

LOGICAL 

=  .TRUE,  if  PEGSUS  restart  file  read 
in 

number  of  meshes  in  restart 

NMESHR 

INTEGER 

NBOUNR 

INTEGER 

number  of  boundaries  in  restart 

NSURFR 

INTEGER 

number  of  surfaces  in  restart 

NBOXR 

INTEGER 

number  of  boxes  in  restart 

NREGNR 

INTEGER 

number  of  regions  in  restart 

NVOLR 

INTEGER 

number  of  volumes  in  restart 

COMMON/RFLAGS/ 

ADDM 

LOGICAL 

=  .TRUE,  if  mesh  added 

ADDB 

LOGICAL 

=  .TRUE,  if  boundary  added 

ADDR 

LOGICAL 

=  .TRUE,  if  region  added 

SUBM 

LOGICAL 

=  .TRUE,  if  mesh  subtracted 

SUBB 

LOGICAL 

=  .TRUE,  if  boundary  subtracted 

SUBR 

LOGICAL 

=  .TRUE,  if  region  subtracted 
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Name 

Type 

Description 

COMMON/SllKF/ 

ISURF(MSSMAX) 

INTEGER 

absolute  index  of  surface  point  in  mesh 

VNX(MSSMAX) 

REAL 

magnitude  of  normal  in  X  direction  of 

VNY(MSSMAX) 

REAL 

surface  point 

magnitude  of  normal  in  Y  direction  of 

VNZ(MSSMAX) 

REAL 

surface  point 

magnitude  of  normal  in  Z  direction  of 

XS(MSSMAX) 

REAL 

surface  point 

X  coordinate  of  surface  point 

YS(MSSMAX) 

REAL 

Y  coordinate  of  surface  point 

ZS(MSSMAX} 

REAL 

Z  coordinate  of  surface  point 

COMMON/ERROR/ 

NUMERR 

INTEGER 

number  of  errors  found 

COMMON/NMUST/ 

AGLOBAL 

CH*2000 

global  namelist  input 

AMESH(MDIM) 

CH*2000 

mesh  namelist  input 

ABOUND(NBDIM) 

,  CH*2000 

boundary  namelist  input 

ASURF(NSDIM) 

CH*2000 

surface  namelist  input 

ABOX(NXDIM) 

CH*2000 

box  namelist  input 

AREGN(NRDIM) 

CH*2000 

region  namelist  input 

AVOL(NVDlM) 

CH*2000 

volume  namelist  input 

COMMON/IWORK/ 

IWORKl(MLEN) 

INTEGER 

integer  work  array 

IWORK2(MLEN) 

INTEGER 

integer  work  array 

IWORK3(MLEN) 

INTEGER 

integer  work  array 
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Name  Type  Description 

COMMON/INVAL/ 


NMESHI 

INTEGER 

number  of  meshes  input 

NBOUNI 

INTEGER  , 

number  of  boundaries  input 

NSURFI 

INTEGER 

number  of  surfaces  input 

NBOXl 

INTEGER 

number  of  boxes  input 

NREGNI 

INTEGER 

number  of  regions  input 

NVOLI 

INTEGER 

number  of  volumes  input 

COMMON/QUALTOT/ 

NUMQUAL(ll) 

INTEGER 

number  of  interpolated  boundary 
points  found  for  different  ranges  of 
minimum  QUALITY 

COMMON/MESH/ 

X(MLEMAX) 

REAL 

X  coordinate  of  grid 

Y(MLEMAX) 

REAL 

Y  coordinate  of  grid 

Z(MLEMAX) 

REAL 

Z  coordinate  of  grid 

COMMON/BOXG/ 

XMAXG 

INTEGER 

maximum  X  of  domain 

XMING 

INTEGER 

minimum  X  of  domain 

YMAXG 

INTEGER 

maximum  Y  of  domain 

YMING 

INTEGER 

minimum  Y  of  domain 

ZMAXG 

INTEGER 

maximum  Z  of  domain 

ZMING 

INTEGER 

minimum  Z  of  domain 

COMMON/BOXM/ 

XMAXM{MDIM) 

INTEGER 

maximum  X  of  mesh 

XMINM(MDIM) 

INTEGER 

minimum  X  of  mesh 

YMAXM(MDIM) 

INTEGER 

maximum  Y  of  mesh 

YMINM(MDIM) 

INTEGER 

minimum  Y  of  mesh 

ZMAXM(MDIM) 

INTEGER 

maximum  Z  of  mesh 

ZM1NM(MDIM} 

INTEGER 

minimum  Z  of  mesh 
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Name 

Type 

Desciiptloa 

COMMON/BOXH/ 

XMAXH(NBDIM} 

INTEGER 

maximum  X  of  hole  boundary 

XMINH(NBDIM) 

INTEGER 

minimum  X  of  hole  boundary 

YMAXH(NBDIM) 

INTEGER 

maximum  Y  of  hole  boundary 

YMINH(NBDIM) 

INTEGER 

minimum  Y  of  hole  boundary 

ZMAXH(NBDIM) 

INTEGER 

maximum  Z  of  hole  boundary 

ZMINH(NBDIM) 

INTEGER 

minimum  Z  of  hole  boundary 

COMMON/MESHVAL/ 

MJMAX(MDIM) 

INTEGER 

maximum  J  for  each  mesh 

MKMAX(MDIM) 

INTEGER 

maximum  K  for  each  mesh 

MLMAX(MDIM) 

INTEGER 

maximum  L  for  each  mesh 

X0(MDIM) 

REAL 

translation  value  of  X 

YO(MDIM) 

REAL 

translation  value  of  Y 

ZO(MDIM) 

REAL 

translation  value  of  Z 

XR(MDIM) 

REAL 

X  coordinate  of  rotation  point 

YR(MDIM) 

REAL 

Y  coordinate  of  rotation  point 

ZR(MDIM) 

REAL 

Z  coordinate  of  rotation  point 

ALPHA(MDIM) 

REAL 

Euler  angle  for  rotation  about  Z  axis 

BETA(MDIM) 

REAL 

Euler  angle  for  rotation  about  Y  axis 

GAMA(MDIM) 

REAL 

Euler  angle  for  rotation  about  X  axis 

SCALE(MDIM) 

REAL 

scaling  value 

DELTA(MDIM) 

REAL 

minimum  diagonal  of  a  cell  in  each 
mesh  multiplied  by  0.0001 

COMMON/INTRPTR/ 

NUINTR(MDIM.MDIM) 

INTEGER 

number  of  boundary  (or  stencil)  points 
for  each  mesh-to-mesh  interface 

COMMON/ILIST/ 

NULIST(MDIM) 

INTEGER 

number  of  points  for  interpolation 

INLIST(MIMAXl) 

INTEGER 

points  for  interpolation 
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Name 

Type 

Description 

COMMON/PTRS/ 

IIPNTS(MDIM) 

INTEGER 

number  of  interpolation  stencils  in  each 
mesh 

IBPNTS(MDIM) 

INTEGER 

number  of  interpolation  boundary 
points  in  each  mesh 

IIEPTR(MDIM) 

INTEGER 

ending  pointer  for  the  interpolation 
stencil  list  for  each  mesh 

IISPTR(MDIM) 

INTEGER 

starting  pointer  for  the  interpolation 
stencil  list  for  each  mesh 

COMMON/INTERP/ 

DXINT(MIMAX2) 

REAL 

X  direction  interpolation  coefficient 

DYINT(MIMAX2) 

REAL 

Y  direction  interpolation  coefficient 

DZINT(MIMAX2) 

REAL 

Z  direction  interpolation  coefficient 

JI(MIMAX2) 

INTEGER 

J  index  of  stencil  point 

K1(MIMAX2) 

INTEGER 

K  index  of  stencil  point 

L1(M1MAX2) 

INTEGER 

L  index  of  stencil  point 

JB(MIMAX2) 

INTEGER 

J  index  of  boundary  point 

KB(M[MAX2) 

INTEGER 

K  index  of  boundary  point 

LB(IVIIMAX2) 

INTEGER 

L  index  of  boundary  point 

IBC{MILEN) 

INTEGER 

interpolation  boundary  point  pointer 

COMMON/BLANK/ 

IBLANK{MLEMAX) 

INTEGER 

iblank  array 
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Figure  C-1  depicts  the  hierarchical  structure  of  PEGSUS  4.0.  Each  box  (with  the  exception 
of  the  PEGSUS  4.0  box,  which  is  the  main  program)  represents  a  subroutine  or  function 
within  PEGSUS  4.0.  An  asterisk(*)  within  a  box  denotes  that  the  subroutine  calls  other 
subroutines  or  functions  which  are  shown  on  subsequent  pages. 
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a.  Program  structure,  page  1 
Figure  C>1.  Program  structure. 


lEDC-TH-SI-B 


c.  Progrtun  structure,  page  3 
Figure  C>1.  Continued. 


d.  Program  structure,  page  4 
Figure  C-1.  Continued. 
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e.  Program  stnicture,  page  5 
Figure  C*l.  Continued. 
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g.  Program  structure,  page  7 
Figure  C*l.  Continued. 
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SETLINK 


h.  Program  structure,  page  8 
Figure  C-1.  Continued. 


1.  Program  strnctiire,  page  9 
Figure  C-1.  Continued. 
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Figure  C-1.  Continued. 
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AEDC'TR-91'a 


ADDHOLES 

ADDLINKS 

BGATHER 

BOUNDNO 

BOUNMODE 

CALCINT 

CHKBPTS 

CHKMBOX 

CHKPARM 

CHKRANGE 

CHKSTEN 


CHKSURF 

CKABOUN 


APPENDIX  D 

SUBROUTINE  DESCRIPTION 

Processes  user  input  that  specifies  added  hole  boundaries. 

Processes  user  input  that  specifies  added  links. 

Performs  gather-type  function.  Target  and  source  arrays  may  have 
different  dimensions. 

Given  a  boundary  name,  returns  the  boundary  number. 

For  a  restart  condition,  processes  the  user  input  based  on  the 
boundary  mode  parameter. 

From  a  list  of  points  requiring  interpolation,  finds  interpolation 
coefficients. 

For  a  restart  condition,  checks  previous  boundary  points  to 
determine  if  they  are  still  valid  in  current  run. 

Given  a  point  and  a  mesh  box  (minimum/maximum  region), 
determines  if  point  is  inside  the  mesh  box. 

For  restart  condition,  checks  current  parameters  against  previous 
parameters. 

Checks  the  J,  K,  and  L  index  ranges  for  surfaces  and  volumes  to 
ensure  that  they  are  between  1  and  the  maximum  values  for  their 
respective  mesh. 

For  restart  condition,  checks  ail  stencils  to  determine  if  previous 
interpolations  are  still  vajid.  If  not,  they  are  removed  from  the  list 
of  points  that  have  been  interpolated  and  are  added  to  the  list  of 
points  that  need  interpolation. 

Checks  surface  inputs  to  determine  if  connected  to  existing  mesh. 
Checks  the  key  words  for  user  input  SBOUNDARY. 
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CKABOX 

CKAGLOB 

CBCAMESH 

CKAREGN 

CKASURF 

CKAVOL 

CKBNAME 

CKNAME 

CKRNAME 

CNTLBOUN 

CNTLMESH 

CNTLREGN 

DEFADD 

DEFBOTH 

DEFBOUN 

DEFNEW 

DEFSUB 

DIAGOUT 


Checks  the  key  words  for  user  input  $BOX. 

Checks  the  key  words  for  user  input  SGLOBAL. 

Checks  the  key  words  for  user  input  $MESH. 

Checks  the  key  words  for  user  input  $REGION. 

Checks  the  key  words  for  user  input  SSURFACE. 

Checks  the  key  words  for  user  input  $VOLUME. 

Given  a  boundary  name,  checks  all  previous  boundary  names  for 
duplication. 

Given  a  mesh  name,  checks  all  previous  mesh  names  for  duplication. 

Given  a  region  name,  checks  all  previous  region  names  for 
duplication. 

Sets  control  flags  based  on  boundaries. 

Sets  control  flags  based  on  meshes. 

Sets  control  flags  based  on  regions. 

Defines  all  interpolated  boundary  points  for  meshes  with  addition. 

Defines  all  interpolation  boundary  points  for  meshes  with  both 
subtraction  and  addition. 

Defines  all  interpolation  boundary  points. 

Defines  all  interpolation  boundary  points  for  new  meshes. 

Defines  all  interpolation  boundary  points  for  meshes  with 
subtraction. 

Writes  out  the  diagnostic  maps  for  a  given  mesh. 
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EXCLPTS 

FFRNGPTS 

FINDPTS 

FNDBOUN 

FNDHOLE 

FNDNEAR 

FNDOUTR 

FNDREGN 

FRINGE! 

FRINGE2 

GETCOOR 

GETINPTS 

INITIA 

INTRP 

ISINSIDE 

ITOJKL 

JKLTOI 


Excludes  specified  interpolation  boundary  points  based  on  user 
inputs  XINCLUPE,  YINCLUDE,  ZINCLUDE,  JINCLUDE. 
KINCLUDE,  and  LINCLUDE. 

Finds  fringe  points  surrounding  hole  points. 

Finds  all  boundary  points  that  require  interpolation. 

Finds  the  hole  points  defined  by  a  boundary. 

Finds  all  hole  points  within  a  mesh. 

Given  a  point,  finds  the  nearest  point  in  the  given  mesh  subset. 

Sets  IBLANK  array  for  outer  boundary  points. 

Finds  the  hole  points  as  defined  by  a  region. 

Finds  a  single  fringe  surrounding  all  hole  points. 

Finds  a  double  fringe  surrounding  all  hole  points. 

From  a  list  of  indices,  gets  the  X,  Y,  and  Z  coordinates. 

Determines  which  boundary  points  identiEed  in  DEFBOUN  require 
interpolation.  Boundary  and  stencil  points  are  corrected  for 
configuration  changes. 

Initializes  all  COMMON  blocks  and  opens  I/O  files. 

Computes  interpolation  coefficients  for  trilinear  interpolation  using 
Newton’s  method.  Performs  secondary  function  of  “stencil¬ 
walking”  during  search  for  interpolation  stencil. 

Returns  .TRUE,  if  a  mesh  point  is  inside  a  surface;  .FALSE, 
otherwise. 

Calculates  J,  K,  and  L  indices  from  a  given  absolute  index  I. 
Calculates  an  absolute  index  I  from  given  J,  K,  and  L  indices. 
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LINKBOUN 

LINKMESH 

LINKREGN 

MESHMODE 

MESHNO 

NEWHOLE 

NEWMESH 

NEWTON 

NORMAL 


PEGSUS40 

PROCBOUN 

PROCBPTS 

PROCGLOB 

PROCIBLK 

PROCIN 

PROCINCL 

PROCMESH 

PROCNEW 


Sets  the  boundary  links  based  on  user  inputs. 

Sets  the  mesh  links  based  on  user  inputs. 

Sets  the  region  links  based  on  user  inputs. 

For  a  restart  condition,  processes  the  user  inputs  based  on  the  mesh 
mode  parameter. 

Given  a  mesh  name,  returns  a  mesh  number. 

Sets  up  surface  normals  and  the  hole  box  for  new  holes. 

Reads  and  processes  new  meshes. 

Performs  trilinear  interpolation  using  Newton’s  method. 

Calculates  norma!  vectors  at  all  points  on  a  surface.  The  surface 
is  assumed  to  be  defined  in  a  right-handed  coordinate  system.  A 
left-handed  system  will  result  in  surface  normals  facing  inward, 
which  will  cause  holes  to  be  generated  incorrectly. 

Main  calling  routine. 

Processes  the  boundary-related  values  of  the  user  inputs'. 

For  a  restart  condition,  reads  and  processes  previous  interpolation 
boundary  points. 

Processes  global  values  of  the  user  inputs. 

For  a  restart  condition,  reads  previous  IBLANK  array. 

Reads  and  processes  all  user  inputs,  meshes  and  restart  files. 

Sets  the  include  values  XINCLUDE,  etc. 

Processes  the  mesh  values  of  the  user  inputs. 

Processes  new  meshes  and  holes. 
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PROCOUT  ' 
PROCREGN 
PROCRST 
PROCSTEN 

RDINBOUN 

RDINBOX 

RDINMESH 

R0INREGN 

RDINSURF 

RDINVOL 

READBOUN 

READIN 

READIPT 

READMESH 

READNMC 

READNMC2 

READNMC3 

READNMC4 

READNMCM 


Processes  and  writes  results  to  appropriate  files. 

Processes  the  region-related  values  of  the  user  inputs. 

Reads  and  processes  the  restart  values. 

Reads  and  processes  the  previous  interpolation  stencils  and 
coefficients. 

Reads  the  boundary  user  inputs. 

Reads  the  box  user  inputs. 

Reads  the  mesh  user  inputs. 

Reads  the  region  user  inputs. 

Reads  the  surface  user  inputs. 

Reads  the  volume  user  inputs. 

Reads  the  boundary-related  values  from  PEGSUS  restart  file. 
Reads  and  processes  all  user  inputs. 

For  a  restart  condition,  reads  pointers  for  previous  interpolation 
points. 

Reads  the' mesh  values  from  PEGSUS  restart  file. 

NAMELIST-like  read  for  a  single-variable-name  character  array. 

NAMELIST-like  read  for  a  two-character  name. 

NAMELIST-like  read  for  a  three-character  name. 

NAMELIST-like  read  for  a  four-character  name. 

NAMELIST-like  read  for  multiple-variable-length-name  character 
array  (double  dimension). 
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READNMCT 

READNMIl 

READNMI2 

READNML 

READNMRl 

READNMR2 

READNMR3 

READREGN 

READRST 

REGNMODE 

REGNNO 

RMSPCS 

ROTATE 

RSETINTR 

RSETLINK 

SCALM 

SETBPTS 

SETCNTL 


NAMELIST-like  read  for  multiple-variable-length-name  character 
array  (triple  dimension). 

NAMELIST-like  read  for  one  integer  value. 

NAMELIST-like  read  for  two  integer  values. 

NAMELIST-like  read  for  a  logical  value. 

NAMELIST-like  read  for  one  real  value. 

NAMELIST-like  read  for  two  real  values. 

NAMELIST-like  read  for  three  real  values. 

Reads  region-related  values  from  PEGSUS  restart  file. 

Reads  PEGSUS  restart  file. 

For  a  restart  condition,  processes  the  user  inputs  based  on  the  region 
mode. 

Given  a  region  name,  returns  a  region  number. 

Removes  the  spaces  from  the  NAMELIST-like  input. 

Rotates  mesh  about  defined  point  in  space. 

Resets  the  interior  blanked  points  back  to  1  before  PEGSUS 
processing. 

Resets  mesh  links  based  on  mesh  boxes. 

Scales  meshes. 

Gets  the  boundary  points  specified  by  the  IBLANK  array  and 
determines  which  values  need  to  be  interpolated.  If  mesh  is  not  new 
some  boundary  points  may  have  already  been  interpolated. 

Sets  control  flags. 
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SETDATA 

SETDEL 

SETFILES 

SETGBOX 

SETHBOX 

SETINPTS 

SETINTR 

SETINTRP 

SETLINK 

SETMBOX 

SETMESH 

SETPINTR 

SETSURF 

SHFBOUN 

SHFBOX 

SHFMESH 

SHFNUINT 

SHFPARM 

SHFREGN 


Reads  and  processes  previous  interpolation  data. 

Calculates  smallest  distance  between  points  in  a  mesh.  The  values 
are  used  to  define  how  close  two  points  from  different  meshes  must 
be  in  order  to  be  considered  coincident. 

Sets  up  direct  access  files  or  work  arrays. 

Calculates  minimum/maximum  box  of  the  domain. 

Calculates  minimum/maximum  box  of  the  hole  creation  boundary. 

Determines  which  boundary  points  require  interpolation. 

Sets  the  interior  blanked  region  equal  to  0. 

Reads  and  processes  orphan  points  from  previous  PEGSUS  run. 

Sets  the  links  based  on  user  inputs. 

Calculates  the  minimum/maximum  box  of  each  mesh. 

For  a  restart  condition,  reads  and  processes  meshes  from  previous 
composite  mesh  file. 

Sets  the  interior  blanked  region  equal  to  0  for  plotting. 

Reads  surfaces  from  PEGSUS  restart  file. 

Shifts  boundary  values  for  subtracted  boundaries. 

Shifts  box  values  for  subtracted  boundaries. 

Shifts  mesh  values  for  subtracted  meshes. 

Shifts  NUINTR  array  for  subtracted  meshes. 

Shifts  all  values  for  subtraction. 

Shifts  region  values  for  subtracted  regions. 
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SHFRST 

SHFSURF 

SHFVOL 

STOREIN 

STRIP 

SUBMESH 

SUMBOUN 

SUMGLOB 

SUMMESH 

SUMQUAL 

SUMREGN 

TESTSTN 

TRANS 

TRILIN 

WRTBOUN 

WRTBPTS 

WRTDIAG 

WRTIBLK 

WRTIDATA 


Shifts  restart  counters  for  subtracted  meshes,  boundaries,  or  regions. 

Shifts  surface  values  for  subtracted  boundaries. 

Shifts  volume  values  for  subtracted  regions. 

Reads  user  inputs  and  stores  them  in  character  arrays  for  later 
processing. 

Removes  comments  from  NAMELIST-like  input. 

Develops  two  different  subsets  of  a  mesh  for  near-point  search. 
Writes  a  summary  of  boundary-related  values. 

Writes  a  summary  of  global  values. 

Writes  a  summary  of  mesh  values. 

Writes  a  summary  of  the  quality  of  the  interpolations. 

Writes  a  summary  of  region-related  values.' 

Tests  interpolation  stencil  to  determine  if  it  is  valid. 

Translates  mesh. 

Finds  an  initial  stencil  reference  point  to  begin  search  for  valid 
interpolation  stencil. 

Writes  the  boundary-related  values  to  PEGSUS  restart  file. 

Writes  the  boundary  points  to  the  interpolation  file. 

Sets  up  array  for  diagnostic  maps  and  writes  IBLANK  array  for 
plotting. 

Writes  the  IBLANK  array  to  the  interpolation  file. 

Writes  the  interpolation  file. 
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WRTINTRP 

WRTIPT 

WRTMESH 

WRTPMESH 

WRTREGN 

WRTRST 

WRTSTEN 

WRTSUM 

WRTSURF 

XSPECS 


Writes  orphan  points  to  the  PEGSUS  restart  file. 

Writes  the  NUINTR  array  to  the  PEGSUS  restart  file. 

Writes  the  mesh  values  to  the  PEGSUS  restart  file. 

Writes  the  phantom  meshes  to  the  PEGSUS  restart  file. 

Writes  the  region-related  values  to  the  PEGSUS  restart  file. 
Writes  the  PEGSUS  restart  file. 

Writes  the  interpolation  stencils  and  coefEdents  to  the  interpolation 
nie. 

Writes  the  summary  of  user  inputs  and  meshes. 

Writes  the  surface  points  and  normals  to  the  PEGSUS  restart  file. 
Calculates  and  writes  parameter  specification  for  flow  solver. 
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APPENDIX  E 
DIRECT  ACCESS  FILES 

The  following  are  stored  as  direct  access  files  during  the  execution  of  PEGSUS  4.0  with 
INCORE =0. 

Mesh  Coordinates 
Unit 

Record  Format 
Index 

Surface  Definitions 
Unit 

Record  Format 
Index 

Blanked  Points 

Unit 

Record  Format 
Index 

Boundary  Points 
Unit 

Record  Format 


Interpolation  Stencil  Reference  Points 
Unit 

Record  Format 
Index 

Interpolation  Coefficients 

Unit  17 

Record  Format  DXINT{MSLEN),  DYINT(MSLEN), 

DZINT(MSLEN) 

Index  Mesh  ID  number 


16 

Jl(MSLEN),  KI(MSLEN).  LI(MSLEN) 
Mesh  ID  number 


15 

JB(MSLEN),  KB(MSLEN),  LB(MSLEN) 
Mesh  ID  number 


14 

IBLANK(MLEN) 
Mesh  ID  number 


13 

XS(MSLEN),  YS(MSLEN),  ZS(MSLEN), 
VNX{MSLEN),  VNY(MSLEN),  VNZ(MSLEN) 
Surface  ID  number 


12 

X(MLEN),  Y(MLEN),  Z(MLEN) 
Mesh  ID  number 
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Points  Requiring  Interpolation 
Unit 

Record  Format 
Index 


19 

INLIST(MILEN) 
Mesh  ID  number 
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APPENDIX  F 

PEGSUS  OUTPUT  FILE  STRUCTURE 

Five  output  files  are  defined:  (1)  composite  mesh  file,  (2)  interpolation  data  file,  (3) 
PEGSUS  restart  file,  (4)  IBLAl^  array  for  plotting  file,  and  (5)  orphan  list  file. 

Composite  Mesh  File 

This  file  is  written  to  file  COMPOUT.  Each  mesh  NM  is  written  separately. 
The  order  in  which  the  meshes  are  written  is  the  order  in  which  the  meshes 
are  specified  in  the  user  input  by  $MESH. 

Record  Number  Format 

1  NM,  MJMAXfNM),  MKMAX(NM).  MLMAX(NM). 
MNAME 

2  (X(I),I  =  11,12), 

(Y(I),1  =  I1,I2). 

(Z(I),I  =  11,12) 

where  11  =  1  and  12  =  MJMAX(NM)*MKMAX(NM)*MLMAX(NM) 

Interpolation  Data  File 

This  file  is  written  to  file  INTOUT.  The  data  for  each  mesh  NM  are  written 
separately. 

Record  Number  Format 

1  IBPNTS(NM),1IPNTS(NM),IIEPTR(NM), 
IISPTR(NM).MJMAX(NM),MKMAX(NM). 
MLMAX(NM) 

2  (JI(1).KI(I),LI(I),DXINT(I),DYINT(I),DZINT(I), 

1=  1,1IPNTS(NM)) 

3  (JB(1).KB(1),LB(1),IBC(I),I  =  1  ,IBPNTS(NM)) 

4  {IBLANK(I).I  =  11,12) 

where  II  =  1  and  12  =  MJMAX(NM)*MKMAX(NM)*MLMAX(NM) 
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PEGSUS  Restart  File 

This  file  is  written  to  file  RSTOUT.  The  PEGSUS  restart  file  is  divided  into  eight  parts: 
(1)  global  values.  (2)  mesh  values,  (3)  boundary-related  values,  (4)  region-related  values,  (5) 
points  requiring  interpolation  (orphan  points),  (6)  interpolation-list  pointer,  (7)  phantom 
meshes,  and  (8)  hole  creation  surfaces. 

(1)  Global  Values 

Record  Number  Format 

1  MDIM,  LNDIM,  ICHAR,  NBDIM,  MHDIM,  NSDIM, 

MSLEN,NXDIM,  NRDIM,  NVDIM,  MLEN, 
MLEMAX,  MILEN 

2  IFRNGGL,  (QUALGL(N),N  =  1,3),  EPSGL, 

(XINCGL(N),  YINCGL(N).  ZINCGL(N),  N  =  1,2), 
INCLDGL,  XMAXG,  XMING,  YMAXG,  YMING, 
ZMAXG,  ZMING 

(2)  Mesh  Values 

Record  Number  Format 

1  NMESH 

The  following  record  is  written  for  each  mesh  NM  from  1  to  NMESH. 

2  (NAME(IC,NM),IC  =  1  ,ICHAR), 

IFRINGE(NM),  (QUALITY(N.NM).N=  1,3), 
EPS(NM).  (XINCLD(N,NM), 

YINCLD(N,NM),  ZINCLD(N,NM). 

J1NCLD(N,NM),  KINCLD(N,NM), 
LINCLD{N,NM),N  =  1,2), 

INCLD{NM),  NLINK(NM), 

LINKALL(NM),  PHANTOM(NM), 
((LINK(IC,NM,NL),IC  =  1,ICHAR),NL  =  1 , 
NLINK(NM)), 

XMAXM{NM),  XMINM(NM), 

YMAXM(NM),  YMINM(NM), 

ZMAXM(NM),  ZMINM(NM), 


139 


AEDC-TR-91-8 


MJMAX(NM),  MKMAX(NM), 
MLMAX(NM),  XO(NM), 
YO(NM),  ZO(NM), 

XR(NM),  YR{NM), 

ZR(NM),  ALPHA(NM), 
BETA(NM),  GAM{NM), 
SCALE(NM).  DELTA(NM) 


(3)  Boundary-Related  Values 

Record  Number  Format 
1  NBOUN 


The  following  record  is  written  for  each  boundary  NB  from  1  to  NBOUN. 

2  (BNAME(IC,NB),IC=1,ICHAR), 
(IPARTMB(IC.NB),IC=  l.ICHAR), 

BTYPE(NB),  NHOLEIN(NB), 
((MHOLEIN(IC.NB.NH),IC  =  1  ,ICHAR), 

NH  =  l,NHOLEIN(NB)), 

XMAXH(NB),  YMAXH(NB),  ZMAXH(NB). 
XMINH(NB),  YMINH(NB),  ZMINH(NB), 
CLOSED(NB) 

3  NSURF 


The  following  record  is  written  for  each  surface  NS  from  1  to  NSURF 

4  (JRANGE{I,NS),  KRANGE(I,NS), 
LRANGE(I,NS),  1  =  1,2), 

NVOUT(NS), 

(IPARTBS(IC,NS),IC  =  l.ICHAR), 
(ISMESH(1C.NS),IC  =  1  ,ICHAR) 

5  NBOX 


The  following  record  is  written  for  each  box  NX  from  1  to  NBOX. 

6  {XRANGE(I,NX),  YRANGE(I,NX), 

ZRANGE{I,NX),  1  =  1,2), 
(IPARTBX(IC,NX),IC=  l.ICHAR) 
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(4)  Region-Related  Values 


Record  Number  Format 

1  NREGN 

The  following  record  is  written  for  each  region  NR  from  1  to  NREGN. 

2  (RNAME(IC,NR),IC  =  I  ,ICHAR), 
(1PARTMR(IC,NR),IC=  l.ICHAR), 

RTYPE(NR) 

3  NVOL 

The  following  record  is  written  for  each  volume  NV  from  1  to  NVOL. 

4  (JRANGEV(I,NV),  KRANGEV(I,NV), 

LRANGEV(I.NV),  1  =  1 ,2), 

(IPARTR(IC,N  V),IC  =  1  .ICHAR) 


(5)  Points  Requiring  Interpolation 

Record  Number  Format 

Records  1  and  2  are  written  for  each  mesh  NM  from  1  to  NMESH. 

1  NULIST(NM) 

The  following  record  is  written  only  if  NULIST(NM)  is  greater  than  0. 

2  (INLIST(I),I  =  I1,I2) 
where  II  =  1  and  12  =  NULIST(NM). 

(6)  Interpolation-List  Pointer 

Record  Number  Format 

The  following  record  is  written  for  each  recipient  mesh  NMR  (outer  loop) 
from  1  to  NMESH  and  each  donor  mesh  NMD  (inner  loop)  from  1  to 
•  NMESH. 

1  NU1NTR(NMR,NMD) 
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(7)  Phantom  Meshes 

Record  Number  Format 

The  following  records  are  written  for  each  mesh  NM  from  1  to  NMESH,  if 
PHANTOM(NM)  is  .TRUE. 

1  NM,  MJMAX(NM),  MKMAX(NM),  MLMAX(NM), 
MNAME 

2  (X(I),I  =  11,12), 

(Y(I),I  =  11,12), 

(Z(I),I  =  11,12) 

where  II  =  1  and  12  =  MJMAX(NM)*MKMAX(NM)*MLMAX(NM). 

(8)  Hole  Creation  Surfaces 

Record  Number  Format 

The  following  records  are  written  for  each  surface  NSURF  that  points  to 
a  hole  creation  boundary  (BTYPE=’HOLE’). 

1  NSPTS(NS) 

2  (XS(I),  YS{I),  ZS(I). 

VNX(I),  VNY(I),  VNZ(I),  1  =  11,12) 

where  II  =  1  and  12  =  NSPTS(NS). 

IBLANK  Array  for  Plotting  File 
This  file  is  written  to  file  IBPLOT. 

Record  Number  Format 

This  record  is  written  for  each  mesh  NM  (except  phantom  meshes)  from 
1  to  NMESH, 

1  (IBLANK(I),I  =  1,MJKL) 

where  MJKL  =  MJMAX(NM)*MKMAX(NM)'»MLMAX(NM). 
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Orphan  List  File 

This  file  is  written  to  rile  ORPHAN. 

Record  Number  Format 

This  record  is  written  for  each  orphan  point. 
1  MESH,  J,  K,  L 
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APPENDIX  G 

EXAMPLE  OF  COMPLEX  CONnCURATION 

This  example  is  taken  from  actual  computations  performed  on  the  Cray  X-MP/12  at 
AEDC.  PEGSUS  4.0  and  the  XMERA  flow  solver  (Ref.  5)  have  been  successfully  executed 
for  each  part  of  this  example.  The  meshes  used  were  generated  with  AEDC  codes  and  sized 
for  the  available  memory  of  the  AEDC  Cray.  An  improved  understanding  of  PEGSUS  4.0 
should  be  accomplished  by  careful  study  of  this  example,  just  as  was  the  case  in  Sec.  4.0 
of  the  main  text. 

The  example  is  presented  in  three  steps  to  demonstrate  the  ability  of  PEGSUS  4.0  to  handle 
complex  aerodynamic  configurations  as  well  as  to  illustrate  many  of  the  advanced  features 
of  PEGSUS  4.0.  The  three  steps  arc  shown  by  the  three  configurations  in  Fig.  G-1 .  Figure 
G- la  depicts  a  rectangular  cavity  mounted  below  a  flat  plate.  The  flow  is  assumed  to  be  above 
the  plate  and  within  the  cavity  only.  Figure  G-lb  shows  the  same  cavity  conHguration  with 
the  addition  of  a  store  mounted  on  a  bent  sting.  The  third  configuration  (Fig.  G-lc)  is  also 
the  same  cavity  and  store  but  with  the  sting  removed.  A  description  of  each  mesh  is  given 
as  well  as  the  user  input  required  to  execute  PEGSUS  4,0.  To  demonstrate  the  addition  and 
subtraction  capabilities  of  PEGSUS  4,0,  the  configuration  with  the  store/sting  (Fig.  G-lb) 
will  be  obtained  by  the  addition  of  the  appropriate  meshes  to  the  conflguration  of  the  empty 
cavity  (Fig.  G-1  a).  Similarly,  the  store  in  a  cavity  without  the  sting  (Fig.  G-lc)  will  be  obtained 
by  subtraction  and  addition  to  the  store/sting  in  the  cavity  (Fig.  G-lb).  These  cases  will  also 
demonstrate  the  use  of  the  double  fringe  for  all  interpolation  boundaries. 

G.l  Empty  Cavity 

To  begin,  the  domain  of  the  configuration  will  be  split  by  a  plane  of  symmetry,  requiring 
only  half  of  the  domain  to  be  gridded.  The  empty  cavity  configuration  is  divided  into  five 
computational  meshes,  two  meshes  defining  the  exterior  region  above  the  flat  plate,  and  three 
meshes  defining  the  cavity  and  an  overlap  region  above  the  cavity  mouth.  The  names  of  the 
five  meshes  are:  ‘EXTERIOR  A’,  ‘EXTERIOR  B’,  ‘CAVITY  A’.  ‘CAVITY  B’,  and  ‘CAVITY 
C’.  Shown  in  Figs.  G-2  through  G-6  are  (a)  the  basic  perspective,  (b)  a  side  view,  (c)  a  top 
view,  and  (d)  a  front  view  of  each  mesh.  The  cavity  meshes  are  shown  to  have  regions  that 
are  omitted  from  the  domain  by  using  interior  blanking.  This  technique  will  be  discussed 
later  in  conjunction  with  Fig.  0-9.  Table  G-1  gives  the  dimensions  of  each  mesh. 

The  user  input  required  to  connect  these  meshes  (see  Sec.  3.2)  is  given  in  Fig.  G-7.  The 
$GLOBAL  user  inputs  (see  Fig.  G-7a)  for  this  example  specify  that  a  double  fringe 
(FRINGE  =  2)  will  be  placed  around  all  holes.  The  QUALITY  user  input  specifies  that  the 
first  attempt  at  interpolation  of  a  point  will  be  done  with  a  QUALITY  of  1.0  and  the  second 
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attempt  will  be  done  with  a  QUALITY  of  0.5  (see  Sec.  2.2.3  for  the  discussion  on  QUALITY). 
Also  set  is  the  value  EPS  equal  to  O.OOOS.  The  user  inputs  XINCLUDE,  YINCLUDE,  and 
ZINCLUDE  are  set  to  reduce  the  number  of  orphan  points  that  can  occur  when  a  hole 
boundary  is  near  the  extreme  boundary  of  the  domain  (see  Sec.  2.3.2).  Also  shown  in  Fig. 
G-7a  are  the  user  inputs  required  for  each  mesh.  For  this  example,  only  the  LINK  values 
are  required. 

The  remainder  of  the  user  input  specifies  the  interpolation  boundary  points.  First,  the 
‘EXTERIOR  A’  mesh  will  interpolate  flow-field  values  from  ‘EXTERIOR  B’  as  well  as  from 
the  cavity  meshes.  To  obtain  interpolated  values  from  ‘EXTERIOR  B’,  an  outer  boundary 
is  specified  for  constant  J  planes  66  and  67  (see  Fig.  G-7b).  In  order  to  avoid  interpolation 
from  the  middle  of  an  expected  high  gradient  region  (the  region  about  the  L  - 1  plane  of 
the  ‘EXTERIOR  A'  mesh),  a  hole  creation  region  is  specified  within  the  ‘EXTERIOR  A’ 
mesh.  The  user  inputs  of  $REG10N  and  $  VOLUME  for  this  hole  creation  region  (see  Sec. 
2. 1 .2.2)  are  shown  in  Fig.  G-7b.  The  resulting  hole  that  is  made  in  the  ‘EXTERIOR  A’  mesh 
is  depicted  in  Fig.  G-8.  The  interpolation  boundary  specification  for  the  ‘EXTERIOR  B* 
mesh  is  similar  to  that  for  the  ‘EXTERIOR  A’  mesh,  with  an  outer  boundary  and  hole  creation 
region  being  specified  (see  Fig.  G-7c),  The  interpolation  boundary  points  within  the  cavity 
meshes  are  all  defined  with  outer  boundary  specifications  (see  Figs.  G-7d  through  G-7f).  For 
the  ‘CAVITY’ A’  and  ‘CAVITY  C’  meshes  five  surfaces  each  are  required.  ‘CAVITY  B’ 
requires  six  surfaces  to  define  its  outer  boundary.  Each  of  these  outer  boundary  surfaces 
is  composed  of  two  constant  planes  in  order  to  provide  a  double  interpolated  boundary  (double 
fringe). 

Another  PEGSUS  4.0  feature  used  for  this  example  is  interior  blanking  (see  Sec.  2.3. 1). 
Interior  blanking  was  used  to  keep  the  corner  of  the  cavity  within  a  single  mesh  that  is  not 
stretched  or  skewed  (i.e.,  this  region  is  not  split  with  one  mesh  defining  the  top  wall  and 
another  mesh  defining  the  side  wall).  Depicted  in  Fig.  G-9  is  a  portion  of  the  front  corner 
of  the  ‘CAVITY  A’  mesh  that  has  been  specified  by  using  the  SREGION  and  S VOLUME 
user  inputs  with  TYPE  =  ’INTR’  (see  Fig.  G-7d).  This  specification  sets  IBLANK=0  within 
the  interior  region.  Solid  wall  boundary  conditions  are  set  on  the  surfaces  of  this  region  within 
the  flow  solver.  Comer  regions  occur  for  the  front,  back,  and  side  walls  of  the  cavity.  This 
requires  the  ‘CAVITY  A’  and  ‘CAVITY  C’  meshes  to  have  two  SVOLUME  inputs  each  and 
‘CAVITY  B’  to  have  one  $VOLUME  input  (see  Figs.  G-7d  through  G-70. 

A  sampling  of  the  diagnostic  maps  for  the  empty  cavity  PEGSUS  execution  is  shown 
in  Fig.  G-10.  Shown  in  Fig.  G-lOa  is  the  legend  for  the  maps  giving  a  key  to  the  letters  and 
symbols  used.  The  maps  are  a  representation  of  the  computational  domain  for  each  L  plane 
of  each  mesh.  Only  selected  L  plane  maps  are  shown  here.  In  Figs.  G-lOb  and  G-lOc  are 
shown  L  =  10  maps  for  the  ‘EXTERIOR  A’  and  ‘EXTERIOR  B’  meshes,  respectively,  with 
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the  hole  and  its  double  fringe.  This  fringe  interpolates  from  the  ‘CAVITY  A’,  ‘CAVITY 
B’,  and  ‘CAVITY  C’  meshes.  The  outer  boundary  specified  for  the  ‘EXTERIOR  A’  mesh 
at  J  lines  66  and  67  interpolates  from  ‘EXTERIOR  B’  and  ‘CAVITY  B’.  For  the  ‘EXTERIOR 
B’  mesh,  the  outer  boundary  specified  for  J  lines  1  and  2  interpolates  from  ‘EXTERIOR 
A*  and  ‘CAVITY  B’.  In  Figs.  G-lOd  through  G-lOf,  the  L=24  diagnostic  maps  for  the 
‘CAVITY  A’,  ‘CAVITY  B’.  and  ‘CAVITY  C’  meshes,  respectively,  are  shown.  The  interior 
blanked  region  (‘#’  signs  with  no  fringe  points)  and  the  double  interpolation  boundaries  are 
clearly  evident. 

G.2  Store/Sting  in  a  Cavity 

Using  the  PEGSUS  results  from  the  empty  cavity  case,  the  store/sting  can  be  added.  The 
store  and  sting  are  each  represented  with  separate  meshes  which  are  shown  in  Figs.  G-1 1 
(‘STORE’)  and  G-12  (‘STING*),  respectively.  Another  mesh  (‘INTERFACE’)  was  also 
required  to  obtain  satisfactory  overlap  between  the  sting  and  exterior  meshes.  Without  the 
‘INTERFACE’  mesh,  the  hole  created  by  the  ‘STING’  mesh  in  the  ‘EXTERIOR  B’  mesh 
did  not  have  sufficient  overlap  between  the  hole  and  outer  boundaries  and  thus  created 
numerous  orphan  points.  The  ‘INTERFACE’  mesh  is  shown  in  Fig.  G-13.  The  dimensions 
for  all  three  meshes  are  given  in  Table  G-1.  In  Fig.  G-14,  the  user  input  for  specifying  this 
example  is  given.  All  three  new  meshes  are  added  by  specifying  MODE  =  ’ADD’  (see  Fig. 
G*14a).  To  adjust  the  LINK  lists  of  the  meshes  from  the  previous  PEGSUS  execution,  the 
ADDLINK  input  is  specified.  For  example,  the  ‘STORE’  mesh  has  ADDLINK  names  of 
‘EXTERIOR  A’,  ‘EXTERIOR  B’,  ‘CAVITY  A’,  ‘CAVITY  B’,  and  ‘CAVITY  C’.  From 
this  input,  the  mesh  name  ‘STORE’-  will  be  added  to  the  top  of  the  LINK  list  of  the  meshes 
specified  by  ADDLINK.  Since  no  hole  creation  boundaries  existed  in  the  previous  execution, 
the  ADDHOLE  input  was  not  required. 

The  interpolation  boundary  definitions  for  each  added  mesh  are  performed  in  the  same 
manner  as  for  an  initial  execution.  The  ‘STORE’,  ‘STING’,  and  ‘INTERFACE’  meshes  each 
make  holes  in  the  meshes  from  the  initial  execution  (example  G-1)  as  well  as  having  outer 
boundaries.  The  user  input  for  the  hole  creation  and  outer  boundary  speciHcations  is  given 
in  Figs.  G-14b  through  G-14d. 

After  creating  the  mesh  for  the  sting,  it  was  determined  that  part  of  the  ‘STING’  mesh 
was  inside  the  solid  back  wall  of  the  cavity.  Instead  of  creating  a  new  mesh,  a  hole  creation 
boundary  was  added  to  create  a  hole  in  the  ‘STING’  mesh.  Shown  in  Fig.  G-14e  is  the  user 
input  specifying  this  hole.  Surfaces  from  the  ‘CAVITY  C’  mesh  were  used  to  create  the  hole. 
This  hole  creation  boundary  was  not  closed  (see  Sec.  2. 1.2.1)  since  enough  of  the  surfaces 
was  available  to  create  the  proper  hole.  The  exterior,  cavity,  and  interface  meshes,  with  the 
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hole  created  by  the  ‘STORE’  and  ‘STING’  meshes,  are  shown  in  Fig.  G-15  for  the  plane 
of  symmetry. 

At  this  point,  a  detail  of  the  use  of  interconnecting  meshes  with  PEGSUS  needs  to  be 
discussed.  In  the  current  example  several  meshes  (e.g.,  ‘EXTERIOR  A’  and  ‘EXTERiOR 
B’  meshes)  are  overlapped  at  their  boundaries.  When  another  mesh  (e.g.,  the  ‘STORE’  mesh) 
crosses  this  boundary,  the  overlap  region  of  the  earlier  two  meshes  must  be  sufficient  to  avoid 
difficulties.  For  this  example,  which  uses  double  interpolation  boundaries,  the  ‘EXTERIOR 
A’  and  ‘EXTERIOR  B’  meshes  must  have  a  five-point  overlap,  i.e.,  five  points  (or  planes) 
that  have  identical  values.  If  a  four-point  overlap  were  to  be  used  (the  minimum  required 
for  a  double  fringe),  difficulties  can  arise  as  shown  in  Fig.  G-16.  First,  in  Fig.  G-16a  a  portion 
of  a  four-point  overlap  for  the  ‘EXTERIOR  A’  and  ‘EXTERIOR  B’  meshes  is  shown  as 
well  as  a  portion  of  the  outer  boundary  of  the  ‘STORE’  mesh.  In  Fig,  G-16b,  the  double 
interpolation  boimdary  points  for  the  ‘EXTERIOR  A’  mesh  are  shown  along  with  interpolation 
points  that  have  valid  or  invalid  interpolation  stencils  (as  far  as  the  ‘EXTERIOR  A’  mesh 
is  concerned)  for  the  ‘STORE’  mesh.  Points  PI  through  P3  are  interpolation  points  with 
valid  stencils  while  the  stencils  for  P4  through  P9  are  invalid.  Shown  in  Fig.  G-16c  are  the 
double  interpolation  boundary  of ‘EXTERIOR  B’,  with  valid  stencils  for  points  P6  through 
P9  and  invalid  stencils  for  points  PI  through  P5  of  the  ‘STORE’  mesh  outer  boundary.  For 
a  four-point  overlap,  then,  points  P4  and  P5  cannot  find  valid  interpolation  stencils  in  either 
niesh.  If  a  five-point  overlap  is  used  as  depicted  in  Fig.  G-17,  points  P4  and  PS  have  valid 
stencils  in  the  ‘EXTERIOR  B’  mesh  (see  Fig.  G-17c).  For  a  five-point  overlap,  then,  all  points 
are  found  to  have  valid  stencils  in  one  mesh  or  the  other  (see  Fig  G-17b  and  G-i7c).  Of  course, 
QUALITY  could  be  used  to  reduce  the  restrictions  for  valid  interpolation,  but  is  not  advised 
when  the  simple  solution  of  increasing  the  overlap  is  available. 

Figure  G-18  presents  a  sampling  of  diagnostic  maps  for  this  case.  In  Fig.  G-18a,  a  key 
for  the  letters  and  symbols  is  given.  In  Fig.  G-18b,  the  L  =  26  plane  is  shown  for  the 
‘EXTERIOR  A’  mesh.  Depicted  in  this  figure  are  the  outer  boundary  interpolation  points 
for  J  lines  66  and  67  and  the  interpolation  stencils  for  the  ‘EXTERIOR  B’  and  ‘INTERFACE’ 
meshes.  The  hole  made  by  the  ‘STING’  mesh  in  the  ‘EXTERIOR  B’  mesh  is  shown  in  Fig. 
G-18c.  Also  shown  are  the  outer  boundary  connections  to  the  ‘EXTERIOR  A’  mesh  and 
the  interpolation  stencils  for  the  ‘EXTERIOR  A’  and  ‘INTERFACE’  meshes.  The  hole  niade 
by  the  ‘INTERFACE’  mesh  in  the  ‘EXTERIOR  B’  mesh  is  depicted  in  Fig.  G-18d  for  the 
L  =  30  plane.  In  Figs.  G-I8e  through  G-18h,  the  hole  made  by  the  ‘STORE’  and  ‘STING’ 
meshes  in  the  ‘CAVITY  A’,  ‘CAVITY  B’,  and  ‘CAVITY  C’  meshes  is  shown.  The  L=29 
outer  plane  of  the  ‘STORE’  mesh  is  depicted  in  Fig.  G-18i.  As  can  be  seen  in  this  figure, 
the  ‘STORE’  mesh  interpolates  flow-field  values  from  the  ‘EXTERIOR  A’,  ‘EXTERIOR 
B’,  ‘CAVITY  A’,  ‘CAVITY  B’,  and  ‘CAVITY  C’  meshes.  In  Fig.  G-I8j,  the  L=29  outer 
plane  for  the  ‘STING’  mesh  is  shown.  The  main  item  to  note  from  this  figure  is  the  hole 


147 


i 


AEDC-TR-91-8  I 

that  appears  at  J  =  18  and  J  =  19.  This  hole  was  made  by  the  hole  creation  boundary  ‘AFT 
CAVITY  HOLE  BOUNDARY*.  Finally,  in  Figs.  G-18k  through  G-18ni,  three  L  planes  of 
the  ‘INTERFACE’  mesh  are  given  to  show  the  hole  created  by  the  ‘STING’  mesh  and  the 
interactions  of  the  ‘INTERFACE’  mesh  with  the  ‘EXTERIOR  A’  and  ‘EXTERIOR  B’  meshes. 

G.3  Store  Alone  in  a  Cavity 

The  final  example  demonstrates,  the  capability  of  PEGSUS  4.0  to  perform  subtraction 
and  addition  during  the  same  execution.  The  ‘STING’  and  ‘INTERFACE’  meshes  will  be 
subtracted  while  the  ‘STORE  CAP’  mesh  will  be  added.  A  description  of  the  ‘STORE  CAP’ 
mesh  is  shown  in  Fig.  G-19.  Table  C^-l  contains  the  dimensions  of  this  mesh.  The  user  input 
for  this  case  is  given  in  Fig.  G-20.  The  mesh  definitions  given  in  Fig.  G-20a  show  subtraction 
of  the  ‘STING’  and  ‘INTERFACE’  meshes  by  using  the  input  MODE  =  ’SUB*.  All  boundaries 
associated  with  these  two  meshes  are  subtracted.  The  ‘STORE  CAP’  mesh  is  added  to  the 
configuration  with  the  associated  LINK  and  ADDLINK  inputs,  as  shown  in  Fig.  G-20a.  In 
Fig.  G-20b,  the  specification  of  the  outer  boundary  for  the  ‘STORE  CAP’  mesh  and  interior 
blanked  region  is  given.  The  interior  blanked  region  for  the  ‘STORE  CAP’  is  used  to  obtain 
the  required  overlap  for  the  double  interpolation  boundary  used  by  the  ‘STORE  CAP’  and 
‘STORE’  meshes.  An  additional  outer  boundary  is  added  to  the  ‘STORE’  mesh  in  Fig.  G-20c 
to  correctly  specify  the  portion  of  the  L  =  1  surface  of  the  ‘STORE’  mesh  that  was  previously 
part  of  the  sting  surface.  Finally,  the  ‘AFT  CAVITY  HOLE  BOUNDARY’  of  the  ‘CAVITY 
C’  mesh  is  subtracted  in  Fig.  G-20c  since  the  ‘STING*  mesh  is  no  longer  present.  The  hole 
that  now  exists  in  the  plane  of  symmetry  for  the  exterior  and  cavity  meshes  is  shown  in  Fig. 
G.21. 


Figure  G-22  is  a  sampling  of  diagnostic  maps  for  this  case.  Figure  G-22a  gives  a  key  to 
the  letters  and  symbols  used.  Most  of  the  diagnostic  maps  have  no  changes  from  those  seen 
for  example  G-2,  with  the  exception  of  the  subtraction  of  the  ‘STING’  and  ‘INTERFACE’ 
meshes.  Figure  G-22b  shows  the  hole  made  by  the  ‘STORE’  mesh  in  the  L  =  56  plane  of  the 
‘CAVITY  C’  mesh.  In  Fig.  G-22c,  the  L  =  27  plane  for  the  ‘STORE  CAP’  mesh  is  shown 
to  highlight  the  interior  blanked  region. 
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a.  Cavity 

Figure  G-1.  Graphical  depiction  of  complex  configuration. 
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Flow 


b.  Cavity/store/sling 
Figure  G-1.  Continued. 


Flow 


c.  Cavily/slore 
Figure  G-1.  Concluded. 
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Flow 


depiction 
RIOR  A’  mesb 
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b.  Side  view 
Figure  G-2.  Continued. 


c.  Top  view 
Figure  G-2.  Continued. 


Concluded. 
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a.  Graphical  depiction 
Figure  G*3.  ‘EXTERIOR  B’  mesh. 
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b.  Side  view 
Figure  G-3.  Continued. 
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c.  Top  view 
Figure  G-3.  Continued. 


d.  Front  view 
Figure  G-3.  Concluded. 
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Flow 


a.  Graphical  depiction 
Figure  G-4.  ‘CAVITY  A*  mesh. 
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J 


with  Interior  Blanking 


b.  Side  view 
Figure  G-4.  Continued. 


AEDC-TR-91-8 


c.  Top  view 
Figure  G-4.  Continued. 


K 


L 


Gridded  Region 
with  Interior  Blanking 


d.  Front  view 
Figure  G-4.  Concluded. 
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a.  Graphical  dcpiclioa 
Figure  G-5.  ‘CAVITY  B’  mesh. 
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AEDC-TR-91-8 


J 


b.  Side  view 
Figure  G-5.  Continued. 


c.  Top  view 
Figure  G-5.  Continued. 
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AEDC-TR-91-8 


K 


L 

Gridded  Region 
with  Interior  Blanking 


d.  Front  view 
Figure  G-5.  Concluded. 
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Flow 


J  =  26 


B.  Graphical  depiction 
Figure  G-6,  ‘CAVITY  C’  mesh. 


c.  Top  view 
Figure  G-6.  Continued. 


d.  Front  view 
Figure  G-6.  Concluded. 
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(;»*****•***  *****•*******'************************•*****  ••*•*** 

C 

C  PEGSUS  4.0  ZHPaT 

C 

C  EXAMPLE  6-1  EMPTY  CAVITY 

C  5  MESH  CASE 

C 

C************************************************************ 


c********** ************************************************** 

c 

C  GLOBAL  DEFZNZTZOHS 
C 

C»1I  ****************** **************************************** 


$GLOBAL 

FRINGE  -  2, 

QUALITY  -  1.0, 0.5, -0.5, 

EPS  -  0.0005, 

XIHCLDDE  -  -30.0,46.999, 

YINCLUDE  =  -30.0,0.0, 

ZIHCLUDE  -  -10.0,15.999, 

$EMD 

C****** ****************************************************** 

C 

C  MESH  DEFINITIONS 
C 

C******  *******  *****•****(******■*********************★********** 


$HESH  NAME  - 
LINK  -> 


'EXTERIOR  A', 
'CAVITY  C', 
'CAVITY  B', 
'CAVITY  A', 

' EXTERIOR  B' , 


$MESH  NAME  = 
LINK  = 


'EXTERIOR  B', 
'CAVITY  A', 
'CAVITY  B', 
'CAVITY  C', 
'EXTERIOR  A', 


$ME5H  NAME 
LINK 


' CAVITY  A' , 
'CAVITY  B', 
'EXTERIOR  A', 
'EXTERIOR  B', 


$MESH  NAME  -  'CAVITY  B' , 
LINK  -  'CAVITY  A', 
'CAVITY  C', 

' EXTERIOR  A'  , 
'EXTERIOR  B', 


$HESH  NAME 
LINK 


'CAVITY  C', 
'CAVITY  B', 
'EXTERIOR  B' , 
'EXTERIOR  A', 


$BND 


$EMD 


$EHD 


$END 


$END 


a.  Page  1 

Figure  G-7.  User  Input  for  empty  cavity. 
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AEDC-TR-91-a 


^**** **«**«***«*«*«**««*«•***«***** ************************** 
C 

C  INTERPOLATION  BOCNDARY  DEFINITIONS 
C 

C************************************************************ 


C++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
C  EXTERIOR  A 


C - 

c  OUTER  BOUNDARY  DEFINITION 

C - 


$BOUNDARY  NAME  =  'EXTERIOR  A  OUTER  BOUNDARY', 
ISPARTOF  -  'EXTERIOR  A', 

SEND 


$ SURFACE  ISPARTOF  =  'EXTERIOR  A  OUTER  BOUNDARY', 
JRANGE  -  66,67, 

KRANGE  2,25, 

LRANGE  °  2,35, 


SEND 

- - 

C  HOLE  REGION  DEFINITION 

C - 


SREGION  NAME  -  'EXTERIOR  A  HOLE  REGION', 
TYPE  =  'HOLE', 

ISPARTOF  -  'EXTERIOR  A', 

Send 


$VOLUME  ISPARTOF 
ORANGE 
ECRANGE 
LRANGE 


'EXTERIOR  A  HOLE  REGION', 
27,67, 

12,26, 

1,20, 

Send 


b.  Page  2 

Figure  G-7.  Continued. 
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o  u  o  u  o  6 


AEDC-TR-01-8 


^BOUNDARY  NAME  *  'EXTERIOR  B  OUTER  BOUNDARY' , 
ISPARTOF  =  'EXTERIOR  B', 

$END 

SSURFACE  ISPARTOF  =•  'EXTERIOR  B  OUTER  BOUNDARY', 
ORANGE  =  1,2, 

KRANGE  -  2,25, 

LRANGE  »  2,35, 

$END 

HOLE  REGION  DEFINITION 


$REG10N  NAME  =  'EXTERIOR  B  HOLE  REGICM' , 
TYPE  -  'HOLE', 

ISPARTOF  =  'EXTERIOR  B', 

^  Send 

SVOLUME  ISPARTOF  -  'EXTERIOR  B  HOLE  REGICXI' , 
orange  =  1,34, 

KRANGE  -  12,26, 

LRANGE  =  1,20, 

Send 

c.  Page  3 

Figure  G-7.  Continued. 
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u  u  u  o  o  o 


AEDC-TR-91-8 


C  CAVITY  A 

C+++ +++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++ 


OUTER  BOUNDARY  DEFINITION 


^BOUNDARY 


ISPARTOF  ■»  'CAVITY  A' 


t 

$EMD 


$SORFACE  ISPARTOF  -  'CAVITY  A  OUTER  BOUNDARY', 
JRANGE  °  1,2, 

KRANGE  °  1,30, 

LRANGE  -  37,85, 

$END 


$SURFACE  ISPARTOF  °  'CAVITY  A  OUTER  BOUNDARY', 
JRANGE  »  29,30, 

KRANGE  °  5,30, 

LRANGE  °  2,56, 

SEND 


SSURFACE  ISPARTOF  =  'CAVITY  A  OUTER  BOUNDARY', 
JRANGE  -  29,30, 

KRANGE  -  1,30, 

LRANGE  •>  57,85, 

SEND 

SSURFACE  ISPARTOF  =  'CAVITY  A  OUTER  BOUNDARY', 
JRANGE  o  3,28, 

KRANGE  »  1,2, 

LRANGE  o  57,85, 

SEND 

$SURFACE  ISPARTOF  -  'CAVITY  A  OUTER  BOUNDARY', 
JRANGE  <=  3,28, 

KRANGE  •>  3,30, 

LRANGE  °  84, 85, 

SEND 


INTERIOR  REGION  DEFINITION 


SREGION  NAME  =  'CAVITY  A  INTR  REGION', 
TYPE  -  'INTR', 


ISPARTOF 

“  'CAVITY 

A', 

SEND 

1 

O 

> 

ISPARTOF 

S) 

'CAVITY 

A 

INTR  REGION' , 

JRANGE 

ta 

1,5, 

KRANGE 

o 

1,31, 

LRANGE 

" 

1/56, 

SEND 

< 

1 

ISPARTOF 

'CAVITY 

A 

INTR  REGION', 

JRANGE 

83 

1,30, 

KRANGE 

- 

1,5, 

LRANGE 

SS 

1,56, 

SEND 

d.  Page  4 

Figure  G-7.  Continued. 
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o  o  o  C* 


AEDC-TR-91-8 


C  CAVITY  B 

C++ ++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++ 


OOTER  BOUNDARY  DEFINITION 


^BOUNDARY  NAME  -  'CAVITY  B  OUTER  BOUNDARY', 
ISPARTOF  =  'CAVITY  B' , 

$EHD 


^SURFACE  ISPARTOF 
JRANGE 
KRAMSE  > 
LRANGE  = 


$ SURFACE  ISPARTOF 
JRANGE  - 
KRANGE  - 
LRANGB  > 


^SURFACE  ISPARTOF 
JRANGE  > 
KRANGE  - 
LRANGE  • 


$ SURFACE  ISPARTOF 
JRANGE  « 
KRANGE  c 
LRANGE 


$ SURFACE  ISPARTOF 
JRANGE  > 
KRANGE  - 
LRANGE  - 


$SURFACE  ISPARTOF 
JRANGE  c 
KRANGE  - 
LRANGE  - 


=  'CAVITY 
1,2, 
5,30, 
2,56, 


■  'CAVITY 
1,2, 
1,30, 
57,85, 


»  'CAVITY 

29.30, 

5.30, 
2,56, 


-  'CAVITY 

29.30, 

1.30, 
57,85, 


■  'CAVITY 
3,28, 
1,2, 
57,85, 


•  'CAVITY 
3,28, 
3,30, 
84,85, 


B  OUTER  BOUNDARY', 

$END 

B  OUTER  BOUNDARY', 

Send 

B  OUTER  BOUNDARY', 

SEND 

B  OUTER  BOUNDARY', 

SEND 

B  OUTER  BOUNDARY', 

SEND 

B  OUTER  BOUNDARY', 

SEND 


INTERIOR  REGION  DEFINITION 


SREGION  NAME  =  'CAVITY  B  INTR  REGION', 
TYPE  -  'INTR', 

ISPARTOF  =  'CAVITY  B', 


$ VOLUME  ISPARTOF 
JRANGE 
KRANGE 
LRANGE 


SEND 

«  'CAVITY  B  INTR  REGION' 

-  1,30, 

-  1,5, 

“  1,56, 


SEND 


, 


e.  Page  5 

Figure  G-7.  Continued. 
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o  n  n 


AEDC-TR-91-8 


C - 

C  OUTER  BOUNDARY  DEFINITION 

- - 

^BOUNDARY  NAME  -  'CAVITY  C  OUTER  BOUNDARY', 
ISPARTOF  •=  'CAVITY  C'V 

$END 

^SURFACE  ISPARTOF  -  'CAVITY  C  OUTER  BOUNDARY', 
ORANGE  -  1,2, 

KRANGE  =  5,30, 

LRANGE  =  2,56, 

$END 

5SURFACB  ISPARTOF  =  'CAVITY  C  OUTER  BOUNDARY', 
ORANGE  -  1,2, 

KRANGE  o  1,30, 

LRANGE  =  57,85, 

$END 

SSURFACE  ISPARTOF  -  ' CAVITY  C  OUTER  BOUNDARY' , 
ORANGE  1  29,30; 

KRANGE  -  1,30, 

LRANGE  >>  57,85, 

SEND 

SSURFACE  ISPARTOF  =  'CAVITY  C  OUTER  BOUNDARY', 
ORANGE  -  3,28, 

KRANGE  =  1,2, 

LRANGE  ==  57,85, 

SEND 

SSURFACE  ISPARTOF  -  'CAVITY  C  OUTER  BOUNDARY', 
ORANGE  -  3,28, 

KRANGE  -  3,30, 

LRANGE  -  84,85, 

SEND 


INTERIOR  REGION  DEFINITION 


SREGION  NAME  -  'CAVITY  C  INTR  REGION', 

TYPE  -  'INTR', 

ISPARTOF  =  'CAVITY  C' , 

SEND 

SVOLUME  ISPARTOF  •>  'CAVITY  C  INTR  REGION', 
ORANGE  -  26,30, 

KRANGE  -  1,31, 

LRANGE  °  1,56, 

SEND 

SVOLUME  ISPARTOF  -  'CAVITY  C  INTR  REGION', 
ORANGE  «  1,30, 

KRANGE  =  1,5, 

LRANGE  -  1,56, 

Send 

C  END  OF  EMPTY  CAVITY  CASE 

f.  Page  6 

Figure  G-7.  Concluded. 
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Figure  G-8.  Hole  creation  region  in  ’EXTERIOR  A’  mesh. 
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Figure  G>9.  Interior  blanking  of  comer. 
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OIAGHQSTIC  MAPS 

LEGEND 

BODNDARY  POINTS  UPDATED  BY  : 

A  -  EXTERIOR  A 
B  -  EXTERIOR  B 
C  -  CAVITY  A 
D  -  CAVITY  B 
E  -  CAVITY  C 


INTERPOLATION  STENCILS  UPDATING  POINTS  IN  : 
a  -  EXTERIOR  A 
b  -  EXTERIOR  B 
c  -  CAVITY  A 
d  -  CAVITY  B 
e  -  CAVITY  C 


.  -  FIELD  POINT 

f  -  HOLE  POINT 

7  -  ORPHANED  BOUNDARY  POINT 

*  -  INTERPOLATION  STENCIL  UPDATING  MORE  THAN  ONE 

INTERPOLATION  BOUNDARY  POINT 

a.  Page  1 

Figure  G-10.  Sample  diagnostic  maps  for  empty  cavity. 
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HhP  FOB.  EXTEBIOR  A 


PLANE  L  -  10 


J123456789 

123456789O123456789O123456789O123456789O123456789O123456789O1234567B9O123456789O12345€709O123<5G7B 

K 

2  i !!!”!!  i !'.!  1  !!!!!!!!!!!!!!*!!!,!!  i  i "  i  i  i  •  . . bb.Bs 

3  -bb-BB 

4  .  .  .  . . . bb.BB 

5  . bb.BB 

g  "liiiiii” . * . bb.BB 

7  !!!!!!!!!!!!!! . . .cccc.*. *.*.*.*.*.*.*.*•*** ***»*-*«*'*-*d**.BB 

B  \\[ . *cccc.*. *.*.*.*.•.*.*.*.****** •*■*•*•*■*<***•88 

10  i . * _ CCOCCCCCCOCCCCCCCCCCCCDDDDDDDDDDDDDDDDDDD 

11  !!!!!!!!!!!!!. .....■.*.. .ccoccccccccccccccccccccddddddddddddddddddd 

12  . *..CCi###«#ftl######H###l########i##i###ti### 

13  . 

14  . *..CC######ll###############»#i#i#il»#«#*i»#* 

15  . *..CCI#####l»###«###########t##fl##fH######it 

16  . *  .  .CC»##4#ll»»##m#f 

17  . *..CCt###«fflf#*lf<lffif##lfifl##«*t«f*«lffi*f 

18  . *  .  .OC#####t»########»####t#»#»l#»t##i###tf  #### 

19  . *-  .CC#####l»i##4#H#####f####»######»#»*»»l»** 

20  . *. 

21  . *.  .CC##4#######4l##################i#ll»###**» 

22  . *. 

23  . * .  .CC#HI#l####«#####f 

24  . *..CC4##########ltt#»#Ht###»########*l#HH»## 

25  . 

26  . ####l###i#tlt##f###i##»########*l###«»#t 

b.  Page  2 

Figure  G«10.  Contioued. 


HAP  FOR  EXTERIOR  B 


PLANE  L  °  10 


QO 


J1234567B9 

123«6789012345678901234Sfi7890123456789012345678901234S67990123«€7B9012345«7e9012345679901234567B 

K 

1  . 

2  AA.aa . . . . . 

3  AA.aa« . . . . . . . 

4  AA.aa . . . 

5  AA.aa . 

6  AA.aa . i . . 

7  AAd**dddd*'***MaaaaMaaMaaMaaaa*.*  . .  - 

8  AAd**dddd****aaaaa«*eeeeee6eeee8a«.* .* . . . . 

9  DD . . . *-.» . 

10  DDDODDDDDDDDEEEEBEEEEEEEEBEEEBEEEE. . . * . 

11  DDDDDDDDDDDDEEEEEEEEEBBEBEEEEEEEEEE . . * . . . 


12 

ttitItfftIHfIfifEE 

13 

ftfillfillittfft 

fIfffIffffIffffllEB 

14 

*#«4«l»«««4tif«t 

ftffl«ftlff««««t«EE 

15 

fflttlfflfflftf# 

16 

fl»#4«4tlff«44tftEE 

17 

««4ilf*flflff«lt 

tffl44lfll«#4ftl*EB 

18 

«444it#«««fl#«l« 

t*««<*tf»ff<44tt*EE 

19 

«444llf#4ftlf«lf 

tf««lf#fff««4f#f#BE 

20 

fItilffifttIHI# 

fffffttIfffifttflEE 

21 

«44lit«<«ttttf4* 

ffifftlfffltiflffEE 

22 

|«ftf»«4ftll*f<f 

#ffl4ltff#«44ilf#EE 

23 

44lttt«««ttH*#* 

«f«4ftt*f««<*ff>#EE 

24 

ffftIHIfItffIff 

«fl44llff«««ffi*«EE 

25 

4<ftl««fftf»flft 

«««««tftfl4fl4#«#EE 

26 

ffitififlltfl*## 

f<l44lffH4fllt#<.. 

c.  Pi^  3 

Figure  G'lO.  Continued. 
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MAP  FOR  CAVITY  A 


PLANE  L  •  24 


j12345^78® 

12345678901234567B9012345678901234S67e901234567e9012345678901234S€7e901234567890123456789012345«78 

2  *«««««**#««tt«t««ftlffi444tti# 

3  *«««flt««ftlf*fft«««««iff*tf** 

4  IfilffffiffifHtfllfllfllltff# 

5  «4l«tt#««tfl«lf««*f<l4f*tddt«« 


6  f«l«t . dd.DD 

7  4 lift . dd.DD 

8  4«4fl . dd.DD 

9  ilftf . <“-DD 

10  . . dd-DD 

11  44#tt . dd-DO 

12  ffttf . dd.DD 

►r  13  4ttl« . dd.DD 

N©  H  tlfll . dd.DD 

15  till* . dd.DD 

16  ##IH . dd.DD 

17  fif«4 . dd.DD  . 

18  ttt«f . dd.DD 

19  f»f«4 . dd.DD 

20  ff««f . dd.DD 

21  «»««! . dd.DD 

22  ##### . dd.DD 

23  I4««t . dd.DD 

24  l«44l . '. . dd.DD 

25  14441 . dd.DD 

26  44444 . dd.DD 

27  44ttl . dd.DD 

28  44444 . dd.DD 

29  44414 . dd.DD 

30  44414 . dd.DD 

31  44444 . 


d.  Page  4 

Figure  G-10.  Continued. 
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MAP  FOR  CAVITY  B 


PLANE  L  -  24 


J12345e7g9 

1234567e90123456789012345«7S901234S678901234567890123456789012345C7e90123456789012345C789012345678 

K 

1  ff#«ttlt*itff«««l#tf#il««tf»H 

2  ff##«ttl«4lf«««*«llff«f44tti«* 

3  fffffHffffftftftllfMItfftfl* 

4  t«»««4ii««#iffl«l««tf##«fittff 

5  tftccftffffttfffttItfflltMtfl 


6  CC.cc . . .ee. EE 

7  CC.cc . m.EE 

6  CC.cc . ee.EE 

9  CC.cc . ■•.EE 

10  CC.cc . e«.EE 

11  CC.cc . •■.££ 

12  CC.cc . ««.EE 

»  13  CC.cc . ......«a.EE 

O  14  CC.CC . ee.EE 

15  CC.cc . ..........ee.EE 

16  CC.cc . ee.EE 

17  CC.cc . ee.EE 

18  CC.cc...... . ...ee.EE 

19  CC.cc . ee.EE 

20  CC.cc . ee.EE 

21  CC.CC . ee.EE 

22  CC.cc . ......ee.EE 

23  CC.cc . ee.EE 

24  CC.cc . ee.EE 

25  CC.cc . ee.EE 

26  CC.cc . ee.EE 

27  CC.cc . ea.EE 

26  CC.cc . ee.EE 

29  CC.cc . ee.EE 

30  CC.cc . ee.EE 

31  . 


e.  Page  5 

Figijure  G<10.  Continued 
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HAP  FOR  CAVITY  C 


PLANE  L  -  24 


J 


K 

1  4iltlli»f*##4ft#f*«itl«*l*l«« 

2  44tt«4«tf««tftlfff«iit«N«l«t 

3  l4tt»l*tf««««tfff««ttff«MI«* 

4  ftl*<4«tf<4fl«fl««4ti«<t*i«tt 

5  4ttdd«ttf«tftttff4«lf*«f«i«*« 


6  DD.dd . fitti 

7  DD.dd . ft4l< 

8  DD.dd . tHf< 

9  DD.dd . Itl#< 

10  DD.dd . tlfli 

11  DD.dd . «###] 

12  DD.dd . <ltt< 

^  13  DD.dd . ###*( 

M  14  DD.dd . Itl#< 

15  DD.dd . .Itif' 

16  DD.dd . *#*# 

17  DD.dd . Mff' 

18  DD.dd . flit 

19  DD.dd . ttft 

20  DD.dd . «!## 

21  DD.dd . Ill# 

22  DD.dd . Iff* 

23  DD.dd . lit# 

24  DD.dd . t«f« 

25  DD.dd . flff 

26  DD.dd . fif# 

27  DD.dd . #t#i 

28  DD.dd . «#«# 

29  DD.dd . ftll 

30  DD.dd . itll 

31  . . . #### 


f.  Page  6 

Figure  G-10.  Concluded. 
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a.  Graphical  depiction 
Figure  G-11.  'STORE’  mesh. 


AEDC-TR-gi-8 


c.  Top  view 

Figure  G-11.  Continued. 


AEDC-TR-9ie 


d.  Front  view 
Figure  G-11.  Concluded. 
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AEDC-TR-91-8 


AEDC-TR-91-8 


d.  Front  view 
Figure  G>12.  Concluded. 
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Flow 


a.  Graphical  depiction 
Figure  G-13.  ‘INTERFACE’  mesh. 


b.  Side  view 
Figure  G*13.  Continued. 


c.  Top  view 

Figure  G-13.  Continued. 
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AEDC-TR-91-8 


d.  Front  view 
Figure  G-13.  Concluded. 
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AEDC>TR-91>8 


****★♦*♦*♦************#*★*★*★*★*♦******«******★********★♦★♦* 
PEGSUS  4.0  INPUT 

EXAMPLE  G-2  EMPTY  CAVITY  WITH  STORE /STING 
8  MESH  CASE 

STORE,  STIHG,  AND  INTERFACE  MESHES  ADDED  TO  5  MESH 
EMPTY  CAVITY 

**4kf«r*<lr'*««******«*«*«***1ir**«r*ir*********1k*********#****Al^l^**i^ 
*************  ********iHt*****iH*fHr*Hr***************»********** 

MESH  DEFINITIONS 

(;**********«*«*********************«***«*****«***•*«•******** 

$MESH  NAME  -  'STORE', 

MODE  '  ADD'  , 

LINK  “  'CAVITY  A', 

' STING' , 

' CAVITY  B' , 

'CAVITY  C', 

' EXTERIOR  A' , 

' EXTERIOR  B' , 

' INTERFACE' , 

ADDLINK  -  'EXTERIOR  A', 

'EXTERIOR  B', 

'CAVITY  A', 

'CAVITY  B', 

'CAVITY  C', 

XO  «  IS. 468,  YO  °  0.0,  ZO  =  0.0, 

$END 

SMESH  NAME  •>  '  STING'  , 

MODE  -  'ADD' , 

LINK  -  'CAVITY  C', 

'EXTERIOR  B', 

' INTERFACE' , 

' STORE' , 

'CAVITY  B', 

'EXTERIOR  A', 

ADDLINK  =  'EXTERIOR  A', 

'EXTERIOR  B', 

'CAVITY  B', 

'CAVITY  C', 

YINCLODE  =  *30.0,30.0, 

KIHCLODE  -  2,32, 

XO  -  IS. 468,  YO  -  0.0,  ZO  =  0,0, 

SEND 

$MESH  NAME  -  ' INTERFACE' , 

MODE  -  'ADD' , 

LINK  =  'EXTERIOR  B', 

' STING' , 

'CAVITY  C', 

'EXTERIOR  A', 

' STORE' , 

'CAVITY  B' , 

ADDLINK  -  'EXTERIOR  A', 

'EXTERIOR  B', 

'CAVITY  B', 

'CAVITY  C', 

SEND 

a.  Page  1 

Figure  G-14.  User  input  for  store/sting  in  a  cavity. 
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Q******<lr*<»it****'****«**************************************'*** 

C 

C  INTERPOLATION  BODHDARY  DEFINITIONS 
C 

Q*  *-************«*'**'**'******'**'***********•  *«•*****•**•*•«*•••• 


C++++++++ ++++++++++++++++++ ++++++++++++++++++++++++++++++++++ 
C  STORE 

C++++++++ ++++++++++++++++++++++++++++ ++++++++++++++++++++++++ 

c - 

C  HOLE  CREATION  BOUNDARY  DEFINITION 

C - 


SBOUNDARY  NAME  -  'STORE  HOLE  BOUNDARY', 
HHOLEIN  =  'CAVITY  A', 

'CAVITY  B', 

'CAVITY  C', 

'EXTERIOR  A', 
'EXTERIOR  B', 

ISPARTOF  -  'STORE', 

$EHD 

$SURFACE  ISPARTOF  »  'STORE  HOLE  BOUNDARY', 
JRANOE  -  1,82, 

KRAN6B  -  2,32, 

LRAN6E  o  5,5, 

NVOOT  -  '  +L'  , 

$END 

$SURFACE  ISPARTOF  -  'STORE  HOLE  BOUNDARY', 
ORANGE  -  82,82, 

KRANGE  -  2,32, 

LRAHGE  -  1,5, 

NVOUT  -  '+J', 

SEND 


C - 

C  OUTER  BOUNDARY  DEFINITION 

C - 


$BOUNDARY  NAME  »  'STORE  OUTER  BOUNDARY', 
ISPARTOF  «  'STORE', 

$END 

$ SURFACE  ISPARTOF  -  'STORE  OUTER  BOUNDARY', 
ORANGE  -  1,82, 

KRANGE  o  2,32, 

LRANGE  -  28,29, 

SEND 

$ SURFACE  ISPARTOF  =  ' STORE  OUTER  BOUNDARY' , 
ORANGE  -  81,82, 

KRANGE  «  2,32, 

LRANGE  -  1,28, 

SEND 

b.  Page  2 

Figure  G-14.  Continued. 
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AEDC-Tn-91-8 


C++++++++++++++++++++++++++++++++++++++++ ++++++++++++++ ++++++ 
C  STING 

C++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


HOLE  CREATION  BOUNDARY  DEE INI T1 OH 


^BOUNDARY  NAME  »  'STING  HOLE  BOUNDARY', 
MHOLBIN  ■>  'EXTERIOR  A', 
'EXTERIOR  B', 

'CAVITY  B', 

'CAVITY  C', 

' INTERFACE'  , 

ZSPARTOF  -  'STING',  ' 

SEND 

$SURFACE  ISPARTOF  -  'STING  HOLE  BOUNDARY', 
ORANGE  a  1,82, 

KRANGE  -  2,32, 

LRANGE  «  5,5, 

NVOUT  =  '+!', 

SEND 

SSURFACE  ISPARTOF  -  'STING  HOLE  BOUNDARY', 
ORANGE  -  1,1, 

KRANGE  °  2,32, 

LRANGE  B  1,5, 

NVOUT  ■  '-O', 

SEND 

SSURFACE  ISPARTOF  -  'STING  HOLE  BOUNDARY', 
ORANGE  B  82,82, 

KRANGE  •  2,32, 

LRANGE  •>  1,5, 

NVOUT  -  '+0', 

SEND 


OUTER  BOUNDARY  DEFINITION 


SBOUHDARY  NAME  »  'STING  OUTER  BOUNDARY', 
ISPARTOF  -  'STING', 

SEND 

SSURFACE  ISPARTOF  -  'STING  OUTER  BOUNDARY', 
ORANGE  -  1,81, 

KRANGE  B  2,32, 

LRANGE  B  28,29, 

SEND 

SSURFACE  ISPARTOF  -  'STING  OUTER  BOUNDARY', 
ORANGE  B  1,2, 

KRANGE  B  2,32, 

LRANGE  B  1,28, 

SEND 

c.  Page  3 

Figure  G-14.  Continued. 
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C+  +4+++++++++++ 4+ 4444 44444444444444444444 44 444444 444444444444 
C  INTERFACE 

c 44 44 44444444444444 44444444444444444444 44 44444444444444444444 


HOLE  CREATION  BOUNDARY  DEFINITION 


$BOUNDARY  NAME  <>  'INTERFACE  HOLE  BOUNDARY', 
HHOLBIN  »  ' EXTERIOR  A'  , 

'EXTERIOR  B', 

ISPARTOF  -  'INTERFACE', 

$EHD 

$ SURFACE  ISPARTOF  »  'INTERFACE  HOLE  BOUNDARY', 
JRANGE  -  10,10, 

KRANGE  «  7,21, 

LRANGE  »  6,35, 

NVOUT  “  ' -J' , 

$EHD 

$ SURFACE  ISPARTOF  »  'INTERFACE  BOLE  BOUNDARY', 
JRANGE  -  10,77, 

KRANGE  -  7,7, 

LRANGE  -  6,35, 

NVOUT  -  '-K', 

$END 

SSURFACE  ISPARTOF  -  'INTERFACE  HOLE  BOUNDARY', 
JRANGE  *  10,77, 

KRANGE  =  7,21, 

LRANGE  -  6,6, 

NVOUT  “  '-L', 

Send 


OUTER  BOUNDARY  DEFINITION 


SBOUNDARY  NAME  -  'INTERFACE  OUTER  BOUNDARY', 
ISPARTOF  °  'INTERFACE', 

SEND 

SSURFACE  ISPARTOF  =  'INTERFACE  OUTER  BOUNDARY', 
JRANGE  -  1,2, 

KRANGE  -  1,20, 

LRANGE  -  1,35, 

SEND 

SSURFACE  ISPARTOF  -  'INTERFACE  OUTER  BOUNDARY', 
JRANGE  -  1,76, 

KRANGE  -  1,20, 

LRANGE  <•  1,2, 

SEND 

SSURFACE  ISPARTOF  -  'INTERFACE  OUTER  BOUNDARY', 
JRANGE  -  1,76, 

KRANGE  =  1,2, 

LRANGE  -  1,35, 

Send 

d.  Page  4 

Figure  G-14.  Continued. 
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C++++++++++++++++++++++++++++++++ ++■♦■+++++++++++++ ++++++++++++ 
C  CAVITY  C  -  ADDED  HOLE 

C++++++++++++++++++++++++++++ ++++++++++++++++++++ ++++++++++++ 


HOLE  CHEAT lOH  BOOMDARY  DEFINITION 


$BOUNDARy  NAME  ■>  'AFT  CAVITY  HOLE  BOUNDARY', 
MODE  »  'ADD', 

HHOLEIH  X  'STING', 

ISPARTOF  -  'CAVITY  C', 

$END 

$ SURFACE  ISPARTOF  -  'AFT  CAVITY  HOLE  BOUNDARY', 
ORANGE  X  26,26, 

KRANGE  -  5,30, 

LRANGE  X  1,56, 

HVOUT  =  '  -  J'  , 

$END 

^SURFACE  ISPARTOF  "  'AFT  CAVITY  HOLE  BOUNDARY', 
ORANGE  X  26,30, 

KRANGE  X  S,30, 

LRANGE  -  56,56, 

NVOUT  -  '+L', 

$BHD 


C  END  OF  CAVITY  WITH  STORE /STING  CASE 

e.  PaiteS 

Figure  G-14.  Concluded. 
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EXTERIOR  B’  Mesh 


EXTERIOR  A’  Mesh 


Overlap  Region 


’STORE’  Mesh 
Outer  Boundar 


^ ^ 

p - 

-  ■  —  _ 

a.  Overlap  region 
Figure  G>16.  Four*point  overlap. 
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b.  ‘EXTERIOR  A’  mesh 
Figure  G-16.  Continued. 
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EXTERIOR  B’  Mesh 


Overlap  Region 


X  Interpolation 
Boundary  Point  for 
’EXTERIOR  A’  Mesh 


□  Interpolation 
Boundary  Point 
with  Valid  Stencil 

■  Interpolation 
Boundary  Point 
with  Invalid  Stencil 


if'. 


’STORE’  Mesh 
Outer  Boundar 


c.  EXTERIOR  B’  mesh 
Figure  G*16.  Concluded. 
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EXTERIOR  B’  Mesh 


EXTERIOR  A’  Mesh 


Overlap  Region 


Outer  Boundai 


t 

p - 

L  * 

a.  Overlap  region 
Figure  G-17.  Five-point  overlap. 
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b.  ‘EXTERIOR  A*  mesh 
Figure  G-17.  Continued. 
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’EXTERIOR  B’  Mesh 


S) 
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c.  ‘EXTERIOR  B’  mesh 
Figure  G-17.  Concluded. 
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DIAGNOSTIC  HAPS 

LEGEND 

BOUNDARY  POINTS  UPDATED  BY  : 

A  -  EXTERIOR  A 
B  -  EXTERIOR  B 
C  -  CAVITY  A 
D  -  CAVITY  B 
£  -  CAVITY  C 

F  -  STORE 
G  -  STING 
H  -  INTERFACE 


INTERPOLATION  STENCILS  UPDATING  POINTS  IN  : 
a  -  EXTERIOR  A 
b  -  EXTERIOR  B 
c  -  CAVITY  A 
d  -  CAVITY  B 
e  -  CAVITY  C 
f  -  STORE 
g  .  -  STING 
h  -  INTERFACE 


.  -  FIELD  POINT 

I  -  HOLE  POINT 

?  -  ORPHANED  BOUNDARY  POINT 

*  -  INTERPOLATION  STENCIL  UPDATING  MORE  THAN  ONE 

INTERPOLATION  BOUNDARY  POINT 

a.  Page  1 

Figure  G-18.  Sample  diagnostic  maps  for  stoie/sting  in  a  cavity 


HAP  FOR  EXTERIOR  A 


PLANE  L  -  26 


LA 


J123456789 

•12345678901234567890123456789012345678901234567090123456789012345678901234567890123456789012345678 

K 

1  . . . 

2  . bb.BB 

3  ’ . bb.BB 

4  bb.BB 

5  I . ...bb.BB 

$  !! . bb.BB 

7  . bb.BB 

8  !!!! . bb.BB 

9  . bb.BB 

10  . bb.BB 

11  . bb.BB 

12  . .......bb.BB 

13  . hh.hhhh*b.BB 

14  . lih.hhhh*b.BB 

15  . hh . bb.BB 

16  . ** . bb.BB 

17  . hh . bb.BB 

18  . ** . bb.BB 

19  . bb.BB 

20  !!” . ** . bb.BB 

21  . ....** . bb.BB 

22  . *  * . bb.BB 

23  . •• . bb.BB 

24  . hh . bb.BB 

25  . hh . bb.BB 

26  . 

b.  Page  2 

Figiire  G-18.  ConUnned. 
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MAP  FOR  EXTERIOR  B 

PLANE  L  =  26 


J123456789 

123456789012345678901234567890123456789012345e7B901234567B901234567890123456789ai234567a9012345«78 

K 

1  . 

2  AA.aa . * . . . 

3  AA.aa. . . 

4  AA.aa . . . 

5  AA.aa . 

6  AA.aa . 

7  AA.aa . . . 

8  AA.aa . 

9  AA.aa . . . 

10  AA.aa . 

^  11  AA.aa . . . 

0\  12  AA.aa . . . 

13  AAh**hhhhlihhhlihhhhh*h*******h*hWih.hhhhhhhhhhhhhhhh.hhhhhhh.hhh.h., 

14  AAh**hhhhhhhhhhhhhh*h*******h*hhbh.hhhl>hhhhhlilihhbhh.hhbhhhh.hlih.b. . 

15  AA.aa . 1. 

16  AA.aa. . . . . 

17  AA.aa . 

18  AA.aa . . . 

19  AA.aa . . . . 

20  AA.aa . HHH . 

21  AA.aa . GGGGG . 

22  AA.aa . ecllles . 

23  AA.aa . SGlftlfce . 

24  AA.aa... . HG«#««tl«GG . 

25  AA.aa . HOfllfIlfGG . 

26  . 

c.  Page  3 

Figure  G-18.  Continued. 
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HAF  FOR  EXTERIOR  B 


PLANE  L  =  30 


s 


J123456709 

123456789O123456789O12345«709O123456789O123456789O1234567e9O1234567e9O123456789O12345€789O1234567B 

K 

1  . 

2  . . 

3  . . 

4  . . 

5  . . 

6  AA.aa . . . 

7  AA.aa . 

8  AA.aa . 

9  AA.aa . 

10  AA.aa . 

11  AA.aa . 

12  AA.aa . 

13  AA********»»»****»******** ********.**********••*••«.*******. ***.*.. 

14  aa********************************.****************'*******'***-*" 

15  AA.aa . 


17  AA. aaHHHHHHHHHHHHHHHUHHHHHHHJlHHHHHKHHHHHHHHUHHUHHKHHKHHHHHHHHHHHHH. 

18  AA . jiHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH . 

19  AA.HHIt###ll##»####t######f*#######»##l##tl»##l####i#l»######*###». 

20  AA.HH############«#t»i####4»l#######»*####»##*###t####ft#i########. 

21  AA.HH«lt«ti4«t«»f««*«ff«*«ftt«*i«ltfffl<«ltl«l*#»t*f«f«l**«««t««»l. 

22  AA.HHttfff<lllfff«tt«ffff4ftllf«f*t#l««<«l*#fl»*t»ilftt<ft««it«*»». 

23  AA.HH«t««#«tl*tffl»t4»*t««llt#fl«*lf«*l«*#*4««**»li«f««tt»*«**«f*t. 

24  AA.HHII«fl««t«ft««t*«»**«#tt*tf*«ttf«f«f»lf4fl**llifttl*f«ll««tft*. 

25  AA.HH«ftf««ttf»flilffllll«ltf*»«ll*f*«t«««i««f«*ttl#ttt«tt»**l4flt. 

26  . 


d.  Page  4 

Figure  G-18.  Continued. 
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MAP  FOR  CAVITY  A 


PLANE  L  -  56 


J  1  2.3  4  5  6  7  8  9 

12345676901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678 

Httifffttftffifttfffifllltf** 

fffttftftttflltfllfffIfKtIlfl 

####»*a*a*a*a*a*a*a*a*a*a*#i»# 

###**a*a*a*a*a*a*a*a*a*a***### 


««4aa . dd.DD 

«<a*f . dd.DD 

««#•« . dd.DD 

«4f*« . dd.DD 

#«##f . dd.DD 

#«i*f . dd.DD 

#«•*• . dd.DD 

#««#* . dd.DD 

»*•** . dd.DD 

#ftff . dd.DD 

#ll*« . dd.DD 

<««i« . dd.DD 

#!«*# . dd.DD 

. . . . dd.DD 

##lf  f . ffffffffFFFFFF 

»»#•#, ....... .FFFFFFFFFFFFFFPF 

##*#* . FFF#####»###H##i 

*f*f f . ..... .FFfIfItIf f««4t««#« 

*f»*» . FFftIfIftffllfIttif 

tiftt . FFt«fl«t*t##H«f«ftf 


t«f*l - FFt«ff*« 

fi««ff...FF«i<«l< 
!«»««... FF4*t«f#« 
!««*«... FFtt«<««< 
f«#a**.FF|»tfff#4 
f«fl« . 
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Figure  G-18.  Continued. 
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MAP  FOR  CAVITY  B 


PliAHE  L  -  56 


K 

1  f»»««ff«f«««#ff#t«fl««tf«»4#l< 

2  I»f44tf«l#tt<fl*t«*l«lff*tff*4 

3  »<«t#ff«tftf«fft*ff*«l«f#ff«f«< 

4  a*a*a*a»a*a*a*a************b*# 

5  #*a***a*a*a*a*a*******»****b## 


6  CC.cc . .....ee.EB 

7  CC.cc.. . ea.EE 

8  CC.cc . ee.EB 

9  CC.cc . ea.EE 

10  CC.cc . ee.EE 

11  CC.cc . ee.EE 

12  CC.cc . ee.EE 

g  13  CC.cc . ee.EE 

S  14  CC.cc . ee.EE 

15  CC.cc . ee.EE 

16  CC.cc . ee.EE 

17  CC.cc . ee.EE 

18  CC.cc . ee.EE 

19  CC.cc . ee.EE 

20  rFFFFFFFFFFFB'EEFFrFrFFFFFFFFPE 

2 1  FFFFFFFFFFFFFFFFFFFFFFFFFFFFFF 

22  f«««t«fffftt««f«l««fl##«»«#«#F 


23  llffiffff 

24  *«««t«ff» 

25  t««M«tt# 

26 

27  *f««t*«t# 

28  llfflltl# 

29  Iffftttf# 

30  llfffttff 

31  . 
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Figure  G-18.  Continued, 
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HAP  FOR  CAVITY  C 

PLANE  L  -  56 


J1  2  3  4  5  6  7  8  9 

12345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678 
K 

1  *«#«f«»«l*f*f««*»«*»«#*fff*f«« 

2  fffft««t«tf*»t«tt*»*lfi*««l«#« 

3  ffiillff#**#**###****###***#** 

4  4««*.***«***********««**b#|«#« 


5  ###**********************b*### 

6  DD.dd . bbfff 

7  DD.dd . *#### 

8  DD.dd . *ft«f 

9  DD.dd . *t#«f 

10  DD.dd . ««4*f 

11  DD.dd . *!«#« 

12  DD.dd . *«fll 

13  DD.dd . «l*tt 

14  DD.dd . *ftlt 

15  DD.dd. .*****g . ««t«t 

16  DD.dd . 

17  DD.dd . «*#«t 

18  DD.dd . *l<l« 

19  DD.dd . *##ft 

20  FFFFd . fff«« 

21  FFFFF . GGGG . *ffl# 

22  fftfFFFFGGGGGGGGGG . If  til 

23  ItfllFFFGGGGGltflGG . If  Ilf 

24  ffflfflftfflllttllGG . *##4f 


25  «i«tflf«f#|ft«||t«lGG....I«#i4 

26  IlfIfilfItlllillllIGG - *1111 

27  #littf#4flf###«t«tfGG....I«ftl 

28  #lil«f4«4l««#l#«»f»GG....f«### 

29  «#«l*f#«#*tf###ttf#fGG...*f«ff 

30  #llllf##4«4f«t«#tt«fGG...bblf# 

31  . If  If  I 
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Figure  G-18.  Continued. 
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MAP  FOR  CAVITY  C 


FLAME  L  =  80 

J1234587B9 

12345«7e9O123456789O12345€789O1234S6789O123456709O123456789O1234S67B9O123456709O12345e7B9OX234567B 

K 

1  bbbbbbbbbbbbbbbbbbbbbbbbbbbbbb 

2  bbbbbbbbbbbbbbbbbbbbbbbbbbbbbb 


3  OD.dd . BB 

4  OD.dd . BB 

5  OD.dd. . . BB 

6  OD.dd . BB 

7  OD.dd . BB 

8  OD.dd . BB 

9  OD.dd . BB 

10  DD.dd . BB 

11  DD.dd . BB 

12  DD.dd . . . BB 

13  DD.dd . BB 

14  DD.dd . BB 

15  DD.dd . . BB 

16  DDf ***f *******9 . *9 . BB 

17  DD.dd . * . gg . bb 

18  DD.dd..*.  .** . BB 

19  DD.dd..* . BB 

20  DD.dd.  .*.......  .GG . BB 

21  DD.dd..* . GGGG . ...BB 

22  DD.dd..* . GGItGG . BB 

23  DD.dd..* . GGtt«fGG.......BB 

24  DD.dd..* . GG#####QG . BB 

25  DD.dd..* . G6<#t#tGG . BB 

26  OD.dd.. *....GGIftl«tGG - g.BB 

27  DD.dd..* _ GGl««ttlGG. .....BB 

28  DD.dd . GGff«f##tGC . BB 

29  DD.dd..* - GGItl*###GG . BB 

30  DD.dd . GGffItfffGG . BB 
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Figure  G-18.  Continued. 
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HAP  FOR  STORE 


PLANE  L  -  29 


K> 
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CCCCCCCCCOCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCDDDDDDDDDDDODDDDDDDDDDDDDDODDEEEEEEEEEEEE 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCDDDDDDDDDDDDDDDDDDDDDDDDDDODDEEEEEEEEEEEE 

CCCCCCCCCCCCCCCCCCCCCCCCCCOCCCCCCCOCCCCCCDDDDDDDDDDDDDDDDDDDDDDDDDDDDDEEEEEEEEEEEE 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCDODDDDDDDDDDDODDDDDODDDDDDDDDEEEEEEEEEEEE 

CCCCCCCCCCCCCCCCCCCOCCCCCCCCCCCCCCCOCCCCCDDDDDDDDDDDDDDDDDODDDDODDDDDDEEEEEEBEEEEE 

CCCCCCCCCCCCCCCCCCCCCCCCCCCOCCCCCCCCCCCCCDDDDDDDDDDDDDDDDODDDDDDDDDDDDEEEEEEEEEEEE 

CCCCCOCCCCOCCCCCCCCCCCCCCCCCOCCCCCCCCCCCCDDDDDDDDDDDDDDDDDDODDDDDDDDDDEEEEEEEEEEEE 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCOCCCCCCCCCCCCDDDDDDDDDDDDDDDODDDDDDDDDDDDDEEEEEEEEEEEE 


K 

1 

2 

3 

4 

5 
8 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


21  CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCDDDDDDDDDDDDDDDDDDDDDDDDDDDDDEEEEEEEEEEEE 

22  CCCCCCX:CCCCCCCCCCCCCOCCCCCCCCCCCCCXX:CXX:CCDDDDDDDDDDDDDDDDODDDDDDDDDDDDEEEEEEEEEEEE 

23  CCCCCCOCCCCCCCCCCCCCOCCCCCCCCCCCCCCCCCCCCDDDDODDDDDDDDDDDDODDDDDDDDDDDEEEEEEEEEEEE 

24  CCCCCCCCCCCCCCCCCCCCCOCCCCCCCCCCCCOXCCCCDDDDDDDDDDDDDDDDDDDDDDDDDDDDDEEEEEEEEEEEE 
CCCCCCCOCCCCOCCCCCCXXOCCCCCCCCCCCCCCCCOCCDDDDOODDDDDDDDDDDODDDDDDDDDDDEEEEEEEEEEEE 


25 

26 

27 

28 

29 

30 

31 

32 

33 


CCCCCCCCCCCCCCCAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABBBBBBBBBBBBBBBBBBBBBBB 

CCCCCCCCCCCCCCAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABBBBBBBBBBBBBBBBBBBBBBB 

ccccccccocccccaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaabbbbbbbbbbbbbbbbbbbbbbb 

CCCCCCCCCCCCCAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABBBBBBBBBBBBBBBBBBBBBBB 

CCCCCCCCCCCCCAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABBBBBBBBBBBBBBBBBBBBBBB 

CCCCCCCCCOCCCAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABBBBBBBBBBBBBBBBBBBBBBB 
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Figure  G>]8.  Continued 
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MAP  FOR  STING 


PLANE  L 


29 


J123456789 

12345678901234567e901234567S9012345fi7e9012345678S012345e7e9012345S78»012345678901234Sfi7e»01234S67S 


KJ 

IM 


K 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 
13 
Ifi 
17 
la 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 


EEEEEBEEEEBEEEEEEEIEBBBBBBBBHHHHHHHHHHBHHHHHHHHHHHHHIIHHHHHIIHHHHHHHHHKHHHHHHHHHHMH. 

EEEEEEEEEEEEEEEEEBIEEEEBBBBBHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH. 

EEEEEEEEEEEEEEEEEEEEEEEBBBBBHHHHHHHHHHHHHlUtHHHHHRimHHHHHHHHHIUUIHHHHHHHHHHHHHHHHHa. 

T;i;BERRi;EKRKEF.KF.RRRgEF.EEBBBaBHHHHHHHHHlUmHHHHHHIltlHllHHHHHIlllHHHHHHHHHHHHHHHHHHHHHHHH. 

RRREEEEEEEBEEEF-RgggEEEEBBBaBHHHHHHHHHHRHHHHHHHHHRHHHHHHHHHHHHHHHHWHHHHHHHHHHHHHHa. 

EEEEEBEEEBEEEEEEEEEEEEEBBBBBHIUlHHHHHHHHHHHKHHlIHIlIlHBililHHHHRHHHHIUUiHHHiUIHHHMHHHHHHH. 

EEEEEEEEEF.RREP.EEREREEEEHBBBBHHHHHHHHHHHHHHHHHHHIlHHHHHHHHIlHHHHmtmHHHHHHIUlHHlIHHaHH. 

REF.REEEEEEEEEEEEEEEEEEEBBBBBHHHHHHHHHHHHHHimHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH. 

EEEEEEEEEEEEEEEEEBEEEEBBBBBBBHHHHHHHHHHHHHHHHHHHHHBBHHHHHHHHHHHHHHHRHHRHHHHHHHHHH. 

EEEEEBBEEEEB^EEEEEEEEB^BB^HHHHHHHHHHHMHHH^HHHHIWHiraHHHHHHHHHHHirailHHHHIUlHHH. 

EEEEEEEEEEEEEEEEEEEEEBBBBBBBBHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH. 

eeeeeebeeeeeeeeeeeebebbbbbbbbbhhhkhhhhhhhhhhhhhhhhhhhhhhbhhhhhhhbhhhhhhhhhhhhhbhh. 


REEREEEEEEEEEEEEEEEEBBBBBBBBBBHHHHHimHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHIfflHHHHHHHHHH. 


EEEEEBBEEEEEEEEEEBBBBBBBBBBBBBBHHHHHHHHHHHHHHHHHHHHIUltlHHHHHlIHHHHHHHHHHHHHHHHHilHHH. 


BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBHHBHHHHHHHHHHHHHHHHHaHRHHHHHHHHHKHHHHHHHHHHHHHHHH. 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH. 

bbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbhhhhhhhhhhhhhhhhhhhhhhhhhhhuhhhuhhhhhhhhhhhhhhhh. 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHKHHHHHHHHHHHH. 

rbbbrrrrrrrrbhrrrrhbbbbbrrbrrhhrhhhhhhhhhhhhbhhhbhhhhhhhhhhiuihhhhhhhhhhhhhhhhhhhh. 


BBBBBBBBBBBBBBBBBB8B6BBBBBBBBBBBBBBHHHHHHIUUIHHHHHHHHHIUU!HHHaHIlHHHHHHHHHKHHHHHHHHH. 
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Figure  G-18.  Continued. 
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MAP  FOR  INTERFACE 

PLANE  L  =  1 


J1234S6789 
1234567890 123456789O123456789O12345S709O123456789O123456789O1234567B9O123456789O123456789O12345678 

^1  DPPpprinnnminnnnf>BBBRREF.RBREEEEEREEEEEEEEEEEEKEEEKBBBBBBBBBBBBBBBBBBBBBBBBBBB. 

2  dddddddddoddddddbeeeeeeeeeeeeeeeeeeebeeeeeebebebebbbbbbbbbbbbbbbbbbbbbbbbbbb. 

3  dddddddddddddoddeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeebbbbbbbbbbbbbbbbbbbbbbbbbbb. 

4  ddddddddddddodddeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeebbbbbbbbbbbbbbbbbbbbbbbbbbb. 

5  PnnT»Pr>PPPPPPpr>I>»Pg^>g^-’^'^EFggBF.RRRREEEEEEEEEEEEEBEEBBBBBBBBBBBBBBBBBBBBBBBBBBB. 

6  dddddddddddddoddeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeebbbsbbbbbbbbbbbbbbbbbbbbbbb. 

7  DDDDDDDDDDDDOOODEEBEEEEEBEEEEEEEEEEEEEEEEEEEEEBEEBBBBBBBBBBBBBBBBBBBBBBBBBBB. 

8  PPPPPPpnnru^nPPPPEBBP-BHEEgggEEgggEBBEEEEEEEEEEEEEEBBBBBBBBBBBBBBBBBBBBBBBBBBB. 

9  dDDDDDDDDDDDDDDDEEEEEE£CEEEEEBGG<3<3GOGCGCGGEEBEEEBBBBBBBBBBBBBBBBBBBBBBBBBBBB  . 

10  FFFF  PFFFFFPFFFFFEEEEEEGGGEBGGeGGGGGGGGGGGGGGGEEEEBBBBBBBBBBBBBBBBBBBBBBBBBBB . 

11  FFFFFFFFFFFFFFFFEBEEBEGGGGGGGGGGOGGGGGGGOGGGGGEEEBBBBBBBBBBBBBBBBBBBBBBBBBBB. 

12  FFFFFFFFFFFFFFFFEEEEEEGGGGGGGGGGOGGGGGGGGGGGGGGEBBBBBBBBBBBBBBBBBBBBBBBBBBBB* 

13  FFFFFFFFFPFFFFFFEBEEEEGCGGGeeCGGGGGGGGGGGGGGOGGGBBBBBBBBBBBBBBBBBBBBBBBBBBBB . 

14  ffffffffffffffffeeeeeeggggggggggggggggggggggoggggbbbbbbbbbbbbbbbbbbbbbbbbbbb . 

15  FFFFFFFFFFFFFFFFEEEEEEGGGGGGGGGGG|###i###GGGGGGGGBBBBBBBBBBBBBBBBBBBBBBBBBBB . 

1 6  FFFFFFFFFFFFFFFFEEEEEEGGGGGGGGGittf f f Ilf ItGGGGGGGBBBBBBBBBBBBBBBBBBBBBBBBBBB  > 

17  FFFFFFFFFFFFFFFFEEEEEEGGGGGGGGGIIttllllf IIIGGCCGGBBBBBBBBBBBBBBBBBBBBBBBBBBB . 

18  FFFFFFFFFFFfFFFFEEEEEEGGGGGGGGIIIfflllilltfCGGGGGBBBBBBBBBBBBBBBBBBBBBBBBBBB. 

1 9  FFFFFFFFFFFFFFFFEEEEEEGGGGGGG6I Ittlittlitt f GGGGGGBBBBBBBBBBBBBBBBBBBBBBBBBBB . 

20  ffffffffffFFFFFFEEEEEEGGgggggcI Itllflllllfl igcgggbbbbbbbbbbbbbbbbbbbbbbbbbBB . 

21  . 
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Figure  G-18.  Continued. 
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HftP  FOR  INTERFACE 


PLANE  L  o  10 


to 

<oi 


J 


123456789 

1234567B9012345678901234S6799012345678901234567B9012345678901234567B90123456789012345678901234S67e 


K 

1 

2 

3 

4 

5 

6 
7 

e 

9 

10 

11 

12 

13 

14 

15 

16 
17 
IB 

19 

20 
21 


AAAAAAABBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB. 

AAAAAAABBBBBBftBBRBnRBaBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB. 


_ *a****a.  .  .  . . 

.  .  .  .  *  *  *<T  ..  .  _ ........... 

**_  _  _ .  .  .fl*** . 

*** 

. CCCGCe.  *  * . *** . 

* 

.  . .  .GGGGGGGGGG . 4 ..  . 

&& . 

.  .  .GGGfIlifIGGGG . g*«»»n*gg.gg.g 

.g-G- -gg-gggg- 

AA . 

..GGlifItillliGGGG . 

AA . 

*****  ***•****. 

AA . 

.Ge4llftt<*li4t«GGGG.. . . 

.  *.*.gg*.****. 

.GG4#tlf##ttf IflGGGG . 

AA . 

GGffittiftftiiHGGGGG . 

AA . 

. *<3- 

GG#iill»»H#l#»#GGGGG. 

.*.*,*•*•****• 
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Figure  G-18.  Continued 
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MAP  FOR  INTERFACE 

PLANE  L  -  15 

'^12345€789O123456709O123456789O12345€189O1234567B9O12345«789O1234567«9O1234567B9O123456709O12345678 


3  AA. 

4  AA 

5  AA 

8  AA 
7  AA 
0  AA 

9  AA 
10  AA 

. GCGCGG . G-G.  .G . 

. GGGGGGGGGGGGeeCGGGGGGGOGGGGGGGCGG . 

, . .GGGGGItlllfGGGGGGGGGGGGGGtGfGGfGGOG. 

. . . .GGGGilf*«ltl*f#«tf«t«l«*«#»tflt»f««t** 

. .GGGG*«»«f*lt»ff«tffttf 

GGGGG####ll»#»#####«ii«#t*##«#*»*««*«*« 

21  . * . 


11  AA. 

12  AA. 

13  AA. 

14  AA. 

15  AA. 
.16  AA. 

17  AA, 
10  AA. 

19  AA 

20  AA 


. g.,g. ***.*** 

. g...g..g . 

•  9 - g . 

gg . 
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Figufe  G-18.  Concluded. 


Flow 
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■.  Graphical  depiction 
Figure  G-19.  'STORE  CAP’  mesh. 
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AEDC-TR-91-8 


AEDC-TR-91-8 


d.  Front  view 
Figure  G-19.  Concluded. 
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**♦*****•******«•»****+*****************************■********* 
PEGSUS  4.0  INPUT 

EXAHPI>E  6-3  STORE  ALONE  IN  A  CAVITY 
7  HESH  CASE 

STING  AND  INTEBEACE  MESHES  SUBTRACTED. 

STORE  CAP  MESH  ADDED. 

*********ik*****1»  »*****»•  ************************************ 


**-****  ^Al^**-*  ***-**-*  *-*««*********  *****l^*4r1^)^**1^**1k*****1^**Alfr^1lt^ 


MESH  DEFINITIONS 

A*******  ********************  ******  ****************** 


$MESH  NAME  -  'STING', 
MODE  =  'SUB', 


$MESH  NAME 

MODE 


' INTERFACE'  , 
=  'SUB', 


SEND 


Send 


SMESH  NAME  >=  '  STORE  CAP'  , 

MODE  =  'ADD', 

LINK  -  ' STORE' , 

'CAVITY  C', 

'EXTERIOR  B', 

'CAVITY  B', 

'EXTERIOR  A', 

' CAVITY  A'  , 

ADDLINK  =»  'STORE', 

'EXTERIOR  A', 
'EXTERIOR  B', 

'CAVITY  A', 

‘  'CAVITY  B', 

'CAVITY  C', 

xO  =  15.468,  YO  =  0.0,  ZO  -  0.0, 

SEND 


a.  Page  1 

Figure  G-20.  User  input  for  store  alone  in  a  cavity. 
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************************************************************* 

INTERPOLATION  BOUNDARY  DEFINITIONS 
************************************************************** 


OUTER  BOUNDARY  DEFINITION 


^BOUNDARY  NAME  =  ' STORE  CAP  OUTER  BOUNDARY' , 
ISPARTOF  =  'STORE  CAP', 

$END 

SSURFRCE  ISPARTOF  =  '  STOEtE  CAP  OUTER  BOUNDARY', 
JRANGE  =  12,39, 

KRANGE  ~  2,32, 

LRANGE  “  26,27, 

SEND 

SSURFACE  ISPARTOF  -  'STORE  CAP  OUTER  BOUNDARY', 
ORANGE  =  38,39, 

KRANGE  -  2,32, 

LRANGE  -  1,27, 

SEND 

SSURFACE  ISPARTOF  -  'STORE  CAP  OUTER  BOUNDARY', 
ORANGE  «  1,39, 

KRANGE  >  2,32, 

LRANGE  =  1,2, 

SEND 

INTERIOR  REGION  DEFINITION 


SREGION  NAME  -  'STORE  CAP  REGION', 

TYPE  =  'IHTR', 

ISPARTOF  -  'STORE  CAP', 

SEND 

SVOLUME  ISPARTOF  -  ' STORE  CAP  REGION' , 
JRANGE  -  1,11, 

KRANGE  =  1,33, 

LRANGE  -  26,27, 

SEND 

b.  Page  2 

Figure  G-20.  Continued. 
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C++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
C  STORE 

C++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
C - 

C  OUTER  BOUNDARY  DEFINITION 

C - 


$BOUNDARY  NAME  -  'AFT  STORE  OUTER  BOUNDARY', 
MODE  =  'ADD', 

ISPARTOF  =  'STORE', 

$EHD 


$ SURFACE  ISPARTOF  »  'AFT  STORE  OUTER  BOUNDARY', 
JRAHGE  -  77,62, 

KRAHGE  2,32, 

LRANGE  -  1,2, 

$END 


C+++++++++++++4++++++++++++++++++++++++++++++++++++++++++++++ 
C  CAVITY  C 

C++++++  4+  ++ ++++++++++++  4+  ++  ++++++++++++++++++++++++++++++++++ 

C - 

C  SUBTRACTED  OUTER  BOUNDARY  DEFINITION 

C - 


Sboundary  name 

MODE 


'AFT  CAVITY  HOLE  BOUNDARY', 
'SUB'  , 

$END 


C  END  OF  STORE  ALONE  IN  A  CAVITY  CASE 


c.  Page  3 

Figure  G-IO.  Concluded. 
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DIAGNOSTIC  HAPS 

LEGEND 

BOUNDARY  POINTS  UPDATED  BY  : 

A  -  EXTERIOR  A 
B  -  EXTERIOR  B 
C  -  CAVITY  A 
D  -  CAVITY  B 
E  -  CAVITY  C 
F  -  STORE 
G  -  STORE  CAP 


INTERPOLATION  STENCILS  UPDATING  POINTS  IN  ; 
a  -  EXTERIOR  A 
b  -  EXTERIOR  B 
c  -  CAVITY  A 
d  -  CAVITY  B 
«  -  CAVITY  C 

£  -  STORE 

g  -  STORE  CAP 


.  -  FIELD  POINT 

#  -  HOLE  POINT 

?  -  ORPHANED  BOUNDARY  POINT 

*  -  INTERPOLATION  STENCIL  UPDATING  MORE  THAN  ONE 

INTERPOLATION  BOUNDARY  POINT 

a.  Page  1 

Figure  G‘22.  Sample  diagnostic  maps  for  store  alone  in  a  cavity. 
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HAP  FOR  CAVITY  C 


PLANE  L  »  56 


J1234S6789 

12345678901234567690123456789012345678901234567890123456789012345670901234567890123456789012345678 

K 

2  llttt«*«*lt*«4ttf«l*««#«#ftllt 

3 


4  #|#*********************b#### 

5  . *.**b*«« 

6  DD.dd . . . .  , . bbtfl 

7  DD.dd . *•## 

8  DD.dd . *#ffi 

9  DD.dd . *l»#' 

10  DD.dd . Ititi 

11  DD.dd . •###( 

12  DD.dd . **!*< 

13  DD.dd . #♦##! 

tin  14  DD.dd . *ttfi 

15  DD.dd . titf 

16  DD.dd . *ttf 

17  DD.dd . ftti 

18  DD.dd . *### 

19  DD.dd . *t»# 

20  FFFFd . t### 

21  FFFFF . *•## 

22  ti#«FFFFG . ttfl 

23  tfIffFFFGG . ttll 

24  »»««««*ttGG . •### 

25  f««ft4ltlGG . *tf» 

26  «««»#|t|lGG . *#»» 

27  f«»#4f#tlGC . i - «f«» 

28  fffllfIttGG..*..., . *#*» 

29  fifHlilIGG . *f#f 

30  *«#«««t*»GG..* . bb## 

31  . #»»» 


b.  P»ge  2 

Figure  G-22.  Cktntinued. 
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MhP  FOR  STORE  CAP 


FLAME  L  -  27 


J  1  2  3  4  5  6 

12345678901234567890123456769012345678901234567690123456789012341 


10  ft 

11  #« 


20  t# 

21  It 

22  It 

23  II 

24  ft 

25  It 

26  II 

27  I# 

28  It 

29  It 

30  If 

31  It 

32  II 

33  ft 


I  #  FFFFFFFFFFFFFFFFFFFFFFFPFEEE 
1 1 ppppFFf pfffppppppffFFFFFFEEE 
IfFFFFFFFFFFFFFFFFFFFFFFFFFEEE 
IIFFFFFFFFFFFFFFFFFPPFFFFFFEEE 
IIFFFFFFFFFFFFFFFFFFFFFFFFFEEE 
flFFFFFFFFFFFFPPPFFFFFFFFFFEEE 
IIFFFFFFFFFFFFFFFFFFFFFFFFFEEE 
IIFFFFFFFFFFFFFFFFFFFFFFFFFEEE 
I IFFFFFFFFFFFFFFFFFFFFFFFFFEKE 
IIFFFFFFFFFFFFFFFFFFFFFFFFFEEE 
IIPFFFPFFFFFFFFFFFFFFFFFFFFEEE 
1 1 FFFFFFFFFFFFFFFFFFFFFFFFFEEE 
IIPFFFFFFFFFFFFFFFFFFFFFFFFEBE 
I IFFFFFFFFFFFFFFFFFFFFFFFFFEEE 
IIFFFFFFFFFFFFFFFFFFFFFFFFFEEE 
I I FFFFFFFFFFFFFFFFFFFFFFFFFEEE 
I  IFFFFFFFFFFFFFFFFFFFFFFFFFEEE 
I  IFFFFFFFFFFFFFFFFFFFFFFFFFEEE 
I  IFFFFFFFFFFFFFFFFFFFFFFFFFEEE 
1 1 FPPPFFFFFFFFFFFFFFFFFFFFFEEE 
I IFFFFFFFFFFFFFFFFFFFFFFFFFEEE 
I  IFFFFFFFFFFFFFFFFFFFFFFFFFEEE 

1 1 ffffffffffffffffffff  fffffeee 

I IFFFFFFF  FFFFFFFFFFFFFFFFFFEEE 

I ifffffffffffffffffffffffffebb 
I  IfffffffFFFFFFFFFFFFFFFFFFBBB 
1 1 ffffffffffffffffffff  fffffbbb 
I J fffffffffffffffffffffffffbbb 
1 1 ffffffffffffffffffff  fffffbbb 
1 1 fffffffffffffffffffffffffbbb 
I  ifffffffffffffffffffffffffebb 
•I.. . . 
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Figure  G-22.  Conduded 
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Table  G-1.  Mesh  Dimensions 


Mesh  Name 

JMAX 

KMAX 

LMAX 

EXTERIOR  A 

67 

26 

35 

EXTERIOR  B 

67 

26 

35 

CAVITY  A 

30 

31 

85 

CAVITY  B 

30 

31 

85 

CAVITY  C 

30 

31 

85 

STORE 

82 

33 

29 

STING 

82 

33 

29 

INTERFACE 

77 

21 

35 

STORE  CAP 

39 

33 

27 
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APPENDIX  H 

SAMPLE  CRAY  JOB  CONTROL  LANGUAGE 

The  Cray  job  control  language  for  each  example  given  in  Appendix  G  is  shown  in  Figs. 
H-1,  H-2,  and  H-3,  respectively. 


228 


AE0C-TR-gi-8 


JOB,  JN-B30276A,  T-900,  HFL. 

ACCOUNT, AC-58000502, US-B30276 . 

FETCH,  DN-PEGSUS,DF-CB, 

TEXT-' DSN-B30276 . PEGSDS4C . FOKT, DISP-SHR' . 

FETCH, DN-GRIDl , DF-TR, * 

TEXT-' DSN-B30276. CAVITY. GE0M16(CAVV1G1) , DISP-SHR' . 
FETCH, DN-GRID2 , DF-TR, " 

TEXT-'DSH-B30276.CAVITY.GEOM16(CAVVLG2) , DISP-SHR'  . 

pcTfif  nvnTD  ^ 

T£XT^'DSH-B30276.CAV1TY.CEOM16(CAVVLG3) , DISP-SHR' . 
FETCH,  ON-GRID4,  DF^TR, '' 

TEXT-' DSN-B30276. CAVITY. GEOM16(CAVVLG4) , DISP-SHR'  . 
FETCH,  DN-GRID5,  DF-TR, 

TEXT-' DSN-B30276. CAVITY. GEOM16  (CAWLG5)  ,  DISP-SHR'  . 


COP  YR, I -GRI D 1 , 0- IHSRI D , NR-3 . 
COE  YR, I =GRID2 , 0= INGRI D ,  NR-3 . 
COPYR, I-GRID3, 0-INGRID, NR-3. 
COPYR, I-GRID4, 0-IHGRID, NR-3. 
COPYR, I-GRID5, 0-INGRID, NR-3. 
REWIND, DN- INGRID . 


CFT77, I-PEGSUS.ON-F. 
L0R,SET-INDEF. 


DISPOSE,  DN-MAPS,DC-ST,  DF-CB, '' 

TEXT-'DSH-B3027e.PECOUT.DATA{MAP5TRl) ,DISF-SHR'  . 

DISPOSE ,  DN-ORPHAH,  DC-ST,  DF-CB, '' 

TEXT-'DSH-B3027e.PEGOUT.DATA{ORPHTRl) ,DISP-5HR'  . 

DISPOSE, DH-COMPODT,DC-ST.  DF-TR,  " 

TEXT-' DSN-B30276. PEGSDS40. GEOM(COMPTRl) , DISP-SHR' . 

DISPOSE, DN-INTODT, DC-ST, DP-TR,  " 

TEXT-' DSN-B30276. PEGSUS40. GEOM(INTRTRl) , DISP-SHR' . 

DISPOSE, DN-RSTOUT, DC-ST, DP-TR, ^ 

TEXT-'DSN-B30276.PEGSUS40.GEOH(RSTTR1> ,0ISP-SHR' . 

DISPOSE,  DN-IBPLOT,  DC-ST,  DP-TR,  '' 

TEXT-' DSN-B30276. PEGSUS40. GEOM(IBLKTfil) , DISP-SHR'  . 


EXIT. 

Figure  H-l.  Cray  Job  control  language  for  empty  cavity. 
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JOB , JH-B3  0276A,T=900,MFL. 

ACCOUNT , AC-58000502, as=B30  276. 

FBTCH,  DH-PEGSUS ,  DF-CB, 

TEXT- ' D8N-B3  0276. PEGSUS4  0 . FORT , D ISP-SHR' . 

FETCH , DN-GRID 1 , DF-TR,  * 

TEXT-' DSH-B30276. CAVITY. GEOH13(STOREPHa) ,DISP=SHR’ . 
FETCH,  DN-GRID2 ,  DF=^TR. '' 

TEXT-' DSH-B30276.CAVITY.GEOK10  <STINGPH6) , DISP-SKR' . 
FETCH ,  ON-GRI D  3 ,  DF=TR, 

TEXT-' DSN-B30276. CAVITY. GEOK16 (CAWLIF)  ,DISP-SHR'  . 

» 

FETCH,  DN-COMPIM,  DF-TR, '' 

TEXT-'DSN-B30276.PEGSUS40.GEOM(COMPTR1) ,DISP-SHR'  . 
FETCH,  DN-INTIN,  DF-TR, 

TEXT-'DSN-B30276.PEG3DS40.GEOM(1HTHTR1) ,DISP-SHR' . 
FETCH, DN-RSTIN, DF-TR, * 

TEXT-'DS1)-B30276.PEGSOS40  .GEOH(BSTTRl)  ,  DISP-SBR'  . 


COFYR, I-GRID1,0-INGRID,NR=3. 
COPYR,  I-GRI02,0-rNGHID,lJR=3. 
COPYR,  I -GRID3, 0-INGRID, NR-3. 
REMIND, DN-IHGRID . 


CFT77, I-PEGSOSl .ON-F , 
LDR, SET-INDEF. 


DISPOSE, DN-MAPS,DC-ST,DF-CB, '' 

TEXT-' DSH-B30276.FEGOUT.DAXA(HAPSTR2)  ,  OISP-SHR'  . 

DISPOSE, DN=ORPHAN,DC-ST,DF-CB, '' 

TEXT-' DSH-B30276.FEGOUT.DATA(ORPHTR2) ,DISF-SHR'  . 

DISPOSE,  DN-COMPODT,  DC-ST,  DF-TR, '' 

TEXT-'DSN-B30276.PEGSDS40.GEOM(COMPTR2) ,DISP-SHR' . 

DISPOSE, DN-IHTOUT, DC=ST, DF-TR, " 

TEXT-'DSH-B30276.PEGSDS40.GEOM(IHTRTR2} ,DISP-SHR' . 

DISPOSE,  DN-RSTOUT,  DC=ST,  DF-TR, 

TEXT-' DSN-B30276 . PBGSDS4  0 .GEOM (RSTIR2) , DISP-SHR' . 

DISPOSE.  DN- IBP  LOT.  DC-ST,  DF-TR, 

TE]CT-'DSN-B30276.PBGSDS40.GEOM(IBLKTR2)  , DISP-SHR'  . 


EXIT. 

Figure  H-2.  Cray  job  control  language  for  store/sting  added  to  empty  cavity. 
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JOB,  JNi=B30276A,  T-9D0,  MFl. 

ACCOONT ,AC>=5 6000502,  US-B3027 6 . 

FETCH,  DN-PEGSUS,DP=CB, 

TEXT=' DSN-B30276 . PEGSUS40 . FORT, OISP^SHR' . 

FETCH, DB-GRIDl,Dr=TR, * 

T£XI='DSK-B30276.CAVITy.GEOM13{STORECAP) ,DISP°5HR' . 
FETCH,  Dll=COMPIN,  DF*=TR, '' 

TEXT-'DSH->B30276.PEGSDS40.GEOM(COMPTR2)  ,DISP-SHR'  . 
FETCH, DH=INTIN, DP-TR, * 

TEXT»'DSH-B30276.PEGSUS40.GEOH(IMTRTR2) ,D1SP=SHR' . 
FETCH, DN=RSTIN, DP=TR, * 

TEXT>'DSN-B30276.PEGSnS4O.GEOH(RSTTR2) ,DISP«SHR' . 


COPYR, I=GRID1, 0-INGRID, NR=3. 
REWIND, DH- INGRID . 


CFT77, I-PEGSOSl, ON=F. 
LDR.SET-INDEF. 


DISPOSE , DN=HAPS, DC-ST, DF-CB, " 

TEXT-’ DSN-b30276.PE60DT.DATA{HAPSTR3) ,DISF-SHR'  . 

DISPOSE, DH=ORPHAH,  DC-ST, DP-CB, 

TEXT-' D5N-B30276 .PEGODT .DATA (ORPHTR3) , DISP-SHR'  . 

* 

DISPOSE, DN-COMPOUT, DC-ST, DF-TR, “ 

TEXT-' DSHf B30276 . PEGSDS40 . GEOM (COMPTR3)  ,  DISP-SHR' . 

DISPOSE, DM-INTODT, DC-ST,  DF-TR,  *■ 

TEXT- '  p  5H-B3  02  7  e . PEGSUS  4  0 . GEOM ( INTRTR3 ) , DISP-SHR' . 

DISPOSE, DN-RSTOUT, DC-ST, DP-TR,  * 

TEXT-' O5N-B30276.PEGSDS40. GEOM {RSTTR3) , DISP-SHR'  . 

DISPOSE, DN-IBPLOT, DC-ST, DP-TR,  " 

TEXT-’DSN-B30276.PEGSU540.GEOM(IBLXTR3) , DISP-SHR'  . 

* 

* 

EXIT. 

Figure  H-3.  Cray  job  control  language  for  subtracted  sting  and  added  store  cap. 
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